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PREFACE 


The  study  reported  herein  is  a  component  of  the  Potamology  Program  (P-1) 
of  the  Lower  Mississippi  Valley  Division  (LMVD) .  The  Potamology  Program  is 
conducted  under  the  direction  of  the  Commander,  LMVD,  and  is  a  comprehensive 
study  of  physical  forces  that  influence  the  flood-carrying  capacity  and  navi¬ 
gability  of  the  lower  Mississippi  River.  The  purpose  of  the  Potamology  Pro¬ 
gram  is  to  define  cause -and-effect  relationships  that  result  in  short-  and 
long-term  changes  in  the  s tage - discharge  relationships  of  the  lower  Missis¬ 
sippi  River  and  to  develop  improved  design  concepts  and  criteria  for  construc¬ 
tion  of  channel  stabilization  works  that  will  improve  flood  control  and 
navigation  along  this  river. 

The  Potamology  Program  is  composed  of  two  major  components:  Sedimenta¬ 
tion,  Mississippi  River  Basin;  and  Aggradation  and  Degradation,  Mississippi 
River.  This  study  is  one  item  under  the  Sedimentation  component.  Future 
studies  will  be  directed  toward  development  and  utilization  of  a  flow  sediment 
model  capable  of  detailed  investigations  of  short-  and  long-term  sedimentation 
trends  locally  and  throughout  the  main  stem  Mississippi  River. 

The  study  reported  herein  was  the  responsibility  of  the  US  Army  Engineer 
District,  New  Orleans  (LMN) ,  New  Orleans,  LA.  The  LMN  contracted  with  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  MS,  for  con¬ 
duct  of  the  study  and  preparation  of  the  report.  The  study  was  conducted  by 
WES  from  1  June  1982  to  31  May  1983. 

This  investigation  was  Phase  IV  of  a  four-phase  study  of  the  sediment 
regime  of  the  Mississippi  River  Basin.  Phase  I,  conducted  under  the  earlier 
LMVD  Potamology  Program  (T-l),  resulted  in  the  publication  of  WES  Technical 
Report  M-77-1,  "Inventory  of  Sediment  Sample  Collection  Stations  in  the  Mis¬ 
sissippi  River  Basin,"  in  March  1977.  The  end  product  of  Phase  II  was  LMVD 
Potamology  Program  (P-1)  Report  1,  "Characterization  of  the  Suspended- Sediment 
Regime  and  Bed-Material  Gradation  of  the  Mississippi  River  Basin,"  published 
in  August  1981.  The  Phase  II  study  identified  a  downward  trend  in  Mississippi 
River  suspended- sediment  loads  that  apparently  began  around  the  middle  of  the 
20th  century.  Phase  III  dealt  with  suspended- sediment  sampling,  analysis,  and 
load-computation  procedures  used  at  key  stations  on  major  streams  in  the  Mis¬ 
sissippi  River  Basin  and  the  possible  correlation  of  these  procedures  with  the 
downward  trend  in  suspended- sediment  loads.  In  this  Phase  IV  study,  existing 
regime  equations  were  examined  in  light  of  their  applicability  to  major 


alluvial  river  systems.  The  principal  focus  of  the  Phase  IV  study  was  the 
examination  of  the  effects  of  changes  in  suspended-sediment  loads  as  they 
relate  to  the  regime  of  alluvial  rivers. 

Publication  has  been  delayed,  with  the  intention  of  later  including  the 
findings  of  this  study  with  those  of  follow-on  efforts.  However,  the  antici¬ 
pated  funding  required  to  conduct  follow-on  efforts  did  not  become  available, 
and  as  a  result,  this  report  has  remained  unpublished  until  now. 
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CONVERSION  FACTORS ,  NON-bl  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Mul t inlv 

Bv 

To  Obtain 

cubic  feet  per  second 

0.02831685 

cubic  metres  per 
second 

Fahrenheit  degrees 

5/9 

Celsius  degrees 
or  Kelvins* 

feet 

0.3048 

metres 

inches 

25.4 

millimetres 

miles  (US  statute) 

1.609347 

kilometres 

pounds  (mass) 

0.4535924 

kilograms 

pounds  (mass)  per 
cubic  foot 

16.01846 

kilograms  per 
cubic  metre 

tons  (2,000  pounds, 
mass ) 

907.1847 

kilograms 

*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  the  following  formula:  C  =  (5/9) (F  -  32).  To  obtain  Kelvin  (K) 
readings,  use:  K  =  (5/9)(F  -  32)  +  273.15. 
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IMPACTS  OF  CHANGES  IN  SUSPENDED  -  SEDIMENT  LOADS  ON  THE 


REGIME  OF  ALLUVIAL  RIVERS 

PART  I:  INTRODUCTION 

Background 

1.  For  nearly  a  century,  fluvial  geomorphologists  and  river  engineers 
have  attempted  to  unravel  and  define  the  complex  natural  and  human  -  induced 
phenomena  that  have  influenced  the  character  of  streams,  particularly  alluvial 
rivers.  Despite  extensive  research,  serious  gaps  in  understanding  these  phe¬ 
nomena  and  a  lack  of  clear  consensus  among  researchers  regarding  their  quanti¬ 
tative  interrelations  still  remain. 

2.  The  concept  of  a  "graded  condition"  in  streams  can  be  traced  to  the 
late  19th  century  writings  of  G.  K.  Gilbert,  who  stated  that  "equilibrium  of 
action  in  streams  consists  of  a  mutual  adjustment  between  velocity,  discharge, 
slope,  and  load"  (Ritter  1986).  In  the  early  part  of  the  20th  century, 

W.  M.  Davis  (1902)  introduced  the  terms  "grade"  and  "graded  slopes"  to 
describe  the  balanced  fluvial  condition.  The  idea  of  a  "graded  river"  was 
essentially  unchallenged  until  Mackin  (1948)  redefined  the  concept  to  include 
slope  adjustment  "over  a  period  of  years."*  Fluvial  geologists  and  geomor¬ 
phologists  prefer  to  use  the  term  "graded"  or  "poised"  to  describe  such  a 
channel  that  is  in  "quasi-equilibrium,"  while  river  engineers  would  say  that 
^uch  a  channel  is  "in  regime"  (sometimes,  "in  regimen").  A  stable  channel 
adjusts  its  slope  to  provide  the  velocity  and  discharge  required  for  transpor¬ 
tation  of  the  sediment  load  with  a  given  particle-size  distribution  derived 
from  the  drainage  basin.  Ritter  (1986)  pointed  out  that  even  actively  down¬ 
cutting  streams  are  still  in  "quas i -equi 1 ibr ium"  as  defined  by  their  hydraulic 
geometry.  The  discharge  associated  with  the  graded  condition  of  a  stream  is 
often  referred  to  as  the  dominant  discharge,  which  is  higher  than  the  average 
discharge  "because  the  sediment  characteristics  and  grade  of  the  stream  are 
particularly  influenced  by  higher  stages"  (Morris  and  Wiggert  1972). 

1.  Any  change,  whether  natural  or  human  induced,  in  the  hydrologic, 
climatic,  or  geologic  characteristics  of  a  drainage  basin  can  cause  a  channel 
(or  a  reach  of  a  channel)  to  lose  equilibrium.  In  most  alluvial  rivers,  such 


*  See  paragraphs  69-/9  for  a  complete-  discussion  of  Markin's  work. 
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changes  occur  continually;  therefore,  to  transport  sediment  with  maximum  effi¬ 
ciency,  the  channel  is  always  compensating  for  these  changes  by  eroding  its 
bed  (degradation),  depositing  part  of  its  sediment  load  (aggradation),  or 
changing  its  course  to  increase  or  decrease  its  length  (meandering) .  By  means 
of  such  adjustments,  a  stream  tends  toward  a  certain  hydraulic  geometry  and 
slope  that  has  been  described  by  some  researchers  as  "dynamic  equilibrium" 
(Winkley  1981). 

4.  For  an  alluvial  channel  to  be  in  equilibrium,  the  energy  input  must 
equal  the  energy  expended.  Morris  and  Wiggert  (1972)  have  expressed  the  rela¬ 
tion  as  follows: 

H  =  Hb  +  HL  +  H„  +  Hs  +  H.  <!>* 


who  re 

//  -=  total  energy  available  after  precipitation  (=  elevat  ion  above 
mean  sea  level) 

Hh  -  energy  dissipated  as  heat  through  friction  of  tlow 

Ht  =  energy  lost  by  turbulence  (e.g.,  rapids,  stream  expansions, 
contractions,  structures,  etc.) 

Hn  =  energy  extracted  by  artificial  means  (e.g.,  turbines) 

Hs  =  energy  used  to  transport  sediment  and  debris  loads 

He  -  energy  expended  in  erosion  or  rocks,  soils,  etc. 

In  an  ideal  situation  where  a  stream  is  graded  and  flowing  at  dominant  dis¬ 
charge,  only  the  sediment  and  debris  supplied  by  the  basin  and  nothing  more 
will  be  transported;  the  He  term  would  then  become  zero.  However,  such  a 
system  does  not  commonly  exist  in  nature. 

Purpose  and  Approach 

5.  Researchers  do  not  specifically  agree  on  the  mechanics  and  effects 
of  suspended- sediment  transport  as  they  relate  to  the  character  of  an  alluvial 
river.  Since  the  end  of  the  last  century,  many  theories  have  been  formulated. 
The  purpose  of  this  study  was  to  survey  pertinent  regime  equations  to  evaluate 
the  effects  of  changes  in  suspended  -  sediment  loads  as  they  relate  to  the 
regime  of  alluvial  rivers.  This  task  was  approached  by  researching  t he 


For  convenience,  symbols  are  listed  and  defined  in  the  Notation 
(Appendix  A ) . 
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hydraulic  engineering  literature,  wherein  pertinent  investigat 
jfied.  These  studies,  which  dealt  with  the  effects  of  changes 
sediment  loads  on  hydraulic  parameters  and  bed- load  transport 
alluvial  streams,  are  summarized  and  presented  herein. 


ions  were  i dent¬ 
in  suspended- 
capabilities  in 
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PART  II:  RESEARCH  AND  DEVELOPMENT  OF  REGIME  EQUATIONS 


6.  A  number  of  researchers,  beginning  with  Robert  Greig  Kennedy  of 
India  in  1895,  have  attempted  to  quantitatively  describe  regime  channels. 

This  part  of  the  report  presents  a  chronological  overview  of  the  21  selected 
investigations  and  provides  discussions  of  the  significant  equations  that 
generally  relate  changes  in  suspended  -  sediment  loads  to  streamflow  regime. 

The  terms  used  in  equations  have  been  standardized  (see  Appendix  A)  insofar  as 
practical  with  those  used  by  the  American  Society  of  Civil  Engineers  (Vanoni 
19/5).  Any  discrepancies  from  the  original  references  are  explained  in 
t  oo  t  no  t.es  . 


Kennedy  (  1  895 ) 

7.  Kennedy  has  been  credited  as  being  the  first  person  to  develop  equa¬ 
tions  for  canals  in  regime.  The  observations  upon  which  Kennedy  based  his 
results  were  made  on  the  lower  reaches  of  the  Bari  Doab  canal  system.  The 
study  area  included  a  90-mile'"  reach  of  the  main  canal  and  12  distributing 
channels.  The  maximum  discharge  rate  in  the  main  canal  was  1,700  cfs,  whereas 
the  12  distributing  channels  had  capacities  ranging  from  30  to  250  cfs.  The 
channels  of  the  Bari  Doab  canal  system  were  thought  to  have  formed  permanent 
cross  sections  by  silting  and  scouring  during  the  many  years  of  steady  flow. 
These  reaches  had  never  been  dredged;  therefore,  for  a  considerable  time  the 
sediment- transporting  capability  of  each  had  been  sufficient  to  handle  the 
sediment  introduced  to  the  system.  Each  channel  "...  could  have  carried 
neither  more  nor  less  silt  under  the  given  conditions  of  width  of  bed,  depth 
of  channel,  and  velocity  of  flow"  (Kennedy  1895).  This  canal  system,  believed 
to  be  in  equilibrium,  was  found  to  have  nearly  rectangular  cross  sections  with 
vertical  sides  of  fine  sediment  and  a  horizontal  bed  of  coarse  sand.  Kennedy 
(1898)  postulated  there  should  be  a  relation  between  bed  width,  depth,  and 
velocity  ensuring  that  a  channel  would  transport  its  load  without  deposition. 

8.  Kennedy  (1895)  derived  the  mean  velocity  of  each  channel  from  the 
known  "full  supply"  discharge  and  the  measured  cross  sections.  The  "full 
supply"  discharge  was  somewhat  greater  than  the  average  discharge  but  less 


*  A  table  of  factors  for  converting  non-Sl  units  to  SI  (metric)  units  is 
presented  on  page  A. 
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than  the  flood  discharge  and  was  generally  a  constant  ratio  to  the  average  in 
all  the  channels.  The  channels  were  very  dissimilar  with  respect  to  width- 
depth  ratios  (b/d)  ,  which  varied  between  about  4  and  15.  By  determining  the 
line  of  best  fit  through  a  logarithmic  plot  of  depth  of  flow  versus  mean 
velocity,  Kennedy  (1895)  developed  the  following  relation: 

Va  =  cdm  (2) 


where 


!7q  =  mean  velocity 
c  =  0.84 

d  =  average  depth 
m  =  0.64 

He  also  concluded  that  bed  width  did  not  affect  the  above  relation. 

9.  Kennedy  (1895)  stated  that  F0  could  be  called  the  "critical 
velocity"  or  that  velocity  at  which  silting  is  just  prevented  for  the  given 
depth.  Thus,  to  avoid  silting  in  any  canal  system,  V0  must  not  be  constant 
in  all  the  channels,  but  rather  must  increase  with  depth  according  to  Equa¬ 
tion  2.  When  the  soil  properties  permit,  a  velocity  greater  than  W0  (but 
never  a  lesser  velocity)  may  be  used  for  design  purposes.  Kennedy  stated  that 
on  different  canal  systems,  the  values  of  c  and  m  could  vary,  but  probably 
by  only  a  small  amount. 

10.  Kennedy's  (1895)  concepts  of  sediment  transport  probably  resulted 
from  the  early  tendency  in  canal  design  to  make  the  channels  narrow  and  deep. 
He  postulated  that 

sediment  in  a  flowing  canal  is  kept  in  suspension  solely  by  the 
vertical  components  of  the  constant  eddies,  which  can  always  be 
observed  in  any  stream,  boiling  up  gently  to  the  surface.  From 
the  sides  also,  some  such  eddies  may  occur  to  a  much  smaller 
degree,  but  any  such  must  be  for  the  greater  part  horizontal,  and 
of  no  silt- supporting  power.  In  order,  therefore,  to  obtain  an 
expression  for  the  silt -supporting  power  of  the  stream,  it  may 
safely  be  assumed  that  the  quantity  of  silt  supported  will  be  pro¬ 
portional  to  the  width  of  the  bed,  all  other  conditions  remaining 
the  same.  It  must  also  vary  with  the  velocity  of  the  stream  V  , 
say  as  Vn_1  ,  since  the  greater  the  velocity  the  greater  must  be 
the  force  of  the  eddies,  which  become  zero  when  the  velocity  is 
zero.  There  is  a  third  variable,  the  depth,  but  this  could  affect 
neither  the  number  nor  the  force  of  the  eddies.  The  amount  of 
silt  supported  in  a  stream  may,  therefore,  be  expressed  by 
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abV11'1  ,  where  a  is  some  constant  unknown.*  But  whilst  thus 
supported,  the  silt  is  being  moved  forward  at  velocity  V  ,  so 
that  the  amount  of  silt  transported  will  be  equal  to  abVn  .  It 
is  here  presumed  that  all  sediment  is  in  suspension,  but  there  is 
doubtless  a  small  portion  of  the  heavier  silt  simply  rolling  along 
the  bed.  This  amount  would  vary  as  bV  ;  instead  of  as  bVn  ;  so 
that  the  value  of  n  to  include  the  rolling  silt  will  be  somewhat 
less  than  it  would  be  if  the  suspended  silt  alone  were  considered. 

If  it  be  assumed  that  the  amount  of  silt  supported  is  proportional 
to  the  upward  pressure  of  the  deflected  currents  of  water,  which 
varies  simply  as  the  square  of  the  velocity,  or  as  V2  ,  the 
expression  n  -  1  =  2  ,  or  n  =  3,  is  arrived  at.  From  these 
considerations,  therefore,  allowing  for  rolling  silt,  a  value 
something  less  than  3  would  be  derived  for  n  . 

The  derivation,  omitted  here,  can  be  found  in  Kennedy  (1895). 

11.  The  formula  Kennedy  derived,  along  with  Kutter's  flow  formula  that 

he  used,  both  represented  only  2  degrees  of  freedom.  Kennedy  did  not  derive  a 

width  formula,  which  meant  that  either  a  deep  narrow  channel  or  a  wide  shallow 

channel  could  be  designed  to  carry  the  same  discharge.  Later  Kennedy  (1904) 

reassessed  his  research  and  developed  a  "rough-rule"  relation  for  width-depth 

ratio  (Graf  1971;  Lindley  1919),  thus  satisfying  three  degrees  of  freedom. 

Kennedy  further  developed  his  "silt  theory"  in  his  "Hydraulic  Diagrams" 

(Kennedy  1907) . 


Lindlev  (1919) 


12.  In  1895,  Kennedy  had  made  the  first  step  toward  determining  regime 
theory.  In  1917,  F.  W.  Woods  followed  Kennedy's  lead  by  publishing  a  list  of 
channel  dimensions  in  "Normal  Data  of  Design  for  Kennedy  Channels"  (Woods 
1917).  Earlier,  Kennedy  (1907)  had  claimed  that  the  degree  of  "silt  charge" 
(or  sediment  concentration)  is  the  same  for  all  channels  in  any  one  system. 
Lindley  (1919)  pointed  out  that  Kennedy's  statements  contained  a  fallacy  and 
that  neither  Kennedy  nor  Woods  (1917)  presented  any  data  to  support  what  they 
had  published;  therefore,  Lindley  (1919)  saw  a  need  for  retesting  these 
earlier  theories.  He  stated  that  he  found  "very  few  small  channels  of  depths 
as  great  as  Mr.  Woods  specifies,  while  some  large  channels  considerably  exceed 
his  limit.  Kennedy’s  rules  were  transgressed  in  design  and  in  nature  because, 
according  to  his  Hydraulic  Diagrams  (Kennedy  1907),  actual  large  channels 
should  silt,  and  most  small  channels  scour"  (Lindley  1919).  Lindley  (1919) 

*  Kennedy  (1895)  used  A  instead  of  a  . 
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then  suggested  that 

channels  carrying  an  average  silt  charge  .  .  .  had  not  attained 

regime,  so  that  (Kennedy's)  .  .  .  selection  of  stable  channels 

showed  unconscious  bias  in  favor  of  regimes  of  low  silt  charges. 

The  writer  has  found  great  local  and  general  variation  of  dimen¬ 
sions  even  on  the  same  channel,  and  it  is  not  certain  to  what 
extent  Kennedy  succeeded  in  choosing  sites  that  represented  fair 
averages  on  the  channels  he  selected.  It,  therefore,  appears 
necessary  to  retest  Kennedy's  theory  for  depths,  and  to  fix  the 
law  of  variation  of  widths  .  .  .  (using  existing  data  and)  .  .  . 

entirely  empiric  methods,  so  as  to  avoid  bias  in  the  choice  or 
interpretation  of  data. 

13.  By  examining  data  from  so-called  regime  channels,  Lindley  (1919) 
found  large  variations  to  which  he  assigned  the  following  probable  causes: 

a.  "Silt-clearance ,  berm  cutting  (as  distinguished 
from  mere  trimming  of  irregularities),  the  use 
of  bushing,  and  the  working  of  silting  tanks.* 

b.  "Use  of  regulators  in  heading  up  more  or  less  at 
different  times. 

c.  "Running  of  feeder  channels  with  varying 
discharges . 

d.  "Breaches  which  scour  out  the  bed  locally. 

e.  "Differences  between  the  silt-drawing  power  of 
offtakes,  not  at  present  measurable,  but  causing 
an  uncertain  modification  of  regime  in  the 
parent  channel  and  setting  an  uncertain  regime 
in  the  distributary. 

f.  "The  apparent  tendency  of  silt  to  move  in  waves, 
so  that  a  channel  contains  a  succession  of 
scouring  and  settling  reaches. 

g.  "Variation  of  silt-charge  at  different  seasons 
in  main  canals,  diminishing  towards  the  lower 
ends  of  the  system  of  channels. 

h.  "Effect  of  different  rugosities  (or  hydraulic 
roughness  of  the  streambed)  in  allowing  regime 
velocities  to  be  attained  with  lesser  or  greater 
gradients . " 

19.  Lindley  (1919)  analyzed  data  from  hydraulic  surveys  made  between 
1915  and  1917  of  the  Lower  Chenab  Canal,  a  system  of  some  2,700  miles  of 
channel  having  a  combined  discharge  of  11,000  cfs.  He  analyzed  widths  between 
point  bar  deposits,  depths,  and  gradients,  averaged  over  1-mile  lengths,  using 


*  "Berm,"  as  used  by  Lindley  (1919),  was  a  point  bar  deposit.  "Bushing" 

probably  refers  to  the  use  of  planted  vegetation  to  control  erosion.  "Silt¬ 
ing  tanks"  would  be  synonymous  with  stilling  basins. 
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1  mile  of  every  3-  to  5-mile-long  reach,  with  the  exception  of  "irregular 
reaches"  (Lindley  1919).  Lindley  then  plotted  average  width  against  average 
depth  on  a  logarithmic  scale  and  noted  the  gradients  of  the  channel  at  each 
point.  Next,  the  widths,  depths,  and  gradients  of  each  branch  system  were 
grouped  and  averaged.  The  average  width-depth  relation  could  be  expressed  as 

=  3.8  d1-61  O)* 


where 

bb  =  bed  width  developed  from  water  surface  width  under  the  assumption 
that  side  slopes  are  IV: 2H 

This  relation  was  valid  for  up  to  9  ft  of  depth  (Lindley  1919)  ,  with  no  indi¬ 
cation  of  a  change  beyond  that  depth.  Although  the  plotted  gradients  did  not 
follow  an  exponential  law,  velocities,  calculated  using  Kutter's  relation  with 
N  =  0.225  ,  were  found  to  follow  the  exponential  law: 

Va  =  0.95  d0-57  (4) 


rather  than  the  relation  developed  by  Kennedy  (1895)  (Equation  2). 

15.  Lindley  (1919)  pointed  out  that  Kennedy's  velocities  are  "  'derived 
from  the  known  full  supply  discharge  and  the  measured  cross-section  area  of 
the  channel,'  while  the  velocities  of  the  present  investigation  are  those 
assumed  in  designing  (i.e.,  design  procedures);  and  as  long  as  Kutter’s 
formula,  with  N  =  0.0225  ,  is  used,  the  new  relation  should  be  followed." 
Although  he  did  not  suggest  that  Kennedy's  relation  was  wrong,  Lindley  (1919) 
did,  however,  point  out  that  all  existing  formulas  "assume  a  uniform  rugosity 
in  all  channels  covered  by  one  value  of  coefficient."  He  contended  that 
silt- formed  channels  probably  demanded  a  regular  variation  of  coefficient  with 
size  and  with  changing  stages. 

16.  By  combining  the  two  relations  between  width  and  mean  velocity  with 
respect  to  depth,  Lindley  (1919)  obtained  the  following  relations  between 
width  and  mean  velocity: 

for  Kutter  velocities  V0  =  0.59  bb  35  (5) 


*  Lindley  used  b  instead  of  bb  . 
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for  Kennedy  velocities 


(6) 


VQ  =  0. 50  bb0  i,° 


Theoretically,  a  certain  velocity  along  the  sides  of  a  channel  is  required  to 
maintain  a  balance  between  erosion  and  deposition.  According  to  Lindley's 
formula,  the  wider  the  channel,  the  greater  the  mean  velocity  permitted.  The 
depth  will  adjust  itself  to  the  velocity.  Lindley  (1919)  defended  this  theory 
by  saying: 

While  it  is  known  that  greater  or  lesser  silt- 
charges  demand  greater  or  lesser  velocities  for  equal 
depth,  the  case  of  widths  is  not  necessarily  similar: 
berms  forming  the  sides  of  channels  are  composed  of  such 
fine  muddy  silt,  so  easily  suspended  in  flowing  water, 
that  it  (the  silt)  flows  equally  all  through  a  canal 
system,  and  berms  formed  of  the  same  silt,  at  the  same 
velocity,  should  be  similar  in  character. 

He,  therefore,  concluded  that  of  the  two  channels  carrying  the  same  discharge, 
the  one  with  the  greater  "silt  charge"  must  be  wider  and  shallower  and  must 
have  a  steeper  gradient.  Lindley  (1919)  had  advanced  the  quantitative  analy¬ 
sis  of  regime  channels  by  defining  a  relation  between  velocity  and  width  that 
added  to  the  already  established  relations  of  velocity  and  depth. 

Lacev  (1930) 

17.  Prior  to  the  time  of  Gerald  Lacey,  several  researchers,  including 
Lindley,  had  continued  the  work  of  Kennedy  (paragraph  8)  by  correlating  mean 
velocity  and  depth  of  flow;  however,  the  coefficient  c  varied  between  0.67 
and  0.95  and  the  exponent  m  between  0.52  and  0.73,  depending  on  the  channel 
studied  (Henderson  1966,  Lacey  1930).  Lacey  (1930)  stated  that  V0  was  a 
function  of  both  depth  and  a  "silt  factor"  and  that  width  was  as  important  as 
depth  in  determining  this  velocity.  He  defined  the  silt  factor  as 


^o/^ok  =  ratio  of  the  regime  velocity  actually  obtained  with  a  given 
depth  to  the  velocity  for  the  same  depth  that  Kennedy's 
formula  (Equation  2)  would  predict  (Lacey  1930) 

Interpreting  d  in  terms  of  R  ,  the  hydraulic  radius,  he  developed  his 

first  general  formula: 

V0  =  1.17  JfR 


This  silt  factor  had  a  very  simple  scalar  significance.  The  hydraulic  radius 
of  the  channel  would  be  1/f  times  the  hydraulic  radius  of  a  standard  silt 
channel  (Kennedy  channel)  flowing  with  the  same  mean  velocity  (Lacey  1930) . 

18.  Lacey  (1930)  stated  that  if  Lindley  had  dealt  with  mean  bed  width 
instead  of  bed  widths,  the  product  of  width  and  depth  (of  flow)  would  have 
given  a  useful  expression  for  the  cross-sectional  area  of  the  channel.  Fur¬ 
thermore,  Lindley  dealt  with  depth  instead  of  the  hydraulic  radius,  and 
because  d  steadily  approaches  R  as  the  width  of  the  stream  increases,  the 
structure  of  Lindley' s  formula  was  useful  for  designed  channels.  However,  the 
formula  did  not  permit  extended  application  by  extrapolation.  Lindley's 
formulas  for  high  discharges  yielded  bed  widths  that  were  admittedly  too  low, 
and,  for  this  reason,  his  suggestion  that  a  channel  of  given  discharge  and 
"silt  charge"  was  uniquely  determined  remained  to  be  proved  (Lacey  1930) . 

19.  Lacey  successfully  developed  a  formula  that  confirmed  Lindley's 
theorem.  The  velocity-discharge  relation  presents  itself  in  the  simplest  form 
possible  if  the  cross-sectional  area  and  the  mean  velocity  are  first  con¬ 
sidered.  As  shown  in  Figure  1,  if  two  different  channels,  abc  and  a'b'c'  , 
have  the  same  mean  velocity  but  different  "silt  factors,"  then 


and 


fR  =  f'R' 

fP  =  f'P' 


where  P  and  P'  are  the  wetted  perimeters  of  channels  abc  and  a'b’c’  , 
respectively.  If  the  above  equalities  are  multiplied  together,  they  can  be 
reduced  to 


Af2  =  A'f'2 


U 


Figure  1.  Channels  abc  and  a'b'c'  have  the  same  mean 
velocity  but  different  silt  factors,  f  and  f'  , 
respectively.  Therefore,  fR  =  f'R'  and  f?  =  f'P' 

(adapted  from  Lacey  1930) 

where  A  and  A'  are  cross-sectional  areas  for  channels  abc  and  a'b'c'  , 
respectively . 

20.  This  relation  assumes  complete  geometric  similarity  in  the  sec¬ 
tions.  Lacey  (1930)  stated  that  the  sediment  distributed  around  the  wetted 
perimeter  does  not  necessarily  have  to  be  of  the  same  gradation  throughout  any 
one  section.  All  that  is  necessary  is  that  the  grading  of  the  sediment  from 
a  to  b  in  the  first  section  is  identically  similar  in  all  respects  to  the 
grading  from  a'  to  b'  in  the  second  section.  To  test  the  relation  of 
area,  silt  factor  and  velocity,  Lacey  (1930)  plotted  Af2  against  VQ  and 
determined  that 

Af 2  =  3 . 8  F05 


This  could  also  be  written, 

Qf 2  =  3 . 8  V* 

where  Q  =  discharge.  Finally  Lacey  (1930)  developed  an  equation  that  related 
wetted  perimeter  to  discharge: 

P  =  2.668  Q°  5 
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21.  Equation  9  shows  that  for  a  given  discharge,  the  wetted  perimeter 
of  a  stable  channel  is  constant  and  independent  of  the  size  of  the  material 
transported,  with  the  silt  size  merely  controlling  the  "shape."  If  channels 
with  horizontal  beds  are  considered,  then  "shape"  is  the  ratio  of  bed  width  to 
depth.  For  channels  having  semielliptical  cross  sections,  shape  is  the  ratio 
of  the  major  axis  to  the  semiminor  axis.  The  wetted  perimeter  is  a  constant 
for  any  given  discharge,  but  the  shape  of  the  cross  section  depends  on  the 
size  of  the  bed  material  (Figure  2).  The  coarser  the  bed  material,  the  flat¬ 
ter  the  semiellipse  and  the  greater  the  width  of  water  surface.  The  finer  the 
bed  material,  the  more  nearly  the  section  approximates  a  semicircle.  Thus, 
Lacey  (1930)  concluded  that  two  elementary  principles  that  govern  the  dimen¬ 
sions  of  stable  channels  in  alluvium  were: 

All  stable  s ilt- transporting  channels  flowing  with  the 
same  mean  velocity  tend  to  assume  the  same  "shape,"  and 
ail  stable  silt- transporting  channels  of  the  same  dis¬ 
charge  have  the  same  length  of  wetted  perimeter  irre¬ 
spective  of  the  particular  grade  of  silt  carried. 

Lacey  (1930)  used  the  term  "silt"  very  loosely,  and  in  the  above  statement, 

"silt"  could  have  also  included  other  particle  sizes. 


Figure  2.  According  to  Lacey  (1930),  cross-sectional  shape 
is  a  function  of  the  size  of  the  bed  material.  Coarser  bed 
material  yields  a  flatter  semiellipse,  whereas  finer  bed 
material  yields  an  approximate  semicircle  (adapted 

from  Lacey  1930) 

22.  The  new  formulas  developed  by  Lacey  covered  a  very  large  range  of 
stable  artificial  channels.  If  truly  general,  these  formulas  should,  when 
extrapolated,  fit  the  maximum  discharges  of  stable  rivers.  When  the  dis¬ 
charges  of  large  alluvial  rivers  are  considered,  the  wetted  perimeter  closely 
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approximates  the  width  of  the  water  surface;  for  the  purpose  of  such  an 
extrapolation,  the  difference  in  wetted  perimeter  and  width  is  negligible, 
regardless  of  the  precise  shape  that  the  channel  assumes.  From  Equation  9, 
Lacey  (1930)  developed  the  following  formula  to  determine  the  approximate 
width  of  a  stable  alluvial  channel: 


bt  =  2.67  Q0-5 


(10)* 


where 

bt  =  top  width,  ft 

23.  Lacey  (1930)  then  attempted  to  relate  the  silt  factor  to  Rutter' s 
rugosity  (roughness)  coefficient  N  ,  developing  the  relation: 

N  =  0.022  f0  2 


He  stated  that  Equation  11  could  apply  equally  to  Manning's  formula,  with  Man¬ 
ning's  n  essentially  being  interchangeable  with  N  .  By  taking  Manning's 
formula  for  velocity  and  substituting  the  appropriate  formulas  for  the  slope 
of  the  energy  grade  line,  he  obtained  the  following  equation: 


f1.51 

2587  Q1/9 


(12)* ** 


where  Se  =  slope  of  the  energy  grade  line. 

24.  Lacey  (1930)  also  researched  the  relation  between  the  sediment  size 
and  the  silt  factor.  He  found  that  the  average  diameter  of  a  sand  grain  in  a 
sample  of  Upper  Bari  Doab  Silt  (standard  f  =  1.0)  from  the  head  reach  of  a 
distributary  was  0.4  mm  (probably  slightly  larger  than  the  size  of  sand  grains 
in  Kennedy's  channels).  The  diameter  of  the  average  large  boulder  transported 
by  heavy  floods  on  the  Song  River  (f  =  40)  was  approximately  25  in.  He  devel¬ 
oped  the  following  formula  to  provide  a  very  rough  estimate  of  the  diameter  in 
inches  of  the  predominant  size  of  sediment  transported: 


*  Lacey  (1930)  used  W  instead  of  bt  for  top  width. 

**  Lacey  (1930)  used  S  instead  of  Se 
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(13)* 


D  = 


f2 

TZ 


Lacey  (1930)  stated  that  if  he  had  been  able  to  collect  data  on  the  size  of  a 
silt  particle  corresponding  to  a  silt  factor  of  about  0.4  (instead  of  40),  he 
could  have  deri-'ed  a  more  reliable  rule.  With  low  silt  factors,  the  silt  has 
a  remarkably  uniform  character,  and  a  scalar  relationship  is,  therefore,  more 
easily  obtained. 

25.  The  general  formulas  developed  by  Lacey  have  a  drawback  that  is 
independent  of  the  gradation  of  sediment.  Decreasing  the  discharge  with  a 
constant  silt  factor  causes  the  cross  section  to  approach  a  semicircle. 

Hence,  at  the  limit 

A  =  6.2832  R2  =  3.8  V$f2 


V0  =  0.882  fps  ( for  any  type  of  sediment ) 


26.  For  velocities  less  than  this  limiting  value,  the  resultant  wetted 
perimeter  cannot  contain  the  cross-sectional  area.  This  fact  can  be  consid¬ 
ered  either  as  a  limitation  of  the  formulas  or,  as  Lacey  (1930)  concluded,  as 
an  argument  that  this  velocity  must  be  one  of  the  physical  constants  of  water 
because  it  is  independent  of  the  sediment  size.  Lacey  (1930)  summarized  his 
conclusions  as  follows: 

a.  "For  a  given  discharge  and  given  silt  factor,  the  Lindley 
theorem  is  true;  that  is,  the  cross-sectional  area,  wetted 
perimeter,  and  slope  of  a  stable  channel  flowing  in  and 
transporting  its  own  silt  are  uniquely  determined. 

b.  "The  wetted  perimeter  of  a  stable  silt- transporting  channel 
varies  as  the  square  root  of  the  discharge  and  is  indepen¬ 
dent  of  the  type  of  silt  transported,  all  stable  channels 
of  the  same  discharge  having  the  same  wetted  perimeter,  and 
the  silt  factor  determining  the  shape. 


*  Lacey  (1930)  used  d  instead  of  D 
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c.  "All  stable  s i It- transpcr t Ing  channels  of  the  same  stable 
mean  velocity  have  geometrically  similar  shapes,  the  silt 
factor  determining  the  only  difference,  which  is  one  of 
scale  .  " 


Lane  (1935; 

27.  In  1935,  E.  W.  Lane  did  extensive  research  in  regime  theory  to  aid 
in  the  design  of  the  All-American  Canal  in  California.  He  classified  the  two 
types  of  formulas  previously  developed  for  stable  channels  in  India  and  in 
Egypt  as: 

a.  Those  giving  an  expression  of  velocity. 

b.  Those  giving  stable  channel  shapes. 

28.  The  first  type  of  formula,  which  is  similar  to  that  developed  by 
Kennedy  (1895)  (Equation  2),  indicates  that  critical  velocity  increases  with 
depth.  However,  later  experience  showed  that,  as  depth  increases,  the 
velocity  at  which  bank  erosion  occurred  is  also  reached.  Kennedy  believed 
that  determining  this  limiting  velocity  was  a  matter  of  experience. 

29.  Several  researchers  contributed  formulas  of  the  second  type.  Woods 
(1917)  developed  relations  of  mean  depth,  velocity,  and  slope  (see  para¬ 
graph  12),  whereas  Lindley  (1919)  provided  a  relation  of  critical  velocity  to 
depth  and  bed  width  (Equations  5  and  6).  Neither  Woods  nor  Lindley  suggested 
any  modifications  for  the  quantity  or  gradation  of  sediment.  Lacey  (1930) 
used  channel  shapes  and  velocities  and  introduced  the  effect  of  particle  size 
but  did  not  consider  the  quantity  of  the  material  transported.  Lane  (1935) 
analyzed  the  regime  equations  that  had  already  been  developed  and  concluded 
that  none  of  the  existing  formulas  would  be  safe  to  adopt  for  the  design  of 
sections  of  the  All-American  Canal.  He  stated  that 

although  these  formulas  no  doubt  provide  workable  relations  for 
the  conditions  for  which  they  were  developed,  these  conditions 
have  not  been  delineated  sufficiently  to  enable  them  to  be  applied 
elsewhere.  In  general,  also,  they  were  developed  empirically  from 
a  very  limited  range  of  conditions,  and  in  most  cases  they  omit 
important  factors  from  consideration  (Lane  1935) . 

30.  Lane  suggested  analysis  of  the  fundamental  interrelations  of  a 
number  of  factors  that  he  attributed  to  controlling  the  shape  of  streams  in 
alluvial  material.  These  factors  included: 

a.  Hydraulic  parameters  (energy  slope,  roughness,  hydraulic 
radius,  depth,  mean  velocity,  velocity  distribution,  and 
temperature) . 
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b.  Shape  of  the  existing  channel  (width,  depth,  and  steepness 
of  side  slopes) . 

c.  Nature  of  material  transported  (size,  shape,  specific  grav¬ 
ity,  dispersion,  quantity,  as  well  as  bank  and  subgrade 
material) . 

d.  Other  factors  (alignment,  uniformity  of  flow,  and 
"aging"*) . 

31.  Of  the  hydraulic  factors,  energy  slope,  roughness,  hydraulic 
radius,  and  mean  velocity  are  interdependent,  and  their  relation  can  be  deter¬ 
mined  quantitatively  through  the  velocity  formulas  available  at  that  time. 

The  effect  of  the  movement  of  material  in  suspension  and  by  traction  upon  the 
roughness  was  not  known;  however,  when  compared  with  the  uncertainty  in  other 
phases  of  the  problem,  the  relations  among  these  four  elements  are  certain 
enough  that  further  research  into  them  would  have  probably  been  unjustified. 
Lane  (1935)  regarded  the  steepness  of  side  slopes  as  relatively  unimportant 
except  that  the  steepest  banks  were  subject  to  sloughing.  He  believed  that 
both  the  velocity  distribution  and  the  mean  velocity  were  of  primary  impor¬ 
tance.  Although  temperature  had  been  suggested  as  having  an  effect  because  it 
influences  the  viscosity  of  the  water,  its  impact  on  cross-sectional  shape  was 
thought  to  be  minimal  in  relation  to  other  factors  influencing  stable  chan¬ 
nels.  Lane  and  other  researchers  prior  to  him  considered  the  size  of  material 
transported  to  be  of  major  importance  in  controlling  cross-sectional  shape 
(with  particle  shape  and  specific  gravity  being  secondary).  In  addition,  he 
postulated  that  the  quantity  of  material  in  motion,  the  scour  resistance  of 
the  material  comprising  the  banks  and  subgrade,  the  channel  alignment,  vari¬ 
able  discharges,  and  "aging"  were  among  other  factors  affecting  cross- 
sectional  shape  of  a  stable  channel  (Lane  1935) . 

32  Lane  (1935)  stated  that 

the  velocity  along  the  bottom  of  a  stable  channel  must  be  suffi¬ 
cient  to  move  the  quantity  of  material  supplied  to  it  but  not  so 
great  as  to  scour  the  subgrade.  .  .  .  The  science  of  hydrodynam¬ 
ics  has  not  yet  progressed  to  the  point  where  the  relation  between 
the  velocity  distribution  adjacent  to  a  surface  can  be  related 
quantitatively  to  the  drag  of  the  water  along  the  surface  or  to 
the  velocity  gradient  adjacent  to  the  surface. 


*  By  "aging,"  Lane  (1935)  could  have  been  referring  to  the  combined  effects 
of  such  processes  as  armoring,  cementation,  and  the  interaction  of  organisms 
with  the  substrate  (i.e.,  the  bottom  sediment  of  the  streambed) .  "Aging" 
could  also  have  referred  to  cross-sectional  shape  and  character  and  develop¬ 
ment  of  the  channel  and  floodplain. 
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He  concluded  that  when  high  velocities  extend  close  to  a  water  surface,  the 
pushing  or  dragging  force  of  water  on  the  bed  was  greater  than  if  the  veloci¬ 
ties  close  to  the  surface  were  low  (Lane  1935)  . 

33.  Researchers  generally  agreed  that  vertical  currents  stronger  than 
the  force  of  gravity  kept  material  in  suspension.  This  led  Lane  (1935)  to 
state: 

A  promising  hypothesis  for  the  capacity  of  a  stream  to  transport 
material  in  suspension,  therefore,  should  be  that  the  capacity  is 
proport’onal  to  its  turbulence,  which,  in  turn,  is  probably 
proportional  to  the  energy  expended.  The  concentration  of  a  given 
quality  [sic]*  of  solids  which  a  stream  could  support,  therefore, 
would  be  proportional  to  the  energy  expended  per  unit  of  volume  of 
the  water.  This  energy  is  proportional  to  the  rate  of  fall  of  the 
water,  which  is  equal  to  the  product  of  the  velocity  and  the 
slope.  This  is  a  kind  of  over-all  relation,  however,  and  silting 
may  occur  in  one  part  of  a  ditch  cross-section  while  other  parts 
may  be  scouring. 

34.  Channels  frequently  carry  a  considerable  amount  of  solid  matter  by 
dragging  it  along  the  bed,  with  either  clear  or  sediment- laden  water  flowing 
above.  The  quantity  of  material  moved  depends  on  the  velocity  near  the  bottom 
of  the  channel.  Therefore,  a  channel  supplied  with  a  heavy  bed  load  requires 
high  velocities  to  move  the  bed  load  through  the  reach  and  to  avoid  sediment 
deposition.  Lane  (1935)  concluded  chat  channels  carrying  bed  load  should  have 
high  ratios  of  bed  width  to  depth  to  prevent  the  sides  of  the  channel  from 
being  eroded  by  the  high  velocity  required  to  carry  a  heavy  bed  load  in  a  nar¬ 
row  channel.  He  stated  that  if  the  banks  are  composed  of  a  scour- resistant 
material,  the  ratio  of  bed  width  to  depth  can  be  reduced. 

35.  Lacey  (1930)  had  concluded  that  stable  channels  would  be  semi¬ 
elliptical,  with  the  major  axis  horizontal;  the  ratio  of  the  major  axis  to  the 
semiminor  axis  would  depend  on  the  nature  of  the  silt  carried,  and  this  ratio 
would  be  greater  for  coarser  material.  Lane's  (1935)  research  did  not  support 
this  conclusion,  however.  He  used  the  "form  factor,"  the  ratio  between  the 
area  of  the  channel  section  up  to  the  water  surface  and  the  area  of  the 
enclosing  rectangle,  for  his  analysis.  For  ellipse,  the  form  factor  would 
be  0./9;  for  a  parabola,  0.67;  for  a  triangle,  0.50;  and  for  a  rectangle, 

1.00.  Those  canals  studied  by  Kennedy  carried  a  heavy  load  of  graded  sediment 
and  were  reported  as  having  nearly  vertical  sides;  therefore,  according  to 
Lane  (1935),  they  would  have  a  form  factor  of  approximately  1.00: 


*  Should  be  "quantity. 
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Just  what  shapes  for  stable  channels  are  produced  by  all  varieties 
of  conditions  has  not  yet  been  determined,  but  the  writer  has 
observed  that  for  channels  carrying  a  heavy  load  of  graded  silt, 
ranging  from  colloids  to  fine  sand,  the  sections  have  nearly 
horizontal  beds  composed  of  the  fine  sand  and  nearly  vertical 
sides  of  silts  and  clays.  Such  channels  have  form  factors  of 
about  0.90.  This  is  the  condition  on  the  canals  of  the  Imperial 
Valley.  This  difference  between  the  composition  of  the  bed  and 
bank  material  has  also  been  observed  in  India.  Similar  conditions 
result  in  channels  carrying  considerable  bed  load  and  a  moderate 
quantity  of  fine  silt.  For  a  channel  carrying  water  containing  a 
small  quantity  of  silt  at  high  velocity  in  a  material  containing  a 
considerable  number  of  cobbles,  the  cross-section  is  distinctly 
saucer-shaped,  most  of  the  section  being  covered  with  cobbles,  but 
with  a  small  silt  berm  at  each  edge.  This  condition  was  observed 
on  some  of  the  ditches  on  the  Uncompahgre  Project,  in  Colorado, 
and  in  the  San  Luis  Valley,  in  Colorado.  One  ditch  that  was  care¬ 
fully  measured,  had  a  form  factor  of  0.85. 

Lane  (1935)  believed  that  further  study  along  these  lines  would  disclose  typi 

cal  shapes  and  the  reasons  for  a  number  of  common  conditions. 


Lacey  (19361 


36.  In  discussing  Lane's  (1935)  article,  Gerald  Lacey  emphasized  that 
channel  conditions  and  the  problems  associated  with  studying  regime  flow  were 
much  simpler  in  India  than  on  the  Colorado  River.  Lacey  contended  that  cer¬ 
tain  elements  of  "flume  traction"*  could  not  be  entirely  dissociated  from 
channel  behavior  on  some  of  the  Colorado  canals.  In  channels  with  "flume 
traction,"  there  was  no  limit  to  the  surcharge  that  could  be  forced  down¬ 
stream,  provided  that  slope  was  available  and  that  the  banks  possessed  some 
degree  of  erosion  resistance.  "In  rivers  flowing  in  deep  alluvium,  the  sur¬ 
charge  is  thrown  down,  and  the  river  moves  to  a  flank  and  picks  up  a  fresh 
charge,  which  it  sweeps  forward"  (Lacey  1936). 

37.  Lane  (1935)  had  stated  that  the  capacity  of  a  stream  to  transport 
suspended  sediment  was  probably  proportional  to  its  turbulence  and  to  the 
energy  expended.  Lacey  (1936),  on  the  other  hand,  pointed  out  that,  accordin 
to  Gilbert  (1914),  the  prime  mode  of  sediment  transport  was  saltation  rather 
than  suspension.  With  "flume  traction"  and  a  heavy  surcharge,  the  turbulence 
necessary  to  drive  the  boundary  layer  forward  would  also  occasionally  result 
in  a  high  suspended  charge.  Lacey  (1936)  believed  that  the  ratio  V2/R  for 


*  For  a  discussion  of  a  laboratory  investigation  of  "flume  traction  and 
transportation,"  see  Chang  (1939). 


any  grade  of  silt,  irrespective  of  the  concentration,  was  proportional  to  the 
turbulence.  He  found  that  fine  silt  with  a  heavy  surcharge  would  probably 
demand  the  same  value  of  V2/R  as  would  a  coarse  silt  with  a  smaller  charge 
in  another  locality. 

38.  Lacey  (1936),  building  on  Lane's  (1935)  research,  correlated  Lane's 
factors  as  dimensioned  variables  using  the  Buckingham  theory  (Buckingham 
1914).  This  method  consisted  of  determining  dimensionless  arguments  from  all 
the  known  independent  variables  and  then  correlating  these  arguments.  A  poor 
degree  of  correlation  could  indicate  that  important  variables  had  been  over¬ 
looked  (Lacey  1936).  In  1"69,  Antoine  Chezy  had  selected  four  independent 
variables  that  resulted  in  the  following  equation  (Clemens  1897): 


the  Chezy  number 


(19)* 


where 

p  =  density  of  water 
R0  =  mean  velocity 

39.  Osborne  Reynolds  added  a  fifth  variable,  the  kinematic  viscosity  of 
water.  The  result  was  the  classic  Reynolds  number  (Reynolds  1883);  thus 


■s.fi 

M  '  0  > 


where 

K,  m  =  constants 

RV0/v  =  Reynolds  number  (dimensionless,  in  terms  of  hydraulic  radius, 
with  v  being  the  kinematic  viscosity  of  water) 

The  acceleration  due  to  gravity  g  does  not  enter  this  relation  as  an  inde¬ 
pendent  variable  because  Reynolds  carried  out  his  experiments  with  pipes. 

40.  In  his  earlier  study,  Lacey  (1930)  had  suggested  that  the  rugosity 
was  implicit  in  the  slope  and  the  hydraulic  mean  depth  for  open  channel  and 


*  Lacey  (1936)  used  i  instead  of  SQ  and  V  instead  VQ 
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that  g  had  to  be  added  as  the  sixth  variable.  From  this  addition,  a  third 
dimensionless  variable,  the  Froude  number,  evolved: 


'  ^ 

m 

2 

RVn 

V 

0 

o 

V 

j 

(16)* 


where 


m,  m'  =  constants 

Lacey  (1936),  addressing  the  nature  in  which  the  ratio  V*/R  ,  enters  the 
Froude  number,  concluded: 

The  silt  factor  is  merely  a  simple  proportion  to  the 
ratio,  V*/R  ,  and  therefore,  as  long  as  dimensional 
analysis  is  used  as  an  instrument  in  hydraulic  analysis 
the  ratio,  V*/R  ,  must  persist  as  a  fundamental  ele¬ 
ment,  and  with  it  the  silt  factor,  which  more  appropri¬ 
ately  could  have  been  termed  "turbulence." 

41.  Substituting  for  the  energy  gradient  in  terms  of  the  geometric 

slope,  density,  and  acceleration  due  to  gravity  yields,  Lacey  (1936)  obtained 

the  following  equation: 


SbS* 


=  K 


RV„ 


-1/3 


gR 


2/3 


where 

Sb  =  bed  slope.** 

If  the  silt  factor  is  treated  as  being  proportional  to  l l\/R  ,  this  equation 
can  be  solved  and  will  prove  the  basis  of  Lacey's  formulas  with  the  exception 
of  those  involving  discharge,  thus  leaving  the  silt  factor  in  the  ratio  form 

as 


Voc  g5/6u-1/6f-1'',i?3/<1Sb1/z 


*  Lacey  (1936)  used  p  instead  of  m'  . 

**  Lacey  (  19  36)  used  S  instead  of  Sb  . 
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This  relation  is  true  for  any  regime  channel  in  unconsolidated  alluvium, 
whether  fine  sand,  coarse  sand,  or  boulders.  The  very  low  power  of  the 
kinematic  viscosity  should  be  noted,  as  should  the  complete  disappearance  of 
the  "silt  factor,"  which  is  inherent  in  the  dimensions  assumed  by  the  channel. 
Rewriting  the  formula,  then  substituting  f1/4  for  n  and  introducing  the 
wetted  perimeter  as  a  variable  gives  (Lacey  1936) : 


P 

P 


(17) 


from  which  P  «  Ql/2  and  P/R  <*  V  where  V  =  velocity. 

42.  Earlier  Lacey  (1930)  had  proposed  the  formula: 

P  =  2.668  Q1/z 


from  which 


P2  =  7.12  PRV0 


Then  dividing  through  by  PR  yields 


£-7-12Vo 


(18) 


This  is  Lacey's  fundamental  shape  formula.  Lacey  (1936)  believed  that  Lane's 
(1935)  work  was  unclear  as  to  whether  the  Lacey  shape  formula  had  been 
confirmed  because  Lane  had  presented  only  shape  versus  particle  diameter. 
Lacey  (1936)  believed  that  a  final  solution  to  the  problem  of  silt  transport 
would  be  found  by  investigation  along  the  lines  of  dimensional  analysis  as  he 
had  outlined. 
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Shulits  (1936) 


43.  Fluvial  morphology  deals  with  the  forms  of  riverbeds  and  their 
shapes  in  plan  and  in  profile.  Knowledge  of  fluvial  morphology  is  generally 
limited  to  qualitative  reasoning  and  some  site-specific  empirical  equations; 
however,  a  few  quantitative  principles  have  been  developed  theoretically.  The 
major  purpose  of  Samuel  Shulits'  (1936)  paper  was  to  examine  these  quantita¬ 
tive  laws;  however,  he  first  made  observations  of  the  qualitative  laws  and 
their  use. 

44.  Morphologic  activity  is  most  reflected  by  streambed  slope,  which 
can  be  either  lengthened  or  shortened  by  the  meandering  processes  of  an  allu¬ 
vial  river.  Meandering  is  usually  more  prominent  in  the  lower  reaches  of  a 
river  where  the  bed-material  gradation  shifts  toward  the  fine  fraction  as  a 
result  of  its  travel  downstream  and  where  bed  slope  is  small,  thus  suggesting 
a  reciprocal  relationship  between  size  of  material  and  slope.  Therefore, 
slope,  abrasion,  and  bed  load  are  the  primary  principles  of  fluvial  morphology 
(Shulits  1936)  . 

45.  Shulits  (1936)  stated  that  a  natural  river  can  decrease  its  slope 
by  meandering  or  lengthening  the  present  meanders.  This  decrease  in  slope  is 
often  caused  by  a  decrease  in  bed  load.  On  streams  in  which  a  dam  has  been 
built,  degradation  will  occur  below  the  dam  in  order  to  decrease  the  slope  to 
account  for  the  loss  of  the  sediment  trapped  in  the  reservoir. 

46.  As  particles  move  downstream,  they  are  reduced  in  size  by  abrasion. 
Collision  also  reduces  particle  size,  but  to  a  much  lesser  extent.  Sternberg 
(1875)  developed  a  theoretical  law  for  abrasion,  which  was  later  verified 
experimentally  by  Schoklitsch  (1926,  1933,  1935).  Sternberg  (1875)  had 
assumed  that  the  reduction  in  particle  weight  was  directly  proportional  to  the 
work  done  in  overcoming  friction  along  the  distance  traveled.  The  Sternberg 
abrasion  law  is  in  the  form: 


-c'l t>  X 


(19)* 


*  Shulits  (1936)  used  P  ,  PQ  ,  and  c  instead  of  I7b  ,  Wbo  ,  and  c' 
respect ive ly . 
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where 


Wh  =  weight  of  bed  sediment 
Wbo  =  mean  weight  of  bed  sediment 
e  =  base  of  natural  logarithm 
c'  =  constant  of  proportionality 
4>  =  coefficient  of  friction 
x  =  distance  traveled  downstream 
Setting  a  =  c4>  yields* 


=  W , 


be 


where  a  =  coefficient  of  abrasion.  The  a  ,  which  is  not  a  constant  for  a 
given  rock  type,  depends  on  the  shape,  tends  to  increase  with  the  one -fourth 
power  of  the  velocity  of  the  particle,  and  is  proportional  to  the  diameter  of 
the  bed  material  (Shulits  1936). 

47.  Because  the  slope  of  a  river  varies  directly  with  the  decrease  in 
bed-material  size,  Shulits  (1936)  stated  that  the  next  logical  step  would  be 
to  search  for  a  relation  between  slope  and  bed-material  size.  Expressed 
mathematically 


S 


a  W, 


bo 


where  S  is  the  slope  and  a  is  a  constant,  and  where  x  =  0  ,  and  S  =  SQ 
so  that  aUho  =  S0  ;  then  the  equation  can  be  rewritten 


S 


S„e 


-ax 


where  Sa  =  slope  at  the  head  of  the  reach.  Shulits  (1936)  found  that  this 
equation  gave  an  approximate  fit  for  the  profile  of  the  Colorado  River  between 
Boulder  Dam  and  the  Arizona- Sonora  Boundary.  If  Z  is  the  elevation  of  any 
point  on  the  profile  above  a  given  datum,  the  following  equation  results: 


*  Shulits  (1936)  used  a  instead  of  a 
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zc  -  z  =  iy  (i  -  e_ax) 


This  equation  fits  the  profiles  of  the  Middle  Rhine,  Maas,  Mur,  and  Enns 
Rivers  in  Europe. 

48.  Sternberg  (1875)  had  arrived  at  a  similar  formula  by  combining  his 
abrasion  law  and  the  theoretical  fact  that  the  mean  velocity  just  capable  of 
moving  a  particle  is  proportional  to  the  one-sixth  power  of  its  weight.  This 
yielded  a  relation  between  mean  velocity,  particle  size,  and  distance  traveled 
by  the  grain.  Substitution  into  the  Chezy  formula  resulted  in 

ZD  -  Z  =  2  (1  -  e-“x/2)  (20) 


This  is  also  claimed  to  approximate  the  profiles  for  the  European  rivers. 

49.  Putzinger  (1919)  had  assumed  that  the  reduction  in  particle  size 
was  proportional  to  the  work  done  by  the  tractive  force  in  the  distance  trav¬ 
eled.  His  final  slope  equation  was 


ZQ  -  Z  =  K,  (1  -  e'Kx) 


(21) 


where  Kj  and  K  were  constants.  Shulits  (1936),  therefore,  concluded  with 
reasonable  certainty  that  the  slope  of  a  river  could  be  expressed  by  the  gen¬ 
eral  law  in  Equation  21.  An  investigation  of  bed  load  led  Shulits  to  modify 
the  Schoklitsch  bed- load  formula  and  combine  it  with  the  abrasion  and  slope 
formulas.  The  final  result  was 


KQ 

f  1 

g ( 4/3 )qx 

[  e“x/3  j 

(22) 


where 

G  =  bed  load,  lb/sec 
Q  =  discharge,  cfs 

b  =  width,  ft,  of  a  rectangular  channel  having  depth  d 


28 


Shulits  concluded  that  if  the  discharge  and  the  width  terms  increase  equally 
downstream,  there  would  be  no  change  in  bed  load;  if  the  discharge  term 
increases  more  rapidly  than  the  width  term,  there  would  be  an  increase  in  bed 
load  and  a  change  in  profile,  probably  degradation,  and  an  increase  in  slope; 
and  should  the  width  term  increase  more  rapidly,  there  would  be  a  decrease  in 
load  and  a  corresponding  change  in  regime. 

Lacey  (1946) 

50.  Lacey  (1946)  attempted  to  advance  the  theory  of  silt  transported  by 
"stable"  channels*  in  alluvium.  When  testing  the  earlier  regime  relations 
using  hydraulic  data,  he  discovered  two  curious  anomalies: 

a.  Despite  the  great  improbability  of  natural  alluvial  chan¬ 
nels  ever  attaining  regime,  alluvial  river  data  and  boulder 
torrent  data  (i.e.,  data  from  steeply  graded  streams)  con¬ 
firmed  the  regime  equations  in  many  respects. 

b.  Although  much  of  the  data  that  had  been  collected  over  a 
period  of  years  from  stable  channels  in  the  Punjab  dis¬ 
played  general  agreement,  certain  observations  from  the 
channels  thought  to  be  stable  proved  to  be  "discordant" 

(i.e.,  conflicting)**  for  no  readily  ascertainable  reason. 

51.  By  applying  the  principles  of  corresponding  speeds  and  dynamic 
similarity  and  by  introducing  the  terminal  velocities  of  silt  particles  into 
his  physical  theory,  Lacey  (1946)  succeeded  in  regrouping  the  fundamental 
variables  of  velocity,  depth,  and  water  surface  slope  to  obtain  two  new  param¬ 
eters  and  to  derive  a  "normal  equation"  applicable  to  all  active  alluvial 
channels  in  steady  flow.  This  "normal  equation"  would  be  able  to  implicitly 
account  for  the  degree  of  aggravation  and  scouring  and  for  departure  from 
"regime  silt  charge"!  in  the  velocity,  hydraulic  radius,  and  water  surface 


*  By  Lacey's  (1946)  definition,  "such  channels  neither  silt  nor  scour." 

Engineers  use  the  term  "in  regime"  (sometimes,  in  regimen)  to  describe  this 
condition,  while  fluvial  geologists  prefer  to  use  "graded"  or  "poised"  to 
describe  a  "stable"  channel. 

**  Lacey  (1946)  pointed  out  in  his  paper  that  much  of  the  "discordant"  data 
could  be  reconciled  when  variables  such  as  "silt  charge"  (concentration) 
and  "shock"  (the  effect  of  roughening  the  sides  and  introducing  irregu¬ 
larities  into  the  bed  of  a  given  channel  with  a  given  discharge,  sediment, 
concentration,  and  bed-material  gradation)  were  taken  into  account. 

+  Lacey  (1946)  defined  "regime  silt  charge"  as  the  minimum  transported  silt 
load  consistent  with  "full  activity  of  the  bed." 
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slope  adopted  by  the  channel,  and  in  the  actual  particle-size  gradation  of  the 
strearabed  (Lacey  1946) . 

52.  Streams  used  to  develop  common  relations  should  have  similar  char¬ 
acteristics.  Because  the  strict  definition  of  dynamic  similitude  cannot  be 
applied  to  nonrigid  boundary  streams  that  transport  not  only  water  but  also 
sediment,  Lacey  (1946)  modified  the  definition.  All  that  can  be  expected  of 
similar  alluvial  channels,  with  respect  to  scale,  is  that  they  should  be  geo¬ 
metrically  similar  in  plan  and  in  cross  section  or  "geometrically  comparable." 

53.  Lacey  (1946)  stated  that 

applying  the  theory  of  corresponding  speeds  to  any  two  alluvial 
channels  of  assumed  infinite  width  and  constant  depth,  and 
employing  the  usual  notation,  the  familiar  expression 


T'  =  L'V " 

T77  I77? 


is  obtained.  This  postulates  that,  in  respect  of  longitudinal 
flow  (L.  L' ,  .  .  .)  ,  the  particles  of  silt  and  water  traverse 

equal  longitudinal  scalar  distances  in  equal  scalar  times 
(T,T' .  .  .)  .  In  addition  to  this  requirement,  it  must  be  ensured 

that  vertical  distances  are  also  traversed  in  equal  scalar  times. 
Particles  of  silt  of  the  same  grade  and  density  fall  through  still 
water  with  the  same  constant  terminal  velocity  v  ,  and  this  can 
be  determined  by  applying  the  modified  Stokes'  law.  This  veloc¬ 
ity,  therefore,  when  silt  of  the  same  grade  is  employed,  is  common 
to  model  and  prototype.  The  author  regards  this  velocity  as  fund¬ 
amental  and  adopts  it  as  a  proportional  to  the  vertical  velocities 
of  his  models.  Hence* 


T'  =  L'V"  _  d'v" 
T"  ~  L"v'  ~  d"v' 


Employing  the  theorem  that  horizontal  planes  in  the  prototype 
remain  horizontal  in  the  model,  and  denoting  the  slopes  by  S' 
and  S"  ,  respectively, 


V"v' 

FT7 


*  Lacey  (1946)  used  D  instead  of  d  . 
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and  finally 


(VS)  =  Kv 


(23) 


where  K  is  a  constant.  Thus,  it  follows  that  in  all  similar 
alluvial  channels,  the  terminal  velocities  of  the  silt  particles 
are  directly  proportional  to  the  product  (VS)  . 

54.  Lacey  (1946)  believed  that  for  a  known  temperature  and  a  constant 
sediment  concentration,  the  mean  velocity  F0  was  a  function  of  the  average 
depth  d  ,  and  the  par  tide  -  size  diameter  D  : 

V0=f{dc,D)  <2Z0 


There  has  been  a  wide  spectrum  of  opinion  related  to  the  evaluation  of  this 
function.  The  most  common  evaluation  is  of  the  form 


Vn  =  Kd" 


Lacey  (1946)  obtained  an  exponent  of  0.5  in  evaluating  this  equation.**  He 
stated  that  for  a  constant  gradation  and  concentration,  mean  velocity  varies 
as  the  square  root  of  the  depth.  Although  the  theory  had  not  been  proved, 
Lacey  accepted  the  dynamic  hypothesis  that  there  is  a  fundamental  equation  for 
the  velocity  of  sediment  transport,  which  is  identical  to  the  equation  for  the 
velocity  of  wave  propagation  when  the  gradation  and  concentration  remain 
constant.  Following  this  theory  and  that  of  Equation  23  yields: 


=  K 


where  K  is  a  constant  (Lacey  1946) . 


*  Lacey  (1946)  used  F"(  )  ,  d  ,  and  m  instead  of  f(  )  ,  dc  ,  and  D  , 
respectively . 

**  Lacey  (1946)  used  n  instead  of  m  . 
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55.  Lacey's  first  parameter  (VS*)  is  directly  proportional  to  the  term¬ 
inal  velocity  of  the  silt  particles,  whereas  the  second  d2/2S*  is  a  function 
of  the  particle  size,  the  density,  and  the  sediment  concentration.  The  first 
parameter  represents  the  forces  tending  to  restore  the  particles  to  the  bed. 
The  second  parameter  represents  the  forces  tending  to  propel  the  particle  for¬ 
ward.  Lacey  stated  that  these  forces  should  be  correlated,  and  he  developed 
the  normal  equation  in  the  form: 


(VS*)  =  f(dc3/2S*) 


He  then  expressed  the  normal  equation  in  the  exponential  form: 


(VS*)  =  lc(dc1/2S*)n 


For  channel  flow,  Lacey  (1946)  substituted  the  hydraulic  mean  depth  d  for 
the  constant  depth  dc  employed  in  his  physical  theory  (in  which  for  constant 
depth,  the  values  of  dc  and  d  are  identical  for  wide  channels),  writing 
the  normal  equation  as 


(VS*)  =  £(d1/2S‘)n 


Lacey  believed  that  his  physical  theory  had  justified  the  form  of  his  equation 
and,  therefore,  evaluated  it  empirically.  The  result  was 


(VS*)  =  1.6(d1/2S*f/3 


(25) 


Lacey  (1946)  called  Equation  25  the  "normal  alluvial  equation"  and  designated 
the  graphical  relation  between  (VS*)  and  d1/zS*  as  the  normal  diagram, 
stating, 

There  is  no  doubt  of  the  general  trend  of  the  diagram,  and  no  sta¬ 
tistical  analysis  is  necessary  in  order  to  obtain  a  sufficiently 
accurate  constant  and  power.  It  must  be  remembered  that  the 
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observations  vary  enormously  in  weight,  and  even  if  they  all  car¬ 
ried  the  same  weight,  a  statistical  relationship,  treating  F  , 
d  ,  and  S  as  independent  variables,  would  be  objectionable  on 
physical  grounds.  Such  a  course  of  procedure  would  be  equivalent 
to  the  abandonment  of,  say,  the  Reynolds  number,  in  the  analysis 
of  data  obtained  from  pipes. 

The  normal  alluvial  equation  developed  by  Lacey  (1946)  applied  only  to  sub- 
critical  flow,  which  yields  upper  limits  of  253.0  and  43.5  for  (VS*)  and 
( d1/2S *)  ,  respectively. 

56.  From  a  study  of  the  normal  diagram,  Lacey  (1946)  concluded  that  the 
new  equation  applies  to  all  active  channels  in  alluvium  and  is  irrespective  of 
the  precise  degree  of  scouring  or  silting,  or  variation  in  silt  charge.  All 
these  variables  are  implicit  in  the  actual  values  of  V  ,  d  ,  and  S  adopted 
by  the  channel  when  associated  with  the  actual  gradation  of  the  bed  material, 
which  must  be  sampled  and  measured  at  the  same  time  as  the  other  hydraulic 
elements . 

57.  Lacey  (1946)  maintained  that  the  "full  activity  of  the  bed"  (see 
footnote,  paragraph  51),  although  random,  must  nevertheless  be  such  that  any 
reduction  would  lead  to  partial  immobility  and,  therefore,  quasi-rigidity  of 
the  channel.  With  regime  silt  charge  defined,  Lacey  stated  that  for  every 
regime  channel  transporting  a  regime  silt  charge,  there  is.  for  a  fixed  diam¬ 
eter  of  the  bed  material,  a  fixed  value  of  the  product  (VS*)  .  For  normal 
nonregime  channels,  the  actual  bed  material  was  not  in  conformity  with  the 
product  (FS*)  .  Therefore,  this  conformity  of  nonconformity  allows  "regime 
channels"  to  be  distinguished  from  "nonregime  channels." 

58.  Data  published  by  the  Lahore  Research  Institute  (1937)  confirmed 
Lacey's  "physical  theory"  and  the  normal  alluvial  equation  (Equation  25). 

Lacey  pointed  out  that  no  suspended  sediment  samples  had  been  taken  in  1937, 
and  no  method  or  device  had  been  developed  to  measure  the  "bed  silt  charge" 
(bed  load).  He  then  made  a  logarithmic  plot  of  (FS*)  against  particle-size 
diameter.  Within  the  entire  range  of  data,  the  modified  Stokes'  Law  showed 
that  the  terminal  velocity  varied  directly  with  the  particle  diameter.  From 
this  relation,  he  derived  the  following  equation: 


D  =  0.73  (VS*) 


(26) 


If  D  denotes  the  diameter  in  inches,  the  equation  becomes 
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D  =  0.03  ( VS *) 


(27) 


This  is  strictly  a  regime  equation,  and  normal  channels  scouring  or  with  a 
heavy  sediment  concentration  would  demand  a  lower  value  for  the  constant,  such 
as  0.025. 

59.  Equation  24  could  then  be  evaluated  for  the  velocity  of  bed  sedi¬ 
ment,  which  should  include  the  sediment  concentration.  Therefore, 


V  =  f(d,  D ,  C) 


where 

C  =  concentration  or  charge* 

Substituting  terminal  velocity  v  as  a  measure  of  gradation  yields 

V  =  f(d,  v,  C)  (28) 

From  Equation  25,  Lacey  derived  the  following  equation: 

V  =  1.42  (yS*)1/Ad1/2  (29) 

Lacey  believed  Equation  29  to  be  the  correct  fundamental  equation.  The 
expression  was  useful  because  (VS*)  was  easy  to  assess  for  any  system  of 
channels.  Substituting  the  terminal  velocity  for  the  product  (V5S*)  yields 
in  terms  of  the  physical  theory, 

V=l Cv1/4d1/2  (3°) 

*  Lacey  (1946)  used  s  instead  of  C  . 
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Lacey  (1946)  noted  that  Equation  30  was  applicable  only  to  a  regime  channel, 
with  a  "regime  silt  charge"  (see  footnote,  paragraph  51)  and  should  he  appli¬ 
cable  throughout  the  entire  range  of  alluvial  channel  bed  material. 

60.  Lacey  (1946)  could  not  improve  on  the  correlation, 

P  =  2.66  Q0-505 


by  the  introduction  of  the  particle  diameters  as  a  variable.  The  original 
data  presented  by  Lacey  suggested  that  the  coefficient  is  a  function  of  the 
gradation  because  the  Madras  observations  in  fine  silt  gave  a  higher  coeffi¬ 
cient  than  Kennedy's  (1895)  observations  for  coarser  material.  Lacey  attrib¬ 
uted  the  scatter  mainly  to  the  fact  that  all  the  channels  were  artificial  and 
excavated  in  soil  of  some  tenacity  and  that  some  of  the  ratios  of  bed  width  to 
depth  first  adopted  were  arbitrary  and  could  not  be  corrected.  Lacey  (1946), 
therefore,  concluded  that  the  sediment  concentration  was  clearly  of  impor¬ 
tance.  He  stated, 

It  has  frequently  been  stated  that  existing  equations  for  the 
velocity  of  flow  in  alluvial  silt  transporting  channels  rest  on  no 
sound  physical  basis.  The  reason  is  not  very  far  to  seek.  The 
problem  is  far  removed  from  that  of  flow  in  pipes,  for  example, 
and  no  laboratory  is  large  enough  to  contain  it.  It  would  appear 
that  the  academic  worker  must  abandon  Lis  laboratory  temporarily 
and  seek  a  wider  field  of  activity.  The  laboratory  worker,  how¬ 
ever,  can  greatly  assist  a  final  solution  by  carrying  out  ad  hoc 
experiments  designed  to  prove  the  truth  or  falsity  of  statements 
such  as  the  author  has  made  in  this  paper. 


Inelis  (1947.  1948) 

61.  Sir  Claude  Cavendish  Inglis,  like  others  who  preceded  him,  believed 
that  alluvial  channels  have  "a  never-ceasing  tendency  to  attain  equilibrium,  . 

.  .  a  condition  in  which  the  load  is  carried  with  the  minimum  expenditure  of 

energy  for  the  existing  conditions"  (Inglis  1947).  This  theory  holds  true, 
whether  flow  conditions  in  the  channel  closely  approximate  equilibrium  or  vary 
widely,  as  in  the  case  of  meandering  channels.  Small  adjustments  to  attain 
equilibrium  are  brought  about  by  changes  in  the  material  exposed  on  the  bed. 
Such  changes  can  occur  either  by  deposition  "where  there  is  an  excess  charge" 
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or  by  scour  "where  Che  charge  is  in  defect"  (i.e.,  where  concentration  is  low) 
(Inglis  1947). 

62.  Inglis  (1947)  stated  that  the  most  efficient  channels  carry  their 
load  with  the  least  expenditure  of  energy.  He  concluded  that  the  most  effi¬ 
cient  shape  for  a  rigid  boundary  channel  is  a  semicircle,  whereas  Lacey  (1930) 
had  suggested  that  for  mobile  bed  channels  the  most  efficient  section  would 
approximate  a  "semiellipse." 

63.  In  an  ideal  channel,  conditions  of  discharge,  concentration,  and 
temperature  remain  constant:  thus,  true  equilibrium  is  achieved.  In  nature, 
however,  the  discharge,  concentration,  and  temperature  of  the  water  all  vary, 
such  that  final  equilibrium  is  never  fully  attained.  Inglis  (1947)  suggested 
that  the  nearest  approach  to  equilibrium  is  a  year-to-year  consistency  for 
similar  discharge  conditions,  which  he  called  an  "annual  stability."  But  even 
this  "annual  stability"  varies  from  year  to  year,  owing  to  variations  of  con¬ 
centration  resulting  from  abnormal  conditions  and  changes  in  the  river.  For 
the  purpose  of  simplification,  Inglis  (1947)  assumed  that  there  was  a  dominant 
discharge  and  an  associated  charge  and  gradient  to  which  a  channel  returns 
annually.  At  this  dominant  discharge,  equilibrium  was  most  closely  approached 
and  the  tendency  to  change  the  least  (Inglis  1947).  Inglis  modified  Lacey’s 
regime  equations  to  meet  dominant  discharge  conditions  in  rivers  and  to  allow 
for  a  "super- regime  charge"  of  material  in  movement.  In  this  way,  he  hoped  to 
visualize  dominant  flow  conditions  and  to  consider  the  effects  of  variations 
from  dominant  discharge. 

64.  Lacey's  (1930)  regime  formulas,  rewritten  by  Inglis  in  terms  of 
discharge  and  "Lacey  silt  factors"  (or,  more  appropriately,  "sand  factors"), 
were : 


P  =  8/3  Q 1/2 


(31) 


R  =  0.4725 


(32) 
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0.000547  f1/6f^2 
- qT7S - 


(33)* 


According  to  Inglis  (1947),  Equation  32  shows  that,  with  the  same  discharge 
and  concentration,  the  depth  varies  inversely  as  the  cube  root  of  the  grada¬ 
tion  of  the  bed  material  and  this  relation  holds  equally  well  in  straight  or 
meandering  channels.  He  also  stated  that  in  Equation  33,  with  a  constant 
"sand  factor"  (i.e.,  silt  factor),  the  slope  required  varies  inversely  as  the 
one-sixth  power  of  the  discharge;  however,  for  constant  slope,  the  same  factor 
would  vary  as  the  one-tenth  power  of  the  discharge  (Inglis  1947).  Therefore, 
Inglis  (1947)  concluded  that  for  any  discharge  less  than  the  dominant  dis¬ 
charge,  "the  slope  will  be  in  defect  for  the  dominant  charge  so  that  material 
will  deposit,  and,  as  a  result,  there  will  be  a  subdominant  charge  of  material 
in  movement."  On  the  other  hand,  when  the  discharge  exceeds  the  dominant  dis¬ 
charge,  the  energy  available  exceeds  that  required  to  carry  the  dominant 
charge,  and  that  excess  energy  is  utilized  in  increasing  the  charge  by  scour¬ 
ing  the  bed,  eroding  the  banks,  and  developing  meanders  (Inglis  1947). 

65.  Inglis  (1947)  also  agreed  with  the  general  theories  that  had  been 
published  to  that  time  stating, 

It  is  desirable  to  emphasize  the  fact  that  the  behavior  of  a 
river  throughout  its  length  depends  on  the  material  washed  into 
it,  chiefly  near  its  head  and  from  tributaries.  Where  the  charge 
is  heavy  or  the  material  coarse,  as  is  generally  the  case  near  the 
hills,  the  slope  is  steep  and  the  channel  is  wide  and  shallow; 
where  the  material  is  fine,  farther  down  the  river,  the  channel  is 
relatively  deep  and  narrow. 

66.  Professor  C.  M.  White  presented  a  set  of  dimensional  formulas  of 
the  type  (Inglis  1948) . 


bgx/ 5 
Q  2/5 


a, 


Q  f’  [DCc'g) 

D  5/2g  1/2  J  yZ/5 


1/3  TV 


(34)* ** 


where 

V  =  velocity  of  water 

Qs  =  "quantity  of  solids  conveyed"  (suspended- sediment  load) 


*  At  regime  f  and  fRS  are  numerically  equal  (Inglis  1947). 

**  Inglis  (1948)  used  m  instead  of  D  ,  C  instead  of  Qs  ,  and  n3  ,  n2 
and  n3  ,  instead  of  ,  m2  ,  and  m3  ,  respectively. 
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Inglis  (1947)  then  worked  on  the  assumption  that  the  charge  C  =  Qs/Q  and  the 
rate  of  fall  of  particles  v  were  linked  but  that  charge  was  not  necessarily 
linked  with  grade.  Therefore,  he  modified  White's  formulas  as  follows  (Inglis 
1947)  : 


bC5  -  «! 

r  1 

Q 

m\ 

Cv 

q2/5  1 

D*/2gl/2 

k.  4 

(vg)l/3 

(35)* 


where 

De  —  effective  diameter  of  bed  sediment 
v  =  kinematic  viscosity  of  water 

Inglis  (1948)  believed  that  th^  justification  for  linking  C  with  v  and 
hence  ug  was  twofold: 

a.  An  investigation  carried  out  at  Poona,  India,  by  Inglis  showed 
that  the  rate  of  deposition  of  material  in  water  was  constant 
as  long  as  the  product  Cv  remained  constant,  as  long  as  C 
varied  inversely  as  v  ,  where  the  material  was  quartz  sand, 
with  D  being  greater  than  about  0.15  mm  and  less  thgn^about 
0.5  mm.  For  grains  coarser  than  about  0.7  mm,  v  a  Du' 

b.  Experiments  initiated  by  Inglis  at  Poona  when  he  was  Director 
of  the  Indian  Waterways  Experiment  Station  had  shown  that  when 
the  gradation  of  material  was  kept  constant  and  the  concentra¬ 
tion  varied,  the  exponents  approximated  closely  to  those  in 
Inglis'  equations. 

67.  Inglis  (1948)  stated  that  the  dimensional  components  of  his  equa¬ 
tions  would  become  either  zero  or  infinite  when  the  charge  became  zero  or  if 
the  bed  slope  and  v  became  infinitely  small.  These  conditions  do  not  occur 
in  practice  because  the  state  of  flow  changes  from  rough  turbulent  to  smooth 
turbulent,  causing  a  change  of  law. 

68.  Inglis  (1948)  compared  the  effect  of  charge  in  the  formulas 
developed  by  Lacey  with  his  own: 

Lacey  Inglis 


V  =  0.7937  Q1/6f I/3  (36)  V  =  a3  Q1/$(De  Cv)1/lz  (39)** 


*  Inglis  (1948)  used  X  and  Vs  instead  of  C  and  v. 

**  Inglis  (1948)  used  m  instead  of  D8  . 
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(40) 


^  =5.65  Q1/6fi/3  (37) 


b  Q1/6(CV)7/12 

~3  7  ^5/18|/7/36£)5/12 


1  1,828  Q1/6 

5  f573 


(38) 


]_  gi/18^5/3601/6 

S  =  a5_ (DeCv)5/12 


(41) 


where 

a3 ,  a7,  a5  =  coefficients 

In  Lacey's  formulas  (Equations  36  and  37),  V  varies  from  P/R  only  in  the 
coefficient.  In  Inglis'  charge  formulas  (Equations  39-41),  the  charge  C  has 
only  a  slight  effect  on  V  ,  but  has  a  very  great  effect  on  b/d  and,  hence, 
PR  .  The  charge  also  has  a  marked  effect  on  slope  so  that  where  charge  is 
ignored,  wide  divergences  are  found  even  in  regime  channels. 

Mackin  (1948) 

69.  J.  Hoover  Mackin,  a  geologist,  defined  the  term  "graded  stream"  as 
"one  which,  over  a  period  of  years,  slope  is  delicately  adjusted  to  provide 
with  available  discharge  and  the  prevailing  channel  characteristics,  just  the 
velocity  required  for  transportation  of  all  of  the  load  supplied  from  the 
drainage  basin"  (Mackin  1948) .  He  pointed  out  that  although  the  slope  usually 
decreases  in  a  downvalley  direction,  the  discharge,  channel  characteristics, 
and  load  do  not  vary  systematically  along  the  stream;  therefore,  the  profile 
of  a  graded  stream  is  not  a  simple  mathematical  curve.  He  stated  that  cor- 
rasive  power  of  a  stream  and  bedrock  resistance  to  corrasion  (i.e.,  mechanical 
erosion  performed  by  moving  agents,  such  as  running  water,  wind,  and  glacial 
ice)  determine  the  slope  of  the  ungraded  stream,*  but  have  no  effect  on  the 
profile  of  a  graded  stream.  The  profile  of  an  aggrading  stream**  differs  in 


*  Mackin  (1948)  defines  "ungraded"  as  out  of  equilibrium  during  a 

transitional  period  when  it  is  adjusting  itself  to  changes  in  the  drainage 
basin. 

**  Mackin  (1948)  defined  an  "aggrading  stream"  as  "upbuilding  approximately 
at  grade"  or  increase  in  the  elevation  of  the  streambed  while  maintaining 
the  same  slope  . 
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form  from  that  of  a  graded  stream,  chiefly  because  of  the  difference  in  rate 
of  downvalley  decrease  in  particle  size  of  the  sediment  load.  Another  differ¬ 
ence  between  the  two  is  that  the  aggrading  profile  is  asymptotic  with  respect 
to  horizontal  line  passing  through  the  base  level;  however,  the  graded  stream 
is  not  (Mackin  1948)  . 

70.  Mackin  (1948)  used  the  term  "stream"  to  include  specific  reaches  of 
concern.  Many  rivers  have  both  graded  and  upgraded  reaches.  By  using  the 
expression  "over  a  period  of  years,"  Mackin  ruled  out  both  short-termed 
(including  seasonal)  fluctuations  and  the  exceedingly  slow  changes  that 
accompany  the  progress  of  the  erosion  cycle.  "Load"  and  "discharge"  are  not 
necessarily  the  most  important  factors  controlling  slope,  but  they  are  the 
only  factors  that  are,  in  origin,  independent  of  the  stream.  Slope  stands 
alone  because  it  is  the  only  factor  in  the  equilibrium  that  is  automatically 
adjusted  by  the  stream  itself  in  such  a  direction  as  to  accommodate  changes  in 
external  controls  that  call  for  changes  in  velocity  (Mackin  1948)  . 

71.  According  to  Mackin  (1948),  a  graded  stream  responds  to  change  in 
accordance  with  Le  Chatelier's  general  law: 

If  a  stress  is  brought  to  bear  on  a  system  in  equilibrium,  a 
reaction  occurs,  displacing  the  equilibrium  in  a  direction  that 
tends  to  absorb  the  effect  of  the  stress. 

He  pointed  out  that  there  are  two  possible  approaches  to  the  study  of  streams 

as  agents  of  transportation: 

a.  In  terms  of  relations  between  slope,  discharge,  channel  form, 
and  the  size  of  grains  comprising  the  load. 

b.  In  terms  of  energy  transformations. 

He  also  believed  that  the  two  should  not  be  combined. 

72.  "The  energy  of  a  stream  between  any  two  points  is  proportional  to 
the  product  of  the  mass  (of  the  water-sediment  mixture)  and  the  total  fall 
between  the  two  points"  (Mackin  1948).  This  "total  energy,"  which  increases 
with  an  increase  in  discharge  or  slope,  can  be  dissipated  by  any  of  the 
following : 

a.  Friction  along  the  wetted  perimeter  of  the  channel  (external 
friction  loss) . 

b.  Friction  between  the  diverse  threads  of  the  turbulent  current 
(internal  friction  losses). 

c.  Transportation  of  load. 
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External  and  internal  friction  losses  occur  whether  or  not  the  stream  is 
engaged  in  transportation  of  load,  and  these  losses  increase  with  increased 
channel  roughness,  with  irregularity  in  the  alignment  of  the  channel,  and  with 
an  increase  in  wetted  perimeter.  Roughness,  alignment,  and  cross-sectional 
form  are  referred  to  as  "channel  characteristics"  and  are  determined  by  the 
"hydraulic  efficiency"  of  the  channel.  Based  on  earlier  work  by  Gilbert 
(1877)  and  Rubey  (1933),  Mackin  (1948)  emphasized  that  the  share  of  the  total 
energy  used  in  the  transportation  of  load  could  be  as  small  as  3  percent.  The 
fact  that  the  energy  dissipated  in  internal  and  external  friction  is  over¬ 
whelmingly  greater  than  that  consumed  in  transportation  means  that  if  total 
energy  determined  by  slope  and  discharge  remain  the  same,  then  relatively 
slight  changes  in  the  channel  characteristics  cause  very  marked  changes  in  the 
transporting  power. 

73.  Most  researchers  have  developed  equations  to  describe  the  condi¬ 
tions  of  a  graded  stream,  but  Mackin  (1948)  believed  that  this  approach  was 
inadequate  because  equations  are  transposable .  As  set  up  in  an  equation,  load 
is  a  function  of  velocity.  Mackin  (1948)  stated  that  most  engineers  agree 
that  velocity  controls  or  determines  the  load  carried  by  a  stream.  However, 
over  a  period  of  years,  the  load  supplied  to  a  stream  is  actually  dependent, 
not  on  the  velocity  of  the  stream,  but  rather  on  the  lithology,  relief,  vege¬ 
tative  cover,  and  erosional  processes  present  in  its  drainage  basin.  A  graded 
stream  will  maintain  the  particular  slope  that  will  provide  just  the  velocity 
required  to  transport  all  of  the  supplied  load.  In  this  sense,  velocity  is 
determined  by,  or  adjusted  to,  the  load.  Therefore,  Mackin  (1948)  stated  that 
"in  the  graded  stream,  load  is  a  cause,  and  velocity  is  an  effect;  this  rela¬ 
tionship  is  not  transposable." 

74.  He  noted  that  several  observations  had  shown  large  graded  streams 
usually  have  lower  slopes  than  smaller  graded  streams,  the  reason  being  that 
with  increased  channel  size,  there  is  usually  an  increase  in  cross-sectional 
area  relative  to  the  wetted  perimeter  and  a  corresponding  relative  decrease  in 
friction.  The  result  is  that  large  streams  commonly  have  higher  velocities 
than  smaller  streams  with  the  same  slope.  The  load  of  a  graded  stream  can 
increase  in  a  downvalley  direction  relative  to  its  discharge  either  because  of 
subsurface  loss  of  water  (i.e.,  to  the  banks  and  bed)  or  as  a  result  of  the 
entry  of  heavily  loaded  tributaries  vMackin  1948) . 

75.  Mackin  (1948)  believed  graded  profiles  decreased  in  a  downvalley 
direction  between  tributary  confluences  chiefly  because  of  a  downvalley 
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decrease  in  gradation  of  load  due  to  processes  within  the  stream.  He  dis¬ 
agreed  with  the  common  theory  that  the  downstream  decrease  in  the  gradation  of 
load  is  due  to  the  downvalley  decrease  in  slope,  stating  that  "downvalley 
decrease  in  caliber  of  load  is  an  important  cause  of  downvalley  decrease  in 
the  slope  of  the  graded  profile"  and  that  "in  the  graded  stream,  the  decrease 
in  caliber  is  due  primarily  to  attritional  comminution  of  particles  comprising 
the  load"  (Mackin  1948)  .  The  downvalley  decrease  in  gradation  of  load  due  to 
attrition  is  more  systematic  than  any  other  factor  affecting  the  slope  of  a 
graded  stream.  However,  the  gradation  of  load  does  not  vary  systematically  in 
graded  streams  joined  by  tributaries  nor  in  graded  streams  in  which  the  min- 
eralogic  composition  of  the  load  differs  greatly  in  resistance  to  attrition. 

76.  Like  slope,  the  other  channel  characteristics  of  a  graded  stream 
are  developed  by  the  stream  itself.  Slope  and  other  channel  characteristics 
vary  from  reach  to  reach,  and  any  change  in  external  controls  usually  results 
in  changes  in  these  variables.  Therefore,  Mackin  (1948)  revised  the 
engineer's  definition  of  a  graded  stream  from  one  "in  which  slope  is 
delicately  adjusted,"  to  one  "in  which  slope  and  (other)  channel  characteris¬ 
tics  are  delicately  adjusted"  (Mackin  1948). 

77.  Griffith  (1927)  and  Lacey  (1930)  had  suggested  that  natural  streams 
with  a  heavy  bed  load  tend  to  flow  in  broad,  shallow  channels.  Mackin  (1948) 
stated  that  the  two  most  important  reasons  for  this  shape  are: 

a.  A  large  bed  load  requires  high  bed  velocity  and  widening  by 
bank  erosion  that  must  continue  until  velocity  at  the  banks  is 
reduced  to  the  point  where  the  resistance  to  erosion  of  the 
bank- forming  materials  equals  the  erosive  force  applied  to 
them. 

b.  Shoaling  by  deposition  will  accompany  widening  of  the  narrow 
channel  by  erosion  because  the  particles  of  the  bed  load  tend 
to  lodge,  and  move,  and  lodge  again,  the  velocity  required  each 
time  to  set  them  in  motion  being  greater  than  that  required  to 
keep  them  in  motion,  and  because  a  higher  velocity  is  required 
to  set  in  motion  a  particle  of  the  bed  than  one  of  the  sloping 

banks . 

The  cross-sectional  shape  most  efficient  for  transport  varies  with  slope,  with 
amount  and  gradation  of  the  load,  and  especially  with  the  proportions  of  the 
total  load  carried  in  suspension  and  moved  along  the  bed  (Mackin  1948)  . 
Cross-sectional  shape  is  also  dependent  upon  the  resistance  of  the  banks  to 
erosion;  therefore,  for  any  set  of  slope-debris  "charge"  (concentration)  fac¬ 
tors,  there  will  be  one  critical  degree  of  erodibility  of  the  bank-forming 
materials,  such  that  a  channel  that  was  originally  too  narrow  will  quickly 
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develop  a  cross  section  with  depth-width  relations  that  provide  "the  maximum 
silt-carrying  capacity"  (Mackin  1948).  In  less  erodible  bank-forming 
materials,  the  stream  could  develop  the  maximum  efficiency  section  over  a 
period  or  time  for  the  gi' an  slope -debris  "oherge’1  factors.  However,  if  the 
bank-forming  materials  are  more  erodible  than  the  critical  degree,  the  stream 
will  adopt  and  maintain  a  section  that  is  wider  and  shallower  than  the  ideal 
transportation  section.  Mackin  (1948),  therefore,  concluded 

A  statement  to  the  effect  that  a  flow  of  water  charged  with 
debris  must  necessarily  develop  for  itself  that  form  of  section 
that  will  give  it  a  maximum  debris -carrying  capacity,  is  an 
invalid  generalization  because  it  ignores  erodibility  of  the 
banks . 

He  also  noted  that  stability  is  no  guarantee  of  maximum  efficiency  (Mackin 
1948)  . 

78.  In  affirming  the  validity  of  Le  Chatelier’s  general  law  (para¬ 
graph  71),  Mackin  (1948)  described  the  manner  in  which  a  system  in  equilibrium 
would  react  to  certain  stress  brought  to  bear  on  it.  He  commented  on  increase 
in  load,  decrease  in  load,  change  in  discharge,  rise  of  base  level,  and  lower¬ 
ing  of  base  level,  as  follows  (Mackin  1948): 

a.  Increase  in  load.  "A  once-graded  stream  responds  to  an 
increase  in  load  primarily  by  steepening  its  declivity  below 
the  point  of  influx.  The  steepening  is  accomplished  by  depo¬ 
sition  of  part  of  the  excess  load  in  the  channel  at  the  point 
of  the  influx  with  a  consequent  upbuilding  of  the  channel  at 
that  point  and  the  formation  of  a  steepened  part  immediately 
below.  Steepening  of  any  segment  permits  increased  transport 
of  load  through  that  segment  to  the  next  segment  which  is  in 
turn  the  site  of  deposition  and  steepening.  Thus  the  effect  of 
an  increase  in  load  is  registered  by  the  downvalley  movement  of 
a  wave  deposition,  large  or  small  depending  on  the  rate  and 
manner  of  addition  of  the  load,  and  deposition  must  continue 
throughout  the  stream  below  the  point  of  influx  until  the  slope 
is  everywhere  adjusted  to  the  transport  of  all  the  debris 
delivered  to  it." 

b.  Decrease  in  load.  "A  decrease  in  load  may  be  thought  of  as 
occurring  at  any  point  in  the  stream.  It  may  be  stated  for 
the  time  being  that  the  stream  simply  makes  up  for  defi¬ 
ciency  in  the  load  supplied  from  above  by  picking  up  addi¬ 
tional  load  from  its  channel  floor.  The  net  result  is 
downcutting,  with  a  consequent  lowering  in  declivity  down- 
valley  from  the  point  where  the  change  occurred.  Downcut¬ 
ting  must  continue  until  the  profile  is  reduced  to  that 
slope  which  will  provide  just  the  velocity  required  to 
transport  the  reduced  load.” 

c.  Change  in  discharge.  "If  a  segment  of  a  graded  stream 
receives  all  or  most  of  its  load  at  the  upper  end,  changes 


in  discharge  call  for  readjustments  in  the  slope  in  much 
the  same  way  as  changes  in  load.  A  decrease  in  discharge 
requires  an  increase  in  declivity  because  the  load, 
remaining  the  same,  must  move  faster  through  a  smaller 
cross  section,  and  because,  as  indicaced  earner,  decrease 
in  the  cross-sectional  area  of  the  channel  involves  a  rela¬ 
tive  increase  in  frictional  retardation  of  flow  and,  hence, 
a  decrease  in  velocity.  The  stream  affected  by  a  decrease 
in  discharge,  being  unable  to  transport  all  the  load  sup¬ 
plied  to  it  on  its  former  slope,  deposits  some  of  the  load 
and  thereby  steepens  the  slope,  the  process  continuing 
until  the  reduced  stream  is  able,  by  reason  of  increased 
velocity  on  the  steepened  slope,  to  transport  all  the  load 
shed  into  it.  The  opposite  adjustment  occurs  in  the  case 
of  increased  discharge." 

d.  Rise  of  base  level.  "A  rise  of  base  level  is  equivalent  to 
the  rising  of  a  barrier  across  the  path  of  the  graded 
stream.  Each  unit  of  increase  in  the  height  of  the  barrier 
trends  to  flatten  the  declivity  immediately  upstream.  The 
stream,  unable  because  of  decreased  declivity  to  carry  all 
of  the  load  through  the  flattened  segment,  deposits  in  the 
segment,  thus  increasing  the  declivity  and  transferring  the 
flattening  upvalley.  Continuation  of  the  process  results 
in  upstream  propagation  of  a  wave  or,  better,  of  an  infi¬ 
nite  number  of  small  waves  of  deposition.  If  the  barrier 
is  raised  slowly  the  stream  may  maintain  itself  in  approxi¬ 
mate  adjustment  during  the  process;  its  rate  of  aggradation 
is  determined  by  the  rate  at  which  the  barrier  is  elevated. 
If  the  barrier  is  raised  rapidly  or  instantaneously  a  lake 
is  formed;  the  distal  part  of  the  delta  is  then  the  'flat¬ 
tened  part'  of  the  profile,  and  the  rate  of  aggradation  is 
determined  by  the  rate  of  delta  building  into  the  lake.  In 
either  case  the  successive  profiles  developed  during  the 
period  of  readjustment  will  differ  markedly  in  form  from 
the  original  profile  and  from  the  eventual  completely  read¬ 
justed  profile.  The  final  readjusted  profile  will  tend 
toward  parallelism  with  the  original  profile,  differing  in 
this  respect  from  the  cases  treated  above.  But,  because  of 
secondary  effects  of  aggradation  .  .  .  precise  parallelism 
will  usually  not  be  achieved;  the  only  generalization  that 
can  be  made  is  that  the  new  profile  will  be  everywhere 
adjusted  to  the  new  prevailing  conditions." 

e.  Lowering  of  base  level.  "Lowering  in  base  level  is,  inso¬ 
far  as  the  response  of  the  stream  is  concerned,  essentially 
the  same  as  the  lowering  of  a  barrier  in  its  path.  Each 
small  lowering  of  the  control  point  steepens  the  gradient 
immediately  upstream.  Accelerated  velocity  in  the  stee¬ 
pened  portion  results  in  downcutting,  and  the  steepening  is 
propagated  upvalley.  Downcutting  must  continue  until  the 
slope  is  again  completely  adjusted  to  supply  just  the  velo¬ 
city  required  to  transport  all  of  the  debris  shed  into  the 
stream;  as  in  the  last  case,  and  with  the  same  qualifica¬ 
tion,  the  final  readjusted  profile  will  tend  to  parallel 
the  original  profile." 
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79.  These  changes  can  be  classified  into  two  categories:  (a)  changes 
in  downvalley  controls  and  (b)  changes  in  upvalley  controls.  In  general,  a 
change  in  upvalley  controls  (load  and  discharge)  results  in  a  change  in  slope, 
whereas  a  change  in  downvalley  controls  (base  level)  will  cause  a  change  in 
the  level  of  the  profile,  but  not  necessarily  in  the  slope.  Mackin  (1948), 
therefore,  had  suggested  many  new  ideas  in  the  study  of  graded  or  regime 
rivers.  He  developed  many  general  statements  but  avoided  manipulation  of  the 
data  into  equation  form. 


Blench  (1951) 

80.  Thomas  Blench  stated  that  a  river  is  "in  regime"  if  its  mean  mea¬ 
surable  behavior  during  a  certain  time  interval  does  not  differ  significantly 
from  its  mean  measurable  behavior  during  comparable  times  before  or  after  the 
given  interval.  He  pointed  out  that  regime  theory  originated  on  the  great 
canal  systems  of  India  and  Pakistan  where  discharges  are  controlled  and  mean¬ 
ders  prevented;  however,  this  theory  is  not  generally  applicable  to  a  "regime  - 
type  river"*  (Blench  1951) . 

81.  Blench  (1951)  believed  that  the  present  generalized  theory  would 
not  exist  without  the  Lacey  theory.  Lacey  (1930)  had  used  a  single  silt  fac¬ 
tor  f  (Equation  9) ,  thereby  effectively  obtaining  results  for  a  physically 
possible  set  of  channels  whose  relative  importance  of  side  action  to  bed 
action  had  been  averaged  out  to  some  arbitrary  value.  This  permitted  him  to 
work  in  terms  of  wetter  perimeter  P  and  hydraulic  radius  R  .  Blench  (1951) 
found  the  exponents  of  Lacey's  regime  formulas  to  be  consistent  with  the  laws 
of  rigid  boundary  hydraulics.  As  width  goes  to  infinity,  depth  approaches  the 
hydraulic  radius  and  width  approaches  the  wetted  perimeter;  therefore,  because 
most  rivers  are  wide,  Blench  replaced  P  with  b  and  R  with  d  and,  leav¬ 
ing  Lacey's  exponents  the  same,  obtained: 


=  4 


(42) 


*  Defined  by  Blench  (1951)  as  "one  that  has  formed  a  major  part  of  every 
cross  section  from  material  that  has  been  transported  or  could  be 
transported  at  some  stage  of  flow," 
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where  fb  =  a  bed  factor  related  to  the  nature  of  the  sediment  load,  which  is 
constant  as  long  as  the  factors  relevant  to  the  bed  load  are  constant.  He 
believed  this  formula  to  be  exact  when  the  discharge  and  bed  load  were  con- 
o'-auL  aiiJ  tegiiue  conditions  had  been  attained.  Therefore,  Equation  42  pro¬ 
vides  the  mean  for  defining  a  dynamically  satisfying  bed  factor  (Blench  1951). 

82.  Blench' s  second  generalized  regime  equation  was 


"5" 


(43) 


where  fs  -  side  factor  that  measures  the  tractive  force.  The  width  term  b 
in  Equation  43  multiplied  by  average  depth  equals  cross-sectionai  area. 
Equation  43,  like  Equation  42,  can  be  taken  as  exact  for  "perfect  regime 
conditions,"  but  is  also  practically  applicable  to  real  regime  conditions 
(Blench  1951)  .  Blench  developed  this  formula  while  attempting  to  determine 
the  reason  for  the  universal  agreement  of  the  exponent  (0.5)  in  the  Lacey 
(iyiO)  formula  (Equation  10) 

P  =  2.67  Q1/2 


83.  There  is,  however,  significant  disagreement  over  the  coefficient 
2.67,  which  Blench  (1951)  contended  should  be  replaced  by  a  variable  dependent 
upon  the  nature  of  the  bed  sediment.  His  line  of  reasoning  was  as  follows 
(Blench  1951): 

Universal  acceptance  of  the  exponent  1/2  (in  Equation  [44]*)  based 
on  data  covering  a  very  wide  range,  suggests  that  the  expression 
PV0R  has  a  physical  significance.  Replacing  P  and  R  by  b 
and  d  ,  respectively,  the  expression  can  be  written  as 

Vo  /d 

Vj7b 


in  which  the  numerator  has  physical  significance.  Therefore, 
there  should  be  physical  significance  in  the  denominator.  Since 
Vl/d  is  a  force  per  unit  mass,  V\/b  is  probably  in  some  way  a 
measure  of  tractive  force  on  the  sides.  PzvVl/b  measures  the 
square  of  the  tractive  force  intensity  on  the  sides  (except  for  a 
multiplying  constant) 


*  Equation  and  figure  numbers  enclosed  in  brackets  within  quoted  material 
refer  to  this  report. 
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provided  that  there  is  a  laminar  film  at  the  sides  (that  is, 
provided  the  sides  are  technically  "smooth"). 

84.  Blench,  therefore,  concluded  that  the  sides  of  the  channels 
observed  must  have  been  smooth  and  that  the  Lacey  formula  represents  nothing 
more  than  an  average  of  the  relative  importance  of  bed  to  sides  as  measured  by 
the  ratio  fb/fs  ■  This  meant  that  all  Lacey  channels  were  regime  channels, 
but  not  all  regime  channels  were  Lacey  channels.  Logically,  but  not  practi¬ 
cally  desirable,  fs  should  be  replaced  by  vfs  (Blench  1951). 

85.  The  third  equation  of  Blench' s  generalized  regime  theory  was 


vl  m  c,  VcP1'* 

gdSc  u 


(45)* 


where 

Sc  =  channel  slope  expressed  as  a  fraction 
C  —  nondimens ional  constant 


V„b 


Reynolds  number  in  terms  of  width  (see  Equation  15) 


Algebraic  manipulation  of  Equations  43,  44,  and  45  yields 


d  =  3 


SQ 

N  fl 


(46) 


b  = 


(47) 


fS/65l/12c,-1/6 

~crg 

„l/4 


(48) 


According  to  Blench  (1951),  the  value  of  C'g/vl/l '  ,  with  v  averaging  10~5  is 
2,080.  The  designer  must,  however,  determine  the  value  of  fb  that  nature 
will  impose  on  his  canal  system,  and  there  is  no  scientific  formula  to  assist 
in  making  the  decision.  There  is  a  rough  empirical  formula,  adapted  from  an 
equation  proposed  by  Lacey  (1936)  (who  was  careful  to  emphasize  its  empiricism) 

*  Blench  (1951)  used  C  instead  of  C'  ,  and  S  instead  of  Sc  . 
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4  =  2^ 


(49) 


in  which  Db  is  tilt;  mean  diameter  of  the  bed  material  (in  millimetres)  to  be 
expected  in  the  canal.  Scientifically,  Equation  49  is  unsatisfactory.  Prac¬ 
tically,  however,  it  fits  many  cases  quite  well,  probably  because  the  quanti¬ 
ties  in  movement  are  very  minute  in  most  canal  applications  (Blench  1951)  . 

86.  When  applying  regime  formulas  derived  on  canals  to  rivers,  an 
allowance  for  "disturbing  effects"  (i.e.,  turbulence)  and  variable  discharges 
must  be  made  "in  terms  of  engineering  common  sense  and  the  formulas  regarded 
with  suspicion  until  supported  by  results"  (Blench  1951).  Some  apparently 
satisfactory  applications  were  found  to  be  (Blench  1951) : 

a.  Determination  of  width  between  incised  banks. 

b.  Determination  of  scour  between  bridge  piers. 

c.  Determination  of  scour  downstream  from  piers,  along  groins,  and 
at  spur  heads. 

d.  Estimation  of  aggradation  upstream  from  reservoirs. 

e.  Estimation  of  degradation  downstream  from  reservoirs. 

f.  Estimation  of  dredging. 

g.  Determination  of  model  scales. 

Leopold  and  Maddock  (1953) 

87.  Geomorphology,  the  branch  of  geology  dealing  with  landforms  and 
their  genesis  and  history,  had  been  classically  treated  almost  exclusively  in 
a  qualitative  manner.  By  the  early  1950s,  the  US  Geological  Survey  had  vast 
amounts  of  streamflow  data  from  all  over  the  United  States,  which  consisted  of 
concurrent  measurements  of  mean  velocity,  width,  shape,  area  of  cross  section, 
and  discharge.  Suspended- sediment  data  covering  a  10-year  period  were  also 
available.  Luna  B.  Leopold  and  Thomas  Maddock,  Jr.,  (1953)  noted  that  these 
data  were  not  collected  for  the  study  of  fluvial  morphology  and,  as  a  result, 
were  incomplete  in  terms  of  water  surface  slope,  bed  roughness,  and  bed  load. 

88.  Leopold  and  Maddock  (1953)  investigated  changes  in  velocity,  depth, 
and  top  width  for  an  increase  in  discharge  at  a  particular  cross  section  of  a 
river.  In  their  paper,  the  term  "depth"  referred  to  mean  depth  resulting  from 
dividing  the  cross-sectional  area  of  flow  by  the  top  width.  "Velocity" 
referred  to  a  mean  velocity  obtained  by  dividing  the  discharge  by  the 
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cross-sectional  area  of  flow.  Although  the  authors  realized  that  this  was  not 
the  most  useful  measure  of  velocity,  especially  when  dealirg  with  sediment 
transport,  it  was  the  only  measure  of  velocity  for  which  a  large  quantity  of 
data  was  available. 

89.  They  found  that  the  relation  of  width,  depth,  and  velocity  to  dis¬ 
charge  could  be  approximated  by  straight  lines  on  logarithmic  plots.  Some 
scatter  was  attributable  to  temporary  scour  and  fill  of  the  streambed  exclud¬ 
ing  any  unnatural  effects;  the  data  analyzed  indicated  that  through  ranges  of 
discharge  up  to  the  bank- full  stage.  The  relations  of  width,  depth,  and 
velocity  to  discharge  could  be  described  by  simple  power  functions: 


d  =  azQm ’ 


(51) 


yo  -  *3<?" 


(52) 


where  m  , 
90. 


mz  ,  m3  ,  ax  ,  az  ,  and  a3  are  numerical  constants.* 

Because  width,  depth,  and  velocity  are  all  functions  of  discharge, 


Q  =  bdV 


therefore , 


Q  =  (a !<?“•)  ( azQ "•)  (a3Qra>) 


or 


*  Leopold  and  Maddock  (1953)  used  W  instead  of  5  ;  a  ,  c  ,  and  k  , 
instead  of  ax  ,  a2  ,  and  a3  ,  respectively;  and  b  ,  f  ,  and  m  instead 
of  ml  ,  mz  ,  and  m3  ,  respectively. 
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Q  -  axa2a3Q 


thus 


m1  +  m2  +  m3  =  1 . 0 


and 


a1a2az  =1.0 


The  constants,  ax  ,  a2  ,  and  s3  represent  the  intercepts  of  the  *  i  i  :*  ■ 
the  logarithmic  plots  and  are,  respectively,  values  of  width,  depth,  ar.d  :r<  .u. 
velocity  at  a  discharge  of  unity.  Because  in  actual  data,  the  discharges  are 
greater  than  unity  or,  better  expressed,  the  product  of  the  values  of  width, 
depth,  and  mean  velocity  at  a  given  discharge  must  equal  that  discharge. 

Width  and  depth  vary  widely  from  one  cross  section  to  another  at  a  given  dis¬ 
charge;  therefore,  more  work  was  necessary  before  suggesting  a  range  of  values 
for  the  coefficients  ax  ,  a2  ,  and  a3  .  The  average  values  for  the 
exponents  mx  ,  m2  ,  and  m3  ,  which  are  biased  toward  semiarid  conditions 
because  they  were  computed  for  20  river  cross  sections  on  a  variety  of  rivers 
in  the  Great  Plains  and  Southwest,  are  (Leopold  and  Maddock  1953): 

mx  =  0.26 
m2  =  0.40 
w3  =  0.34 

From  these  computations,  Leopold  and  Maddock  (1953)  noted  that  depth  increases 
faster  than  width  for  an  increase  in  discharge.  The  relative  rates  of 
increase  of  width  and  depth  are  functions  of  the  channel  shape  (Leopold  and 
Maddock  1953)  . 

91.  Next,  they  plotted  width,  depth,  and  velocity  corresponding  to  the 
mean  annual  discharge  at  cross  sections  from  all  parts  of  a  river  system  and 
showed  that  discharge  increases  in  the  downstream  direction  because  of  the 
increase  in  drainage  area.  The  depth,  width,  and  velocity  all  tend  to 
increase  as  power  functions  of  the  discharge.  The  general  alignment  of  points 
on  the  downstream  graphs  indicates  that,  in  a  given  river  basin  where  all 
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cross  sections  are  experiencing  the  same  frequency  of  discharge,  the  corre¬ 
sponding  values  of  depth,  width,  and  mean  velocity  at  different  cross  sections 
(even  on  different  tributaries)  having  the  same  discharge  tend  to  be  similar, 
regardless  of  the  location  within  the  watershed  (Leopold  and  Maddock  1953) . 

92.  Leopold  and  Maddock  (1953)  point  out  that  most  geomorphologists  are 
under  the  impression  that  the  velocity  of  a  stream  is  greater  in  the  head¬ 
waters  than  in  the  lower  reaches  because  the  river  is  steeper  in  the  upper 
reaches.  However,  Manning's  equation, 


V  =  ^  ^2/3£t/2 

n  e 


(53)* 


shows  that  the  depth  of  flow  (or  hydraulic  radius)  is  more  influential  than 
slope  on  change  in  velocity.  Leopold  and  Maddock  (1953)  concluded  that  the 
rate  of  increase  of  depth  downstream  tends  to  overcompensate  for  the  decreas¬ 
ing  slope  and  to  provide  a  net  increase  of  velocity  at  mean  annual  discharge 
in  the  downstream  direction  of  a  river.  When  they  plotted  width,  depth,  and 
velocity  against  the  mean  annual  discharge  at  various  cross  sections,  the 
result  was  a  simple  power  function.  The  exponents  and  coefficients  varied 
from  basin  to  basin,  but  the  average  exponent  values  for  the  basins  studied 
we  re  : 

m1  =  0.5 
m2  =  0.4 
m3  =  0.1 

93.  Thus,  the  changes  in  hydraulic  characteristics  downstream  in  a 
river  system  have  been  considered  using  only  the  mean  annual  discharge.  They 
then  made  comparisons  between  the  data  already  presented  and  similar  data  for 
infrequent  flows.  At  each  gaging  station,  in  addition  to  the  mean  annual  dis¬ 
charge,  they  studied  flows  that  were  equaled  or  exceeded  1,  4,  10,  30,  and 
50  percent  of  the  time.  For  each  discharge  frequency,  curves  representing  the 
change  of  hydraulic  factors  with  discharge  in  the  downstream  direction  were 
plotted.  Leopold  and  Maddock  (1953)  stated 


*  Leopold  and  Maddock  (1953)  used  1.5  for  1.49  and  substituted  d  for  R 
and  S  for  Se  in  Manning's  equation. 
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The  channel  characteristics  of  natural  rivers  are  seen  to 
constitute,  then,  an  interdependent  system  which  can  be  described 
by  a  series  of  graphs  having  simple  geometric  form.  The  geometric 
form  of  the  graphs  describing  these  interactions  suggests  the  term 
"hydraulic  geometry."  Channel  characteristics  of  a  particular 
river  system  can  be  described  in  terms  of  the  slopes  and  inter¬ 
cepts  of  the  lines  in  the  geometric  patterns.  .  .  .  Though  the 

listing  of  slopes  and  intercepts  of  these  lines  may  not  provide 
any  visual  picture  of  a  river  basin,  comparison  of  the  values  of 
these  factors  among  rivers  has  useful  aspects. 

94.  Leopold  and  Maddock  (1953)  also  studied  the  effect  of  suspended 
sediment  on  hydraulic  geometry.  By  the  early  1950's,  suspended- sediment  sam¬ 
ples  were  being  taken  on  a  regular  basis  at  a  number  of  locations  on  major 
streams.  The  samples  in  any  given  cross  section  could  then  be  analyzed  for 
concentration  and,  with  the  mean  discharge  for  that  same  day,  used  to  estimate 
the  passing  daily  suspended-sediment  load.  Plots  of  the  relation  of 
suspended- sediment  load  versus  discharge  could  then  be  made.  These  plots  were 
characterized  by  a  large  scatter  of  points;  however,  for  long  periods  of 
record,  this  relation  can  be  approximated  by  a  simple  power  function  of  the 
form  (Leopold  and  Maddock  1953)  : 


Qs  - 


(54)* 


where 

Qs  =  suspended-sediment  load,  tons  per  day 

a,  j  =  numerical  constants,  with  the  range  of  j  being  on  the  order  of 
2  to  3 

95.  They  concluded  that  typically  at  a  given  cross  section,  the 
suspended-sediment  load  should  increase  more  rapidly  than  the  discharge.  Con¬ 
sideration  of  the  physical  characteristics  of  surface  runoff  and  observations 
of  bed  scour  led  to  the  tentative  conclusion  that  the  observed  increase  in 
sediment  concentration  resulted  primarily  from  erosion  of  the  watershed  rather 
than  from  scour  of  the  bed  of  the  main  stream  in  the  reach  where  the  measure¬ 
ment  was  made.  There  were  not  enough  suspended- sediment  stations  on  any  given 
stream  to  analyze  downstream  changes  in  concentration.  Leopold  and  Maddock 
(1953)  stated  that  in  long  reaches  of  a  river  through  more  or  less  homogeneous 

*  Leopold  and  Maddock  (1953)  used  L  instead  of  Qs  and  p  instead  of  a  . 
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topography  or  lithology,  the  Input  of  various  tributaries  did  not  appear  to 
have  any  great  effect  on  the  turbidity  of  the  water  (Leopold  and  Maddock 
1953)  . 

96.  Leopold  and  Maddock  (1953)  contended  that  certain  aspects  of  the 
effect  of  suspended  sediment  on  channel  shape  were  best  exemplified  by  con¬ 
sideration  of  the  condition  of  constant  discharge.  For  each  station,  width 
was  plotted  against  suspended  load  on  a  semilogarithmic  graph  (Figure  3)  (with 
suspended  load  being  plotted  on  the  logarithmic  axis) ,  and  at  each  point  the 
corresponding  velocity  was  entered.  Velocity  isopleths  slanted  steeply  down¬ 
ward  to  the  right,  while  depth  isopleths  slanted  gently  downward  to  the  right. 
The  limitations  of  the  graph  were: 

a.  The  velocity  field  was  only  roughly  defined  by  values  in  the 
available  data. 

b.  Only  the  suspended- sediment  load  was  considered. 

c.  There  was  a  large  scatter  of  points  on  any  average  curve  show¬ 
ing  relation  of  suspended  sediment  to  discharge. 

Figure  3  shows  that  there  was  an  increased  capacity  for  suspended  load  at  con¬ 
stant  discharge  by  decreasing  width  at  a  constant  velocity  or  by  increasing 
velocity  at  constant  width. 

97.  The  slope  j  of  the  curve  relating  suspended  load  to  discharge 
(Equation  54)  depends  on  the  rate  of  increase  of  velocity  with  discharge 
defined  by  the  exponent  m3  (Equation  52),  and  because  mx  +  m2  +  m3  =  1  , 


Figure  3.  Average  relation  of  suspended-sediment  load  to  channel  width  for 
a  constant  discharge  of  500  cfs.  Velocity  V  isopleths  are  in  feet  per 
second,  while  depth  d  isopleths  are  in  feet  (adapted  from  Leopold  and 

Maddock  1953) 
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when  m1  and  m3  are  specified,  then  mz  is  also  specified.  Therefore,  if 
the  ratio  of  m3  to  m2  is  known  and  the  value  of  m1  is  known,  then  the 
values  of  m3  and  m2  are  also  known.  From  these  relations,  Leopold  and 
Maddock  (1953)  developed  Figure  4. 

From  Figure  [4],  a  number  of  useful  concepts  may  be  inferred. 

The  first  concerns  the  interaction  of  velocity,  depth,  and  the 
suspendedsediment  load.  For  any  given  rate  of  increase  of  width 
with  dischargethat  is,  for  any  given  value  of  m1  ,  the  diagram 
indicates  that  the  value  of  j  increases  with  an  increase  in  the 
m3  to  m2  ratio.  That  is,  the  rate  of  increase  of  suspended- 
sediment  load  with  discharge  is  a  function  of  the  ratio: 


rate  of  increase  of  velocity  with  discharge 
rate  of  increase  of  depth  with  discharge 


Slopes  of  the  lines  representing  depth,  width,  and  velocity 
respectively,  in  relation  to  discharge  are,  according  to  the 
implications  contained  in  Figure  [4] ,  functions  of  the  slope  of 
the  line  representing  the  relation  of  suspended  sediment  to 


Figure  4.  Relation  of  the  rates  of  increase  of  width, 
depth,  and  velocity  and  suspended-sediment  load 
expressed  by  values  of  m1  ,  m3/m2  ,  and  j  (adapted 
from  Leopold  and  Maddock  1953) 
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discharge.  For  a  given  value  of  b  ,  the  steeper  the  slope  of  the 
suspended  sediment  discharge  line,  the  greater  will  be  the  slope 
of  the  velocity-discharge  line  (Leopold  and  Maddock  1953)  . 

98.  According  to  Leopold  and  Maddock  (1953),  the  suspended- sediment 
concentration  should  decrease  in  the  downstream  direction.  A  value  of 

j  =  1.0  would  mean  that  suspended  load  and  discharge  are  increasing  at  the 
same  rate.  Therefore,  the  value  of  j  for  natural  rivers  should  be  somewhat 
less  than  1.0.  Using  the  mean  values  for  mx  ,  m2  ,  and  m3  ,  Figure  4 
correctly  yields  j  =  0.8  .  Therefore,  Figure  4  appears  to  apply  to  average 
conditions  both  at  a  station  and  downstream.  By  their  analysis,  Leopold  and 
Maddock  (1953)  stated,  "For  a  given  width,  and  at  a  given  discharge,  an 
increase  in  suspended- sediment  load  requires  an  increase  in  velocity  and  a 
reduction  in  depth."  From  the  study  of  graphs  developed  by  Gilbert  (1914) 
from  flume  data,  they  'Concluded, 

At  constant  discharge,  an  increased  velocity  at  constant  width 
is  associated  with  an  increase  of  both  suspended  load  and  bed  load 
in  transport.  At  constant  velocity  and  discharge ,  an  increase  in 
width  is  associated  with  a  decrease  of  suspended  load  and  an 
increase  in  bed  load  in  transport  (Leopold  and  Maddock  1953)  . 

99.  Past  researchers  had  placed  too  much  emphasis  on  bed  load,  and  the 
role  of  suspended  load  in  determining  the  observed  characteristics  of  natural 
streams  had  been  underestimated  (Leopold  and  Maddock  1953) .  From  analysis  of 
rivers  in  the  western  United  States  during  floods,  Leopold  and  Maddock  (1953) 
determined  that  the  changes  in  the  riverbed  occurred  simultaneously  with 
changes  in  the  rate  of  change  of  suspended-sediment  concentration. 

It  is  the  thesis  of  this  paper  that  the  observed  changes  in  the 
stream  bed  resulted  from  changes  in  the  sediment  load  brought  into 
the  measuring  reach  from  upstream.  It  is  postulated  that  the 
hydrodynamic  factors  involved  tend  to  promote  a  mutual  adjustment 
between  channel  shape  and  the  sediment  load  carried  into  the 
reach.  The  change  in  sediment  load,  which  results  in  a  change  of 
channel  shape,  involves  both  bed  load  and  suspended  load.  How¬ 
ever,  because  only  the  suspended  load  is  measured,  it  is  necessary 
to  use  the  data  on  the  suspended  fraction  of  the  load  as  an  index 
to  how  the  total  sediment  load  interacts  with  the  hydraulic  vari¬ 
ables.  That  the  suspended  fraction  is  a  meaningful  index  is  dem¬ 
onstrated  by  the  fact  that  the  relations  between  the  channel  shape 
factors  and  suspended  load,  which  are  derived  from  measurements  of 
a  number  of  different  rivers,  appear  to  apply  in  principle  to 
channel  changes  at  a  given  station  during  an  individual  flood. 
Specifically,  with  no  change  of  channel  width,  a  decrease  of  sus¬ 
pended  load  at  a  given  discharge  was  accompanied  by  an  increase  in 
depth  by  bed  scour  that  resulted  in  a  decrease  in  velocity 
Figure  [4],  The  decrease  in  velocity  provides  the  adjustment  of 
capacity  for  carrying  the  load  of  the  particular  size  distribution 
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supplied  by  the  watershed.  In  response  to  a  decrease  in  load,  the 
channel  shape  became  adjusted  through  scour  to  the  lower  capacity 
required  for  quasi -equilibrium . 

100.  Changes  in  channel  shape  occur  in  response  to  changes  in  sediment 
load  brought  into  the  reach  (Leopold  and  Maddock  1953)  .  This  postulate 
required  proof  that  the  observed  change  in  suspended- sediment  load  was  not  the 
result  of  the  observed  change  of  velocity  in  the  reach,  but  rather  the  cause. 
If  the  high  suspended-sediment  concentrations  resulted  from  the  scouring 
action  of  high  velocities,  the  implication  was  that  high  velocity  in  a  given 
reach  scours  the  channel  in  that  reach.  The  increase  in  sediment  in  transport 
resulting  from  the  local  bed  scour  should  then  account  for  the  observed 
increase  in  sediment  concentration.  Under  such  an  assumption,  increasing 
velocity  should  be  associated  with  bed  scour  and  decreasing  velocity  with 
aggradation.  Leopold  and  Maddock  (1953)  concluded  that  the  scour  and  fill  of 
the  bed  of  the  main  stem  of  an  alluvial  river  during  a  flood  appeared  to  be  a 
rapid  adjustment  of  channel  shape  in  response  to  a  varying  sediment  load. 
Vanoni  (1941)  had  shown  that  an  increase  in  suspended  load  tended  to  decrease 
channel  roughness  and  cause  an  increase  in  velocity.  Buckley  (1922)  had  also 
found  similar  results.  The  effect  of  suspended- sediment  concentration  on 
channel  roughness  was  also  discussed  by  Thomas  (1946),  who  concluded  that  an 
increase  in  concentration  resulted  in  lower  values  of  Mannings's  n  . 

101.  Any  large  changes  in  suspended- sediment  concentration  that 
occurred  with  changes  in  discharge  in  a  given  stream  should  be  reflected  by 
substantial  changes  in  roughness.  Leopold  and  Maddock  (1953)  had  found  a 
change  in  velocity-depth  relations  that  coincided  with  a  change  in  suspended 
load-discharge  relations  at  a  discharge  to  be  well  below  the  peak  on  the  Colo¬ 
rado  and  San  Juan  Rivers.  Such  local  adjustments  in  velocity-depth  relations 
resulted  primarily  from  changes  in  roughness  associated  with  changes  in  sedi¬ 
ment  load,  rather  than  from  changes  in  the  slope  of  the  water  surface.  They 
found  that  a  change  in  roughness  caused  by  a  change  in  suspended- sediment  load 
was  also  evident  350  miles  downstream  from  Hoover  Dam  at  Yuma,  AZ.  Since  con¬ 
struction,  a  quasi-equilibrium  has  been  established  at  Yuma  by  change  in  chan¬ 
nel  characteristics  and  not  by  changes  in  slope.  Changes  in  slope  that  have 
occurred  in  upstream  reaches  near  dams  could,  in  time,  affect  the  Yuma  reach. 

102.  Comparing  regime  canals  to  natural  rivers,  Leopold  and  Maddock 
(1953)  stated: 

The  comparison  of  regime  canals  and  natural  rivers  is  valid 
only  in  the  downstream  direction.  Canals  carrying  water  high  in 
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sediment  concentration  would  not  normally  be  expected  to  maintain 
the  regime  condition  except  at  the  discharge  for  which  they  were 
designed.  The  lack  of  a  varying  discharge  at  any  particular  canal 
cross  section  precludes  a  proper  comparison  with  the  at-a-station 
relation  of  rivers.  Failure  to  recognize  this  would  lead  to  fal¬ 
lacious  reasoning,  as  in  a  discussion  by  Stevens  (1937) .  (See 
also  Lane  (1935),  in  which  he  compared  the  width-discharge  formula 
of  a  particular  river  gaging  station  with  the  Lacy  (1930)  width- 
discharge  formula  for  canals.) 

Another  fact  that  demonstrates  the  analogy  between  regime 
canals  and  rivers  in  the  downstream  direction  lies  in  a  considera¬ 
tion  of  slope.  In  the  design  of  canals,  it  has  usually  been  found 
necessary  to  provide  greater  slope  in  the  small  laterals  than  in 
the  large  trunk  canals.  Thus,  canals  in  regime  have  a  longitudi¬ 
nal  profile  convex  upward  because  the  flow  is  reversed  as  compared 
with  natural  rivers;  that  is,  the  discharge  in  a  canal  system 
decreases  downstream.  This  is  analogous  except  for  the  direction 
of  flow  to  the  condition  in  natural  rivers,  which  are  steeper  in 
the  headwaters  where  the  discharge  is  small  than  downstream  where 
the  flow  is  great. 

Concurrent  changes  in  m3  ,  m2  ,  and  j  during  individual  floods  (Figure  4) 
define,  at  least  in  general,  the  mutual  adjustments  that  occur  locally  when 
suspended- sediment  concentration  changes.  Because  Figure  4  describes  rela¬ 
tions  both  local  and  downstream  for  average  values  of  the  factors  observed  in 
rivers  and  approximately  fits  regime  canals,  any  point  in  the  graph  was  con¬ 
sidered  by  Leopold  and  Maddock  (1953)  to  represent  the  concurrent  values  of 
m3  ,  m2  ,  and  j  required  for  equilibrium.  Leopold  and  Maddock,  however, 
noted  that  the  exact  form  and  positions  of  the  lines  in  Figure  4  must  be  con¬ 
sidered  subject  to  adjustment  as  more  data  become  available. 

103.  Leopold  and  Maddock  (1953)  also  found  that  graded  reaches  of  a 
river  showed  a  width-depth-velocity-discharge  relation  similar  to  that  of 
reaches  not  known  to  be  graded.  They  stated  that  this  similarity  indicated 
that  the  same  factors  operating  to  maintain  equilibrium  in  a  graded  river  are 
also  sufficiently  active  in  nongraded  reaches  to  provide  specific  consistent 
patterns  in  the  hydraulic  variables,  both  locally  and  in  a  downstream  direc¬ 
tion.  Stream-gaging  data,  however,  were  not  available  in  sufficient  quantity, 
in  a  close  geographic  network,  or  over  a  long  enough  period  to  allow  differen¬ 
tiation  between  graded  and  nongraded  reaches  by  hydraulic  geometry  alone, 
though  such  differentiation  could  theoretically  be  possible. 

104.  From  their  research,  Leopold  and  Maddock  (1953)  concluded  the  fol¬ 
lowing  with  regard  to  the  longitudinal  profile  of  a  river: 

Channels  of  different  shapes  have  characteristic  distribu¬ 
tions  of  velocity  in  the  cross  section.  With  the  same  discharge, 
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a  very  shallow  and  wide  channel  tends  to  have  a  greater  average 
shear  in  the  vertical  plane  than  in  the  horizontal.  Conversely  a 
very  deep  and  narrow  channel  tends  to  have  greater  average  shear 
on  the  bank  than  on  the  bed.  The  relative  erodibility  of  bed  and 
bank  allows  the  channel  shape  to  develop  a  velocity  distribution 
which  is  in  approximate  equilibrium  with  the  particular  erodibil¬ 
ity  characteristics.  In  this  manner  the  increase  of  width  dis¬ 
charge  (value  of  mx)  probably  is  determined.  Why  this  exponent 
(n^)  should  so  consistently  be  nearly  equal  to  0.5  in  the  down¬ 
stream  direction  for  widely  different  rivers  is  not  known  and  con¬ 
stitutes  an  important  unsolved  problem. 

Different  concavities  of  longitudinal  profile  are  associated 
with  different  relations  of  velocity  and  depth  to  discharge,  as 
required  by  the  hydraulics  of  open  channels.  However,  a  particu¬ 
lar  rate  of  increase  of  both  velocity  and  depth  downstream  is 
necessary  for  maintenance  of  approximate  equilibrium  in  a  channel, 
inasmuch  as  the  drainage  area  produces  sediment  and  water  in  a 
characteristic  manner. 

Roughness  varies  with  sediment  load  and  with  particle  size 
but  the  interaction  of  opposing  tendencies  results  in  a  small 
range  of  values  of  roughness  downstream  in  most  rivers.  The  sus¬ 
pended  load  and  its  change  downstream  characteristic  of  natural 
rivers  require  a  particular  rate  of  increase  of  velocity  and  depth 
downstream.  Under  the  conditions  of  a  nearly  constant  roughness, 
to  provide  the  required  velocity- depth  relations  slope  must  gener¬ 
ally  decrease  downstream;  it  is  for  this  reason  that  the  longitu¬ 
dinal  profile  of  nearly  all  natural  streams  carrying  sediment  is 
concave  to  the  sky. 

On  the  basis  of  their  work,  Leopold  and  Maddock  (1953)  agreed  with  Mackin 
(1948)  ,  stating  that  the  wording  in  the  definition  of  regime  should  be  changed 
from  "in  which  slope  is  delicately  adjusted"  to  one  in  which  slope  and  all 
other  channel  characteristics  are  delicately  adjusted. 


Lane  (1955) 


105.  Fluvial  morphology,  as  developed  by  geomorphologists,  had  been 
largely  a  qualitative  or  descriptive  science  that  suffered  somewhat  from  a 
lack  of  quantitative  relations.  Lane  (1955)  wanted  to  present  geomorphology 
from  an  engineering  viewpoint  because  the  concept  of  equilibrium  in  streams  as 
developed  by  the  geomorphologist  was  very  useful  to  engineers. 

106.  Lane  (1955)  was  a  proponent  of  the  energy  concept;  that  is,  when  a 
stream  is  capable  of  carrying  more  sediment  than  is  supplied  to  it,  bank  ero¬ 
sion  and  scour  occur.  Scour  increases  the  slope  of  the  tributary  streams, 
thus  allowing  the  streams  to  introduce  more  material  to  the  main  stem.  The 
more  a  stream  scours  into  the  surrounding  land,  the  greater  the  available 
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load.  Lane  stated  that  most  streams  eventually  reach  either  the  ocean  or  some 
other  base  level  as  they  cut  down  their  beds,  decreasing  their  slopes  and 
diminishing  their  ability  to  transport  the  sediment  brought  to  them.  As  the 
amount  of  material  brought  down  to  the  stream  increases  and  the  stream's  abil¬ 
ity  to  carry  it  decreases,  the  rate  at  which  sediment  is  brought  to  the  stream 
becomes  equal  to  the  rate  at  which  the  stream  carries  it  away,  and  equilibrium 
is  attained.  Because  of  the  rapid  variations  of  flow  and  sediment  supply,  a 
condition  of  equilibrium  in  which  the  load  exactly  equals  the  capacity  of  the 
stream  rarely  exists,  except  for  a  short  duration.  For  most  practical  engi¬ 
neering  purposes,  however,  the  beds  of  a  large  number  of  streams  are  in  equi¬ 
librium  (Lane  1955) . 

107.  Lane  (1955)  defined  equilibrium  as  follows: 

For  engineering  purposes,  a  section  of  a  stream  may,  there¬ 
fore,  be  said  to  be  in  equilibrium  if,  although  it  may  continually 
fluctuate  between  aggradation  and  degradation  over  a  long  period 
of  years  in  terms  of  human  history,  the  net  amount  of  change  is 
not  sufficiently  large  to  be  detected  by  quantitative 
measurements . 

Most  alluvial  streams  not  affected  by  the  works  of  man  are  in  equilibrium  for 
most  practical  engineering  purposes  (Lane  1955).  Lane  also  contended  that  if 
equilibrium  existed  in  a  reach  and  one  of  the  conditions  were  changed,  then 
over  time  the  conditions  of  equilibrium  would  change  throughout  the  entire 
reach.  In  long  reaches,  there  is  frequently  more  than  one  change  at  a  time; 
however,  one  of  the  changes  commonly  overshadows  the  effects  of  the  other 
changes  with  regard  to  the  new  equilibrium  conditions. 

108.  In  most  cases,  there  was  not  a  sufficient  amount  of  data  to  make  a 
quantitative  analysis  of  equilibrium  change;  therefore,  Lane  (1955)  developed 
the  following  expressions  as  a  qualitative  guide: 

QbD  ~  QSb  <55)* 


Here  Qb  is  the  quantity  of  sediment,  D  is  the  particle 
diameter  or  size  of  the  sediment,  Q  is  the  water  discharge,  and 
Sb  is  the  slope  of  the  stream.  This  is  an  equation  of  equilib¬ 
rium  and  if  any  of  the  four  variables  is  altered,  it  indicates  the 
changes  which  are  necessary  in  one  or  more  of  the  others  to 
restore  equilibrium.  For  example,  if  a  stream  with  its  sediment 
load  is  flowing  in  a  condition  of  equilibr ium  and  its  sediment 


*  Lane  (1955)  used  Qs  instead  of  Qb  ,  d  instead  of  D  ,  Qw  instead  of 
Q  ,  and  S  instead  of  Sb  . 
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load  is  decreased,  equilibrium  can  be  restored  if  the  water  dis¬ 
charge  or  the  slope  is  decreased  sufficiently  or  if  the  diameter 
of  the  sediment  is  increased  the  proper  amount.  This  equation  is 
not  an  exact  mathematical  equation  as  it  will  not  give  the  quanti¬ 
tative  values  of  the  variables  involved  which  bring  about  equilib¬ 
rium,  but  it  is  helpful  to  indicate  qualitatively  the  changes 
which  will  take  place  in  a  stream  when  a  change  of  any  one  of  the 
variables  occurs. 

The  sediment  discharge  Qb  in  this  equation  is  the  coarser 
part  of  the  sediment  load  or  more  exactly  the  bed-material  load, 
since  this  is  the  part  of  the  load  which  largely  molds  the  bed 
formation.  In  most  cases,  the  quantity  of  the  fine  load  of  silt 
and  clay  sizes  can  change  almost  indefinitely  without  materially 
affecting  the  river  profile  (Lane  1955). 


Chien  r 1955 ^ 


109.  Until  the  mid-1950's,  the  basic  equations  as  proposed  by  Lacy 
(1946)  remained  essentially  the  same,  although  others  had  attempted  to  express 
them  in  different  forms  with  various  interpretations.  The  formulas  as  pro¬ 
posed  by  Lacey  involved  factors  that  were  based  essentially  on  experience 
alone,  making  use  of  t-he  theory  difficult  for  those  who  lacked  this  experi¬ 
ence.  Ning  Chien,  therefore,  decided  to  evaluate  these  factors  quantitatively 
in  terms  of  the  physical  conditions  imposed  on  the  channel  and  to  show  the 
possible  limitations  of  Lacey's  (1946)  regime  theory. 

110.  Inglis  (1947)  had  rewritten  Lacey's  (1946)  three  regime  formulas 
(Equations  31-33),  and  later,  Chien  (1955)  expressed  the  Lacey  formulas,  as 
follows : 


V2 

~  =  1-325 


(56) 


£l/3S2/3  =  0.0052  fRS 


(57) 


P  =  aQl/2  =  2.67  Q1/z 


where 

^VR.  ^RS 
a 


silt  factors* 
coefficient 


(58) 


*  Inglis  (1947)  had  used  f  instead  of  fTO  . 
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He  stated  that  the  above  questions  "can  be  applied  only  to  channels  which  con¬ 
tain  a  bed  of  loose  sediment  of  the  same  type  that  is  moved  along  the  bed  and 
which  are  essentially  at  equilibrium"  (C hien  1955). 

111.  Chien's  (1955)  study  was  structured  around  four  premises,  the 
first  two  concerning  "the  behavior  of  alluvial  channels,  in  general,"  and  the 
second  two  dealing  with  "Lacey's  theory,  in  particular." 

a.  "The  equilibrium  state  of  an  alluvial  channel  is  attained  by 
adjusting  the  dimensions  of  the  cross  section  and  the  slope  of 
the  channel  in  accordance  with  the  conditions  imposed  on  the 
channel  by  the  drainage  basin.  These  conditions,  which  are 
independent  of  the  channel  itself,  include  the  sediment 
inflow,  the  discharge,  and  the  gradation  and  size  of  the  mate¬ 
rial  which  composes  the  channel  bed  and  the  banks.  It  has 
long  been  recognized  that  the  sediment  inflow  is  a  primary 
factor  governing  the  stability  of  an  alluvial  channel. 

b.  "The  alluvial  channel  differs  from  the  fixed-boundary  channel 
by  the  fact  that  the  boundary  is  molded  by  the  flow  and  may 
change  with  the  flow.  The  frictional  resistance  of  an  allu¬ 
vial  channel  flow  cannot  be  described  generally  by  one  formula 
with  universal  constants,  such  as  the  Manning's  formula.  It 
tr’S*"  be  interpreted  as  a  composite  effect  of  the  following 
parts : 

(1)  Resistance  of  the  sediment  grains  which  compose  the 
channel  bed. 

(2)  Resistance  of  banks. 

(3)  Resistance  of  sand  bars  and  other  irregularities. 

Each  of  these  follows  different  and  independent  laws. 

c.  "Equations  [56-58]  are  empirical  relationships  based  on  field 
observations  only.  They  should  not  be  visualized  as  fixed 
laws  with  universal  application. 

d.  "Lacey's  regime  theory  is  the  product  of  sixty  years'  experi¬ 
ence  gathered  in  India  and  Pakistan,  and  it  serves  today  as  a 
useful  tool  in  designing  irrigation  canals  in  that  area." 

112.  Chien  (1955)  hoped  to  find  a  philosophical  basis  for  Lacey's 
regime  theory  determining  why  the  sediment  load  was  not  included.  He  believed 
that  Lacey’s  omission  of  sediment  load  as  an  explicit  variable  could  be 
explained  by  either  of  the  following  two  possibi lities : 

a.  "The  silt  factors  .  .  .  actually  include  implicitly  the  sedi¬ 

ment  load  together  with  the  bed  material  size. 

b.  "Although  the  dimensions  and  the  slope  of  an  alluvial  channel 
must  depend  on  the  sediment  inflow,  a  certain  combination  of 
these  variables  may  depend  very  little  on  the  sediment  load" 
(Chien  1955) . 


61 


113.  Chien  (1955)  used  graphical  solutions  relating  channel  depth  and 
slope  with  a  unit  discharge  and  sediment  load  taken  from  his  earlier  work 
(Chien  1954)  as  the  basis  for  his  analysis.  These  graphical  solutions  were 
constructed  according  to  the  Einstein  bed-load  function,  the  only  working  tool 
available  at  that  time  to  completely  describe  alluvial  channel  flow  and  sedi¬ 
ment  transport.  He  first  ignored  bank  friction  and  investigated  velocity, 
hydraulic  radius,  slope,  discharge,  bed-material  size,  and  silt  factors. 
Plotting  the  silt  factors  against  sediment  concentration,  he  developed  the 
following  equations  (Chien  1955) : 


LVR 


0.061 


qT 


0.715 


(59) 


■^RS  ~ 


1.18 


(60) 


where 


qT  =  total  sediment  transport  rate  (including  both  the  bed  load  and  the 
suspended  load)  per  unit  width 

q  =  discharge  per  unit  width 

The  silt  factor  fm  depends  strongly  on  the  sediment  load  and  is  the  com¬ 
bination  of  slope  and  hydraulic  radius  and  is  practically  independent  of  the 
sediment  transport  rate  (Chien  1955). 

114.  Chien  (1955)  then  stated 

In  the  wider  range  of  variables,  especially  of  the  sediment 
size,  the  trend  of  fRS  remains  the  same.  With  sediment  size  D 
as  a  variable  and  expressed  in  millimetres,  Equation  [60]  becomes 


^rs  -  2.2  D 


0.45 


„  \0 . 052 

9t[ 

.  q 


(61) 


for  qT/q  <  200  ppm.  The  change  of  fTO  with  sediment  concentra¬ 
tion  not  only  follows  a  somewhat  different  pattern,  but  also  is 
quite  erratic.  For  other  than  conditions  in  India  and  Pakistan, 
fw  depends  also  on  the  hydraulic  characteristics  of  the  channel. 
This  makes  the  general  application  of  Lacey's  regime  theory 
extremely  difficult.  It  is  not  just  the  question  of  whether  one 
has  the  proper  training  or  experience  in  estimating  the  values  of 
silt  factors,  but  rather  the  experience  gathered  in  northern  India 
and  Pakistan  may  fail  to  apply  when  the  conditions  are  substan¬ 
tially  different  from  those  existing  in  that  area. 
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He  found  that  bank  friction  had  little  effect  on  the  overall  trends  according 
to  which  the  silt  factors  vary  with  the  sediment  concentration.  He  did,  how¬ 
ever,  note  a  slight  shift  in  relative  position  of  the  fTO  and  fRS  curves 
for  small  channels  (Chien  1955) . 

115.  Chien  (1955)  noted  two  important  points: 

a.  "The  sediment  load  used  in  this  study  refers  only  to  the  bed- 
material  load,  the  material  which  has  been  found  abundant  in 
the  bed.  The  actual  sediment  concentration  of  the  canal  flow, 
including  that  of  the  wash-load,  may  be  considerably  higher. 

b.  "In  using  Einstein's  bed-load  function,  the  effect  of  bar 
resistance  is  determined  from  river  measurements.  There  is 
some  evidence  indicating  that  below  a  certain  limit,  the  bar 
resistance  is  affected  by  the  width  of  the  channel.  In  a 
narrow  channel,  the  bars  do  not  develop  to  their  full  sizes, 
due  to  the  confinement  by  the  banks.  This,  in  effect,  will 
bring  the  values  of  fTO  and  fRS  closer  together  at  small 
sediment  concentrations,  but  will  not  materially  change  the 
conclusions  derived  from  the  present  analysis." 

116.  The  findings  of  Chien' s  (1955)  study  can  be  summarized  as  follows: 

a.  "The  silt  factors  in  Lacey's  regime  theory  have  been 
correlated  with  the  conditions  imposed  on  the  channel  by  the 
watershed  according  to  the  Einstein  bed-load  function. 

b.  "For  the  channels  in  the  Indian  plain  area  from  which  the 
regime  theory  was  originally  derived,  fTO  depends  on  the 
sediment  concentration  of  the  flow  while  fRS  varies 
primarily  with  the  bed  material  size.  Using  these  relation¬ 
ships,  both  the  Einstein  bed- load  function  and  Lacey’s  regime 
theory  give  the  same  depth  and  slope  of  an  alluvial  channel  in 
equilibrium. 

c.  "The  general  application  of  Lacey's  regime  theory  to  other 
than  conditions  in  India  and  Pakistan  is  questionable  as  fra 
then  depends  also  on  the  characteristics  of  the  flow,  making 
its  selection  extremely  difficult." 

Schumm  (1960) 

117.  Lack  of  a  simple  quantitative  expression  for  the  physical  proper¬ 

ties  of  alluvium  had  hindered  the  study  of  hydraulics  and  morphology  of 
stream'  St-^nley  A.  Schumm  (1960)  pointed  out  that  studies  of  the  behavior  of 
sediment  grains  under  different  physical  environments  revealed  a  great  change 
in  sediment  character  within  the  0.05-  to  1.0-mm  range.  Schumm  (1960)  dis¬ 
cussed  the  effect  of  sediment  character  on  the  shape  of  alluvial  streams.  He 
pointed  out  "to  relate  any  cr.e  aspect  of  stream  morphology  to  one  other 
variable  .  .  .  was  an  oversimplification  of  the  problem  and  that  the  reader 
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should  be  aware  of  the  other  important  factors  in  stream  morphology"  (Schumm 
1960)  . 

118.  Schumm  chcsc  the  No.  200  sieve  (0.074  mm)  as  the  upper  limit  for 
silt  and  clay.  He  stated 

Assuming  that  the  shape  of  a  stream  channel  depends  on  the 
resistance  of  sediment  composing  the  perimeter  of  the  channel  and 
the  erosion  potential  of  stream  discharge,  then  sampling  of  bank 
and  channel  sediment  at  stable  channel  cross  sections  is  neces¬ 
sary.  Data  were  assembled  for  90  cross  sections,  most  of  which 
were  at  or  near  Geological  Survey  gaging  stations.  Generally,  the 
gaging  stations  are  located  at  stable  reaches  of  the  river.  How¬ 
ever,  they  often  are  located  at  bridges  which  might  be  assumed  to 
affect  the  shape  of  the  cross  section.  In  such  situations,  wher¬ 
ever  practical,  the  samples  were  collected  some  distance  upstream 
or  downstream  from  the  bridge.  Sometimes,  however,  the  depth  of 
water  necessitated  obtaining  samples  and  water  depth  from  the 
bridge  by  using  a  small  clam-shell  type  dredge  (Schumm  1960). 

Generally,  the  channel  and  banks  were  sampled  only  to  a  depth  of  1  in. ,  with 

10  and  20  sample  points  taken  across  the  stream;  in  some  cases,  a  ^-in.-deep 

sample  was  collected.  In  the  laboratory,  the  samples  were  subjected  to  a 

standard  grain-size  analysis.  The  samples  were  first  sieved,  and  if  they 

contained  more  than  20-percent  srlt-clay,  they  were  prepared  for  hydrometer 

analysis  (Schumm  1960). 

119.  Schumm  (1960)  described  the  character  of  the  sediment  comprising 
the  perimeter  of  each  channel  in  terms  of  its  weighted  mean  percent  silt-clay, 
which  he  calculated  as  follows: 


(sc)  (bf)  +  (sb)  (2dm) 

- 5^-7^ - 


(62)* 


where 

M  =  weighted  mean  percent  silt-clay 
sc  =  percentage  of  silt-clay  found  in  the  channel  alluvium 
bt  =  bank- full  width 

sb  =»  percentage  of  silt-clav  found  in  the  bank  alluvium 
The  shape  of  each  cross  section,  expressed  as  a  dimensionless  width-depth 
ratio,  was  plotted  against  M  .  The  result  of  this  plot  was  the  equation: 

Schumm  (I960)  used  S,  ,  W  ,  Sb  ,  and  D  instead  of  sc  ,  b£  ,  sb  , 
and  dm  ,  respectively. 
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(63)* 


^  =  255  M'108 

Neither  the  percentage  of  silt-clay  in  the  banks  nor  in  the  channel  alone 
shows  a  correlation  with  the  width-depth  ratio.  Equation  63  assumes  that 
channels  containing  little  silt-clay  are  relatively  wide  and  shallow,  whereas 
those  composed  predominantly  of  silt-clay  are  relatively  narrow  and  deep.  A 
plot  of  the  median  grain  size  against  percentage  of  silt  and  clay  in  the  chan¬ 
nel  samples  showed  no  correlation,  indicating  that  median  grain  size  has  no 
relation  to  channel  shape  for  the  range  of  channels  investigated.  Schumm 
(1960)  stated 

It  is  hazardous  to  attempt  to  extrapolate  the  relationships 
of  (width-depth  and  weighted  mean  percent  silt-clay)  as  a  general 
law  for  alluvial  channel  shape,  but  it  is  interesting  to  note  that 
Rubey's  (1952)  comparison  of  the  Illinois  and  Mississippi  Rivers 
reveals  differences  similar  to  those  noted  between  the  sections 
characterized  by  high  M  and  low  M  in  this  study.  .  .  .  The 

reason  for  the  relation  between  channel  shape  and  percentage  of 
silt  and  clay  is  found  in  studies  made  by  hydraulic  engineers. 

Lane  (1935)  states,  '.  .  .  the  greater  the  width-depth  ratio  the 
greater  will  be  the  ratio  of  the  velocity  acting  on  the  bottom  to 
that  acting  on  the  sides  of  a  channel.'  Where  a  narrow  trench  is 
cut  in  an  alluvial  valley,  the  tractive  forces  acting  on  the  chan¬ 
nel  sides  will  be  great,  causing  widening  of  the  channel.  Widen¬ 
ing  will  continue  until  the  resistance  of  the  banks  to  scour 
prevents  it.  If  the  material  in  which  the  channel  is  cut  is 
highly  cohesive  (has  a  high  percent  of  silt-clay),  the  channel 
will  be  narrow,  but  if  the  alluvium  lacks  cohesion  (has  a  small 
percent  of  silt-clay),  the  channel  will  widen  to  a  greater  extent. 

120.  Schumm  (1960)  also  noted  that  the  shape  of  the  channels  seemed  to 
be  independent  of  the  discharge,  and  the  absolute  size  of  the  channel  (i.e., 
width  and  depth)  was  related  to  mean  discharge  (Leopold  and  Maddock  1953) ,  but 
the  ratio  of  width  to  depth  was  apparently  determined  by  sediment  type  for  the 
channels  Schumm  sampled.  When  he  plotted  depth  versus  M  ,  Schumm  found  a 
large  scatter  of  points.  Because  of  the  lack  of  correlation  between  channel 
depth  and  M  Schumm  (1960)  concluded  that  perhaps  channel  width  is  more  sen¬ 
sitive  to  changes  in  M  than  to  changes  in  depth. 

121.  Schumm  (1960)  also  investigated  the  effect  of  M  on  downstream 
changes  in  channel  dimensions.  He  found  that  if  M  remained  constant 
downstream,  the  width  and  depth  of  the  channel  increased  at  a  uniform  rate 


*  Schumm  (1960)  used  F  instead  of  5f/c/m  for  the  channel  shape  factor 
(width-depth  ratio). 
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with  discharge,  excluding  the  influence  of  other  variables,  and  the  width- 
depth  ratio  remained  constant.  If  H  increased,  downstream  channel  depth 
increased  more  rapidly  than  width  (which  might  even  decrease);  therefore,  the 
width-depth  ratio  decreased.  If  M  decreased,  downstream  channel  width 
increased  more  rapidly  than  depth  (which  might  even  decrease) ,  and  width- depth 
ratio  increased. 

122.  This  study  by  Schumm  (1960)  suggests  that  M  represents  the 
resistance  to  erosion  or  general  behavior  of  sediment  in  a  stream  channel 
containing  only  small  amounts  of  gravel.  Neither  the  mean  annual  discharge 
nor  the  mean  annual  flood  affects  this  relation  significantly,  at  least  for 
the  channels  that  Schumm  sampled. 


Schumm  (1968) 

123.  Schumm  (1968)  conducted  a  study  of  the  Murrumbidgee  River  in  Aus¬ 
tralia  and  its  associated  paleochannels  in  the  Riverine  Plain 

to  obtain  additional  data  that  would  be  relevant  to  the  hypothesis 
that  the  sediment  load  largely  determines  the  shape  and  pattern  of 
alluvial  channels  and  would,  in  turn,  support  the  proposed  clas¬ 
sification  of  alluvial  channels. 

Ten  readily  accessible  representative  reaches  of  the  Murrumbidgee  River  were 
selected  for  study.  Five  of  these  reaches  were  at  or  near  existing  gaging 
stations.  At  each  location,  a  cross  section  was  surveyed,  and  bed  and  bank 
sediment  samples  were  collected. 

124.  Particle-size  determinations  were  made  by  seiving  and  by  hydro¬ 
meter  analysis.  The  percentage  of  sediment  smaller  than  0.074  mm  (i.e.,  silt- 
clay)  in  the  perimeter  can  be  expressed  by  Equation  62.  The  modern 
Murrumbidgee  River  channel  has  a  low  width-depth  ratio,  a  low  gradient,  a 
moderately  high  senuosity,  and  a  small  sediment  load. 

123.  Dury  (1965)  had  shown  that  meander  wavelength' was  a  better  crite¬ 
rion  of  bank- full  discharge  than  the  meander  belt  width.  He  subsequently 
presented  the  following  equation: 


£  =  30  Qq'5 


(64)* 


*  Schumm  (  1  968)  has  used  Qh  instead  of  Qfl  . 
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whe  re 


2  =  meander  wavelength,  ft 
Q0  =  bank-full  discharge,  cfs 

As  is  usually  the  case  with  empirical  equations,  there  is  a  scatter  of  data 

points.  Schumm  (1968)  pointed  out  that 

Dury  (1965)  recognized  that  this  scatter  may  reflect  the  influence 
of  other  variables,  such  as  sediment  type,  and  Hack  (1965)  and 
Alexandre  (1962)  suggested  that  the  sediment  transported  by 
streams  is  a  significant  factor  which  can  influence  the  dimensions 
of  the  meander  pattern. 

To  test  this  hypothesis,  Schumm  (1968)  performed  a  multiple  regression  analy¬ 
sis*  to  determine  the  relation  of  meander  wavelength,  bank- full  discharge,  and 
percentage  of  silt-clay  in  the  channel  perimeter  using  data  from  28  US  rivers, 
three  Murrumbidgee  River  sections,  two  "prior-stream"  sections,  and  an  "ances¬ 
tral  river"  section  of  the  Murrumbidgee,**  which  restilted  in  the  following 
equation : 


2 


438 


Qe  43 


(65) 


126.  Another  multiple-regression  analysis  of  the  data  collected  from 
the  rivers  in  the  United  States  and  the  Murrumbidgee  River  yielded  the  follow¬ 
ing  relation  among  meander  wavelength,  mean  annual  discharge,  and  percentage 
of  channel  silt-clay: 


2  =  1,890 


M°-  7‘ 


(66) 


where  Qm  =  mean  annual  discharge,  cfs.  The  same  set  of  data  also  yielded  the 
following  relation  (Schumm  1968): 


*  Schumm  (1968,  1969)  provided  statistical  data  for  all  variables  that  he 
i nves t igated . 

**  Schumm  (1968)  stated  that  the  Riverine  Plain  of  the  Murrumbidgee  River 
contains  traces  of  "old  aggraded  and  abandoned  river  channels."  He  termed 
the  youngest  set  of  this  as  "ancestral  river  channels,"  while  the  older  set 
has  been  referred  to  as  "prior-stream"  channels. 


6/ 


(67) 
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where  Qma  =  mean  annual  flood,  cubic  feet  per  second 

127.  To  determine  the  influence  of  discharge  and  the  type  of  sediment 
load  on  channel  morphology,  Schumm  (1968)  combined  earlier  data  from  US  rivers 
with  those  collected  at  the  Wagga  Wagga,  Narrahdera,  and  Darlington  Point 
cross  sections  on  the  Murrumbidgee  River.  A  multiple  regression  of  that  data 
yielded  the  following  relations: 


bf  =  37 


,,0.38 

'mD 


M 


(68)* 
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128.  The  same  set  of  data  was  again  used  to  determine  the  effect  of 
discharge  and  type  of  sediment  load  on  the  channel  depth.  When  these  two 
variables  were  used  in  a  multiple  regression  analysis,  the  following  equation 
was  obtained: 


dm  =  0.59  M°-3t,Q°'29 


(70) 


When  the  mean  annual  flood  was  substituted  for  the  mean  annual  discharge,  the 
following  equation  resulted: 


d  =  0.09  rt°'34Q°a*Z 


(71) 


The  same  data  showed  a  relation  between  width-depth  ratio  and  sediment  load  of 
the  form: 


*  Schumm  (1968)  used  b  instead  of  bc 
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^£  =  106M-°-78  (72)* 

4 


This  was  just  slightly  different  from  Equation  63.  The  addition  of  the  mean 
annual  discharge  and  the  mean  annual  flood  resulted  in  the  following 
relations : 


(73) 


(74) 


Schumm  (1968)  concluded  that 

although  the  relation  between  width-depth  ratio  and  M  differs 
from  that  obtained  previously,  the  general  relationship  that 
bed-load  channels  are  wide  and  shallow,  whereas  suspended- load 
channels  are  narrow  and  deep,  pertains. 

129.  Schumm  (1968)  then  introduced  the  following  equation: 


ut 


-.0.10 


56 


M 


0.74 


,0.18 


=  21 


M 


0.74 


Sb  =  59.5  W° ■  38Qr 


-0.32 


(75)** 


If  the  valley  slope  was  also  used,  the  equation  became 


,0 . 94 


Sk  =  1  ■  3 


M 


0.23o0.02 


(76) 


where  Sv  =  valley  slope.  Equation  76  demonstrated  that  on  a  given  valley 
slope,  a  straight  channel  would  begin  to  meander  if  the  mean  annual  discharge 
of  M  increased  (bed-load  decreases).  Schumm  (1963)  had  earlier  found  M 


*  Schumm  (1968)  used  F  instead  of  br/dm  . 

*  Schumm  (1968)  used  S(  instead  ot  . 
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and  sinuosity  p  (ratio  of  channel  length  jo  valley  length)  to  be  related  as 
follows : 


p  =  0.94  M 


0.25 


(77)* 


However,  the  data  in  his  1968  study  changed  the  relation  to  the  form: 

p  =  1,05  M0-18  (78) 


Based  on  his  analysis,  Schumm  (1968)  concluded  that  the  morphology  of  a  stable 
alluvial  river  channel  reflects  the  hydrologic,  climatic,  and  geologic  char¬ 
acteristics  of  the  drainage  basin. 

Schumm  (1969) 

130.  One  objective  of  Schumm' s  (1969)  work  was  to  demonstrate  that, 
with  time,  river  regulation  could  cause  a  complete  metamorphosis  of  river 
morphology.  Although  the  immediate  result  of  dam  construction  is  local  and 
generally  manifested  by  channel  degradation,  geologic  evidence  suggests  that 
the  influence  of  regulatory  and  diversion  structures  can  cause  a  complete 
change  of  river  channel  morphology  throughout  the  length  of  the  river  system 
(Schumm  1969)  . 

131.  Schumm  (1969)  stated  that  the  dimensions,  slope,  gradient,  and 
pattern  of  stable  alluvial  rivers  could  be  controlled  by  the  quantity  of  water 
and  quantity  and  type  of  sediment  moved  through  their  channels.  To  investi¬ 
gate  these  interrelations,  he  collected  data  on  the  channel  geometry  and  sedi¬ 
ment  characteristics  of  36  stable  alluvial  rivers  in  the  field  and  obtained 
hydrologic  data  from  the  records  of  nearby  gaging  stations.  The  channels  he 
selected  were  located  in  semiarid  to  subhumid  regions  of  the  Great  Plains  in 
the  United  States  and  of  the  Riverine  Plain  in  New  South  Wales,  Australia. 
These  channels  had  shown  no  progressive  adjustments  during  the  previous 

10  years  of  record. 

132.  Schumm  (1969)  found  no  relation  between  channel  morphology  and  the 
particle  size  of  the  bed  and  bank  sediments,  possibly  reflecting  the  small 
range  in  bed-material  size  of  the  channels  studied.  Although  data  on  the 


*  Schumn  (1968,  1969)  used  P  instead  of  p 
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total  sediment  load  were  available  for  only  five  cross  sections,  the  percent¬ 
age  of  silt  and  clay  (particles  less  than  0.074  mm  in  diameter)  in  the  sedi¬ 
ments  forming  the  perimeter  of  the  channels  was  found  to  be  inversely  related 
to  the  percentage  of  the  total  sediment  load  that  was  bed-material  load  at 
mean  annual  discharge  using  the  following  relation: 


M  =  21  (79) 

Qt 

where  Qt  =  percentage  of  total  sediment  load  that  is  bed  load.  Schumm  (1969) 
stated  that  M  could  be  used  as  an  index  of  the  ratio  of  bed-material  load  to 
total  sediment  load  and  that  the  percentage  of  silt  and  clay  can  be  the  same 
for  rivers  of  greatly  diverse  size  and  discharge. 

133.  The  width-depth  ratio  and  the  sinuosity  p  were  found  to  be 
significantly  related  to  the  type  of  sediment  load  M  as  follows: 

_kf  =  225  W1  08  (80) 

p  =  0.94  M025  (81) 


The  results  showed  that  both  channel  shape  and  sinuosity  are  determined 
largely  by  the  type  of  sediment  load  moved  through  the  channels.  Schumm 
(1969)  stated  that 

.  .  .  with  data  acquired  in  the  field  and  from  gaging  station 

records  multiple  regression  equations  were  obtained  for  channel 
width  bf  ;  depth  dm  ;  width-depth  ratio  bt/dm  ;  and  meander 
wavelength  I  ;  as  a  function  of  channel  silt-clay  M  ;  and  mean 
annual  discharge  Qm  ;  or  mean  annual  flood  Qroa  (Schumm  1969).* 

After  reexamining  the  equations  he  developed  in  1968  (Equations  68-76),  Schumm 

(1969)  then  presented  the  following  revised  equations: 


bt  =  37 


,0,38 


Jr'S? 


(82) 


*  Schumm  (1969)  used  w  instead  of  bt  for  channel  width. 
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(83) 


<4  =  0.6  tfO.34qO.29 

4,  -0.09  w5-3^*2 
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(84) 


(85) 


(86) 


(87) 


(88) 


(89) 


(90) 


134.  Schumm  (1969)  then  pointed  out 

The  empirical  equations  presented  above  tell  several  things 
of  general  interest.  The  shape  (width-depth  ratio)  and  sinuosity 
of  an  alluvial  channel  are  primarily  determined  by  the  type  of 
sediment  load  or  the  ratio  of  bed-material  load  to  total  load 
moved  through  the  channels,  whereas  channel  width,  depth,  gradi¬ 
ent,  wavelength,  and  amplitude  are  significantly  influenced  by 
both  discharge  and  type  of  sediment  load.  The  equations  demon¬ 
strate  that  for  most  changes  of  hydrologic  regimen,  which  involve 
both  a  change  in  discharge  and  type  of  sediment  load,  many  aspects 
of  channel  morphology  will  change.  In  addition,  the  reaction  of  a 
channel  to  a  change  of  discharge  and  type  of  load  may  result  in 
changes  of  channel  dimensions  contrary  to  those  indicated  by  the 
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standard  regime  equations.  That  is,  it  is  conceivable  that  under 
certain  circumstances  with  a  decrease  of  discharge,  depth  will 
decrease  and  width  will  increase. 

Schumm  made  this  statement  because  Equations  82-90  were  empirical  (i.e.,  they 
yielded  directions  rather  than  magnitudes  of  change) . 

135.  Equations  82-85  and  88-89  show  that  channel  width,  depth,  and 
meander  wavelength  were  directly  related  to  discharge,  while  Equation  90  indi¬ 
cates  that  the  slope  is  inversely  proportional  to  discharge.  From  this  Schumm 
concluded: 


Q 


bf,  <4,  i 
S 


(91)* 


Mean  annual  discharge  or  mean  annual  flood  can  be  substituted  for  Q  . 

136.  Equations  82-83  and  88-90  demonstrate  that  channel  width,  meander 
wavelength,  and  slope  are  inversely  related  to  the  type  of  sediment  load  M  , 
while  Equations  81,  84,  and  85  show  a  direct  re1  onship  between  channel 
depth,  sinuosity,  and  M  . 

The  percentage  of  silt-clay  in  the  perimeter  of  a  channel 
reflects  the  nature  of  the  sediment  load  moving  through  that  chan¬ 
nel  expressed  as  the  percentage  of  bed-material  load  in  total 
load,  Qt  ,  but  for  a  given  channel  or  for  channels  of  equal  dis¬ 
charge,  M  probably  can  be  related  inversely  to  the  quantity  of 
bed-material  load  that  moves  through  the  channel;  therefore,  1/Qb 
can  be  substituted  for  M  in  Equations  80  through  90  if  discharge 
is  constant  (Schumm  1969).  Schumm  then  developed  the  following 
relation: 


Qb 


bt,  t,  S 

<4.  p 


(92)** 


The  width-depth  ratio  was  not  included  in  Equation  92,  although  this  variable 
was  highly  dependent  on  M  because  both  width  and  depth  appear  separately  in 
the  equation,  nevertheless,  the  relation  between  bt/dm  and  M  proved  useful 
for  interpreting  changes  of  channel  width  and  depth  in  subsequent  relations. 


*  Schumm  (1969)  used  Qw  instead  of  Q  for  discharge,  w  instead  of  b{ 
for  width,  and  L  instead  of  (  for  meander  wavelength. 

**  Like  Lane  (1955),  Schumm  (1969)  used  Qs  instead  of  Qb  . 
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137.  Schumm  (1969)  then  stated  the  following: 

To  discuss  in  more  detail  the  effects  of  changing  discharge 
and  sediment  load  on  channel  morphology,  a  plus  or  minus  exponent 
will  be  used  to  indicate  how,  with  an  increase  or  decrease  of 
discharge  or  bed-material  load,  the  various  aspects  of  channel 
morphology  will  change.  For  the  relatively  straightforward  cases 
of  an  increase  or  decrease  in  discharge  or  bed-materia1  load 
alone,  Equations  [93]  through  [96]  are  obtained. 
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An  increase  or  decrease  of  discharge  Q  alone  could  be 
caused  by  diversion  of  water  into  or  out  of  a  river  system.  An 
increase  of  Qb  can  result  from  increased  erosion  in  the  catch¬ 
ment  area  which  can  be  induced  by  deforestation  or  by  an  increase 
in  area  under  cultivation.  A  decrease  of  Qb  could  result  from 
improved  land  use  or  a  program  of  soil  conservation. 

An  increase  or  decrease  of  discharge  changes  the  dimensions 
of  the  channel  and  its  gradient,  but  an  increase  or  decrease  of 
bed-material  load  at  constant  mean  annual  discharge  changes  not 
only  channel  dimensions  but  also  the  shape  of  the  channel  (width- 
depth  ratio)  and  its  sinuosity. 

In  nature,  however,  it  would  be  rare  that  a  change  in 
discharge  or  sediment  load  could  occur  alone.  Generally,  any 
change  in  discharge  will  be  accompanied  by  a  change  in  the  type  of 
sediment  load  and  vice  versa.  This  is  because  these  two  variables 
are  dependent  on  other  factors  which  influence  both  runoff  and 
sediment  yield. 

When,  however,  the  effects  of  changing  water  discharge  and 
sediment  load  are  combined,  we  can  no  longer  assume  that  M  is 
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related  to  2/Qb  because  discharge  is  changing.  Therefore,  M 
can  be  used  only  as  an  index  of  the  ratio  of  bed-material  load  to 
total  sediment  load  or  of  the  percentage  of  bed-material  load  Qt 
in  the  equations  that  follow. 

Using  the  plus  or  minus  exponents  to  indicate  an  increase  or 
decrease  in  a  variable,  four  combinations  of  changing  discharge 
and  sediment  load  can  be  considered.  For  example,  if  both  dis¬ 
charge  and  the  percentage  of  bed-material  load  increase,  perhaps 
as  the  result  of  diversion  of  water  from  a  bedload  channel  into  a 
suspended- load  channel,  Equation  [97]  suggests  the  nature  of  the 
resulting  channel  changes 


Q'Qt 
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Equation  [97]  indicates  that,  with  an  increase  in  both  dis¬ 
charge  and  percentage  of  bed-material  load,  width,  meander  wave¬ 
length,  and  width- depth  ratio  should  increase  and  sinuosity  should 
decrease.  The  influences  of  increasing  discharge  and  percentage 
of  bed-material  load  on  channel  depth  and  gradient  are  in  opposite 
directions,  and  it  is  not  clear  in  what  manner  gradient  and  depth 
should  change.  However,  by  including  width-depth  ratio  in  Equa¬ 
tion  [97],  an  estimate  of  the  direction  of  change  of  depth  can  be 
obtained.  Width-depth  ratio  is  predominantly  influenced  by  type 
of  load  [Equations  80,  86,  and  87]  and,  therefore,  it  increases  in 
Equation  [97].  This  suggests  that  depth  will  remain  constant  or 
decrease  because  both  width  and  width-depth  ratios  increase. 

Channel  gradient  will  probably  increase  because  sinuosity 
decreases,  thereby  straightening  the  channel  and  increasing  its 
slope . 

When  both  Qt  and  Q  decrease,  as  could  result  from  dam  construction,  the 

reverse  of  Equation  97  pertains  as  follows 


r  ■» 
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When  as  common  in  nature,  the  changes  in  Q  and  Qt  are  in  opposite 
directions,  the  following  relations  are  obtained: 
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The  situation  expressed  in  Equation  [99]  could  result  from  a 
combination  of  controls.  For  example,  dam  construction  with 
impoundment  of  sediment  and  diversion  of  water  into  the  channel 
from  another  source.  The  situation  of  Equation  [97]  could  result 
from  increased  water  use  and  increased  land  use,  thereby,  decreas¬ 
ing  discharge  but  increasing  the  percentage  of  bed-material  load 
(Schumm  1969)  . 

138.  Schumm  (1969)  stated  that  the  above  interrelations  qualitatively 
suggest  that  channel  metamorphosis  could  occur  with  changes  of  discharge  and 
the  gradation  of  sediment  load.  The  magnitude  of  changes  could  be  calculated 
using  Equations  80  through  90  if  the  magnitudes  of  the  changes  in  discharge 
and  the  percentage  of  bed-material  load  in  the  total  sediment  load  are  known. 


Blench  (1970) 


139.  In  discussing  regime  theory,  Thomas  Blench  (1970)  pointed  out  that 
the  most  publicized  equations  to  evolve  from  India  had  been  those  developed  by 
Lacey  (1930,  1934).  In  recent  years,  the  Lacey  equations  had  "undergone 
repeated  assault  and  test  by  research  institutes  and  by  individuals"  (Blench 
1970).  A  number  of  modifications  had  been  made  to  the  original  equations,  but 
Blench  (1970')  believed  that  "the  only  alterations  justified  by  fact  are  in 
interpretation  and  in  the  unmeasurable  overall  silt  factor  by  resolution  into 
components  " 

140.  Blench  (1970)  emphasized  that  his  equations  were  restricted  to  the 
conditions  from  which  they  were  derived.  He  split  the  overall  silt  factor  in 
Lacey's  equations  into  parts  that  expressed  the  effects  of  bed  material,  cohe¬ 
sive  sides,  concentration  of  bed  material,  and  effective  viscosity  of  water. 
Additionally,  Blench  replaced  width  and  depth  by  the  wetted  perimeter  and 
hydraulic  radius  because  they  gave  slightly  better  correlations  and  simplified 
the  calculations. 
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The  findings  inherent  in  the  Lacey  original  analysis  are 
stated  herein,  along  with  three  independent  equations,  which  are 
applied  to  the  flow  of  a  steady  water- sediment  discharge  in  a 
canal  (with  the  following  constraints) : 

a.  The  Froude  number  in  terms  of  a  representative  depth  d 

(which  may  be  hydraulic  radius  or  depth  from  water  surface  to 
mean  sand  bed)  will  tend  to  acquire  a  value  dependent  on  the 
nature  of  the  water- sediment  complex.  Therefore,  a  practical 
bed  factor,  fb  ,  can  be  defined  by 


in  which  fb  depends  on  D  ,  CB  ,  ps  ,  pt  ,  v  and  factors 
defining  the  nature  of  the  bed-load  material  and  quantity  and 
nature  of  the  suspended  load.  In  Equation  [101],  Va  =  the 
mean  velocity  of  flow;  D  =  particle  diameter;  CB  =  the  bed¬ 
load  charge  defined  as  a  multiple  of  weight  per  second  per 
unit  of  weight  of  fluid  per  second;  ps  =  a  mean  particle  mass 
density;  pf  =  the  fluid  mass  density;  and  v  =  the  kinematic 
viscosity  of  fluid. 

b.  Cohesive  erodible-depositable  sides  appear  to  act  as  hydrau¬ 
lically  smooth,  and  equilibrium  can  be  achieved  within  a  range 
of  values  of  a  practical  side  factor  fs  defined  by 


= 


~F 


(102)* ** 


The  upper  limit  is  set  by  erodibility,  and  the  lower  by  depos- 
itability  from  suspension.  The  term  b  =  width  defined  by 
b  =  A/d  in  which  A  =  cross-sectional  area  of  flow.  This 
definition  ensures  that  Vbd  =  Q  in  which  Q  =  the  discharge. 

c.  The  Manning  equation  should  be  amended  to  read 

V0  =  n~1d3/t‘S1/2  ( 103 )  f 


in  which  n  =  an  absolute  rugosity  dependent  largely  oa  fb 
and  somewhat  on  fs  (Blench  1970) . 


*  Blench  (1970)  used  Fb  instead  of  fb  and  V  instead  of  F0  ,  C 
instead  of  CB ,  which  is  expressed  in  parts  per  hundred  thousand  by 
we i gh  t . 

**  Blench  (1970)  used  Fs  instead  of  fs  . 
t  Blench  (1970)  used  N  instead  of  n  . 
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141.  From  these  findings,  Blench  derived  the  following  set  of 


equations : 


(104) 


in  which 


=m 


(105) 
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He  reported  that  there  was  a  strong  parallel  between  determination  of 
Manning's  n  for  a  rigid  boundary  and  fb  and  fs  for  the  noncohesive 
mobile  sand  bed  and  erodible  cohesive  sides  of  a  canal.  The  evaluation  of 
Manning's  roughness  coefficient  and  the  silt  factors  of  Lacey,  although 
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guides,  are  given,  based  largely  on  experience.  Blench  (1970)  concluded  that 
his  method  offered  the  following  advantages: 

a.  'Equations  derived  from  data  of  self-adjusting  channels  of  the 
type  and  size  for  which  design  is  required. 

b.  "Provisions  of  the  three  necessary  and  sufficient  independent 
equations  for  finding  equilibrium  values  at  b  ,  d  ,  and 

5e  • 

c.  "Consolidation  of  the  unmeasurables  of  the  water  -  sediment 
complex  and  of  the  cohesive  erodible - depos i table  sides  into 
two  overall  parameters  fb  and  fs  that  can  be  calculated 
from  the  runn'"ng  conditions  of  existing  channels  in 
equilibrium  and  used  by  engineers  with  field  experience. 

d .  " S impl i c i ty . " 

Maddock  (1972) 

1m2.  Maddock  (1972)  introduced  two  equations  relating  regime  velocity 
to  hydraulic  roughness  and  transport  of  sediment  in  alluvial  channels  that 
move  appreciable  discharges  of  sediment, 
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-  tali  velocity  of  a  characteristic  particle  size  D  ,  fps 
f(D)  =  function  of  the  characteristic  particle  size  D 
C  =  concentration  of  the  total  sediment  discharge,  ppm 
The  values  for  f(D)  (Figure  5)  were  derived  from  flume  experiments  (Maddock 

19/2)  . 


f  (0) 


Figure  9.  Relation  between  mean  diameter  of  moving 
sediment  particles  and  1(D)  in  liquation  111 
derived  from  representative  flume  experiments 
(adapted  from  Maddock  1972) 

1.1  Dynamic  equilibrium  is  usually  the  result  of  opposing  tendencies. 

'•!  .ddo.-k  (  19/2)  concluded  that  in  alluvial  channels  the  opposing  tendencies  are 

:  •  u  t  he  tri'-tion  tact  or  and  the  bed  shear  stress  to  be  constant  or  for  the 

•  < , :  : .  r  v.,r  i  ance  of  the  friction  factor  and  the  bed  shear  stress  to  be  mini  - 

.  ;  ■  d  Tlie  I  i  m  i  i  of  this  re  1  •  t  i  on  is 


Its-:;,  u!  of 


where  Se  =  slope  of  the  energy  grade  line.  This  relation  had  been  common  to 
many  earlier  studies  of  the  hydraulics  of  alluvial  channels.  Maddock  (1972) 
stated 


The  relation  is  not  unique  because  there  is  a  coefficient, 

6  ,  associated  with  a  suite  (or  group)  of  bed  forms  having  similar 
energy-dissipation  characteristics.  If  velocity,  depth,  and  slope 
are  dependent  variables,  if  discharge  of  water  and  sediment  are 
independent  variables,  and  if  the  bed  is  formed  in  the  same  way-- 
i.e.,  by  deposition  in  a  channel  of  uniform  width- -then  ft  is  a 
constant.  If  the  response  of  the  dependent  variables  to  the  inde¬ 
pendent  variables  changes,  then  fi  is  not  a  constant  but  may  have 
a  systematic  variation  that  changes  with  the  independent 
variables . 

Coefficient  B  is  equal  to  4.7,  and  an  average  or  regime 
value  of  B  is  from  4  to  5.0;  some  obviously  unstable  streams 
have  values  as  low  as  2,  and  the  Marala-Ravi  Link  Canal  in 
Pakistan  has  a  value  of  about  7.  Streams  having  variable  dis¬ 
charges  of  sediment  and/or  water  that  are  associated  with  degra¬ 
dation  and  aggradation- -  as  at  the  heads  of  canals  and  in  the 
headwaters  of  drainage  basins- -usually  have  values  of  B  of  about 
3.5  to  4.  Conversely,  the  lower  ends  of  long  canals  of  uniform 
width  and  discharge  tend  to  have  high  values  of  B  .  As  the  area 
of  a  drainage  basin  increases,  the  distance  between  major  tribu¬ 
taries  usually  increases.  As  a  result,  the  variation  in  the  sedi¬ 
ment  discharge  tends  to  decline,  perturbations  are  dampened,  and 
the  value  of  B  increases  to  about  5  or  5.5.  The  value  of  B 
declines  slightly  as  the  sinuosity  increases.  In  Equation  [110], 
the  value  of  4.7  represents  a  straight  channel,  self-formed  by 
aggradation,  so  that  it  has  an  above  average  value  of  B  . 

144.  Equation  110  represents  a  dynamic  equilibrium  that  should  not  be 
regarded  as  a  precise,  but  rather  as  a  most  probable  relation  that  has  an 
expected  statistical  variation  about  a  mean  value.  Maddock  (1972)  reported 
that  there  was  a  maximum  concentration  of  sediment  for  a  given  sediment  size 
and  for  a  given  rate  of  expenditure  of  energy  per  unit  mass  (Equation  111) . 
Figure  5  indicates  that  at  least  in  the  smaller  particle  diameters,  mixtures 
having  a  median  particle  size  D  could  have  a  lower  value  of  f (D)  than  that 
of  equi granular  material.  For  sediment  having  median  particle  sizes  greater 
than  about  0.5  mm,  t(D)  is  proportional  to  the  square  root  of  the  charac¬ 
teristic  particle  diameter.  Maddock  (1972)  pointed  out  that  Equation  111, 
iike  Equation  110,  should  not.  he  "considered  deterministic.  Both  equations 
represent  tendencies  toward  minima." 

145.  Maddock  (1972  )  slated  that  the  relations  are  not.  affected  hy  the 
complexities  of  the  discharge  of  the  sediment  and  water  spectrum  as  long  as 

; he  bed  is  deformed  and  the  channel  is  sell -formed.  He  noted  that  Equa- 
i  on  111  was  not  applicable  to  shallow  flows  having  low  sediment 
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concentrations.  Equations  110  and  111  can  be  combined  to  yield  two  other 
equations 
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In  Figure  5,  f(D)  “  fD  for  particle  sizes  greater  than  0.5  mm;  assuming  that 
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in  Equations  106  and  107 
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166.  One  of  the  complications  in  using  the  fall  velocity  as  a  parameter 
is  the  effect  of  temperature  on  the  fall  velocity  of  spheres  and  the  effect  in 


st  reams 


Verv  strong  temperature  effects  have  been  observed  in  data  from  the 


Misso"ri  River 


1  lesser  effects  in  other  data.  Maddock  (1972)  commented: 


Equation  [1111  is  the  most  important  equation  in  this  paper, 
and  therefore,  its  limitation  must  he  understood- - that  is,  there 
must  tie  a  prescribed  amount  of  sediment  moving  on  the  bed  before 
there  ran  he  a  maximum  amount  of  sediment  transported  for  a  given 


87 


rate  of  energy  expenditure  per  unit  of  mass.  This  is  a  basic,  yet 
unrecognized,  limitation  of  the  regime  theory  and  is  also  a  limi¬ 
tation  on  the  application  of  the  equations  presented  herein. 

Blench  (1969)  had  stated  that  (US)  is  constant  for  a  given  particle  size  and 

sediment  concentration;  he  was  not  concerned  with  the  relation  among  (US)  , 

size  of  sediment,  and  concentration.  Inglis  (1967)  did  not  emphasize  the 

(US)  relations  in  any  way,  but  a  combination  of  his  equations  yielded 


(US)  oc  £) i/2c i/2v  1/2 


Lacey  (1946)  had  introduced  (US)  as  a  parameter,  in  which  (US)  -  D  ,  but 
did  not  recognize  the  importance  of  sediment  concentration  or  the  variability 
of  (US)  with  different  particle  sizes  and  temperatures. 

147.  If  /3  is  held  constant  in  Equation  114,  then  the  equation  is  a 
regime  equation.  In  an  earlier  paper,  Maddock  (1969)  stated  that  for  the 
equations  of  Blench  (1969),  Einstein  and  Barbarossa  (1952),  Engelund  (1966), 
Inglis  (1968),  and  Lacey  (1946),  this  would  also  hold  true;  however,  only 
Inglis  (1968)  had  used  the  same  form  in  his  equation  as  that  used  in  Equa¬ 
tion  114  (Maddock  1972). 

148.  Equation  114  was  derived  from  Equations  110  and  111  and  if  p  is 
constant,  Equation  114  can  be  expressed  in  the  form 
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Equation  117  confirms  Blench' s  statement  that  the  Froude  number  in  terms  of  a 
representative  depth  remains  the  same  for  all  channels  having  the  same  concen¬ 
trations  and  par  tide  -  size  distribution  (Maddock  1972).  Maddock  then  restated 
Equation  113,  in  the  following  form: 
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*  Maddock  (1972)  used  b  instead  of  bt  . 
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Commenting  on  Equation  118,  he  stated 

The  equation  is  descriptive  of  the  behavior  of  a  channel  at 
a  given  point.  The  fact  that  width  is  proportional  to  the  cube  of 
the  slope,  if  other  factors  are  constant,  was  shown  by  Maddock  in 
studies  of  the  Rio  Grande  near  Vinton,  TX.  The  equation  describes 
the  instability  of  a  channel,  as  represented  by  a  change  in  p  , 
when  a  change  in  discharge  is  not  accompanied  by  a  change  in  sedi¬ 
ment  concentration  or  size. 

Note  that  if  the  slope  remains  constant,  in  order  for  p  to 
remain  constant,  the  width  must  be  proportional  to  Q  unless 
there  is  a  compensating  change  in  the  size  or  concentration  of  the 
sediment.  Streams  have  a  tendency  to  widen  because  the  sediment 
concentration  frequently  fails  to  increase  sufficiently  during 
floods.  After  a  flood,  P  becomes  larger  for  the  normal  dis¬ 
charge  in  the  widened  channel  having  a  preflood  slope  and  concen¬ 
tration.  This  generally  leads  to  an  unstable  channel  that  tends 
toward  aggradation  or  bar  building.  Therefore,  in  most  instances 
the  reduction  of  width  is  an  effective  procedure  in  channel  recti¬ 
fication.  Using  Equation  [111],  it  can  be  determined  that  normal 
flow  in  the  rectified  channel  will  have  the  same  velocity  as  that 
in  the  wide  channel  regardless  of  the  increase  in  depth.  The 
narrow  channel  is  subject  to  attack- -too  low  a  value  of  p  --at 
high  discharges  and,  therefore,  must  be  protected  (Maddock  1972). 

149.  Maddock  (1972)  found  that  the  greatest  effect  of  a  change  in  sedi¬ 
ment  size  occurred  in  the  coarse  silt  and  fine  sand  sizes.  A  change  in  parti¬ 
cle  size  from  0.05  to  0.10  mm  decreases  the  sediment  concentration  about 
16  times  for  the  same  value  of  the  velocity-slope  product  (see  Figure  5).  In 
contrast,  a  change  in  particle  size  from  5  to  10  mm  decreases  the  sediment 
concentration  about  1.6  times  for  the  same  value  of  the  velocity- slope  prod¬ 
uct.  Secondly,  he  reached  some  unexpected  conclusions  regarding  the  effect  of 
temperature  on  Equations  110  and  111. 

In  Equation  [110],  the  only  changing  variable  is  v  ,  which 
is  affected  by  temperature.  Therefore,  the  presumption  would  be 
that  all  the  variations  in  v  are  taken  up  by  p  and  Se  .  In 
the  actual  experiment,  however,  a  temperature  reduction  from  86°  F 
to  35°  F  increased  velocity  about  9  percent  and  slope  54  percent. 
Assuming  that  the  fall  velocity  of  the  median  grain  size  is  pro¬ 
portional  to  a  sphere  of  the  same  size,  the  fall  velocities  for  a 
0.16-mm-diam  grain  (the  size  of  the  sediment  in  the  flume)  are 
about  2  cm  per  second  and  1.03  cm  per  second  at  temperatures  of 
86°  F  and  35°  F,  respectively,  or  a  difference  of  about  two  times. 
Introducing  the  observed  variations  into  Equation  [110],  the 
change  in  the  coefficient,  ft  is  1.05. 

The  small  change  in  the  coefficient,  P  ,  is  not  unexpected 
because  the  slope,  in  a  sense,  is  unconstrained  in  its  formation. 
However,  in  Equation  [111],  where  V0  is  constant,  Se  is 
related  to  f(D)  and  the  sediment  concentration.  The  unknown  is 
which  change  controls  the  adjustment- -a  change  in  ft  or  a  change 


in  f(D)  .  In  this  experiment,  most  of  the  change  was  in  f(D) 
(Maddock  1972). 

150.  From  the  temperature  effects  studies  on  the  Missouri  River  at 
Omaha  (US  Army  Engineer  Division,  Missouri  River  1969),  Maddock  (1972) 
observed 


If  during  falling  temperatures,  the  observed  increase  in 
sediment  concentration  exceeds  the  increase  in  sediment  concentra¬ 
tion  required  to  maintain  a  given  velocity,  then  there  is  an 
increase  in  velocity  of  the  stream.  If  the  observed  increase  is 
less  than  the  required  increase,  the  velocity  falls,  and  the  bed- 
form  factor  li  declines. 

Maddock  (1972)  then  concluded 

One  of  the  reasons  that  little  progress  has  been  made  in  under¬ 
standing  the  behavior  of  alluvial  channels  is  that  many  commonly 
accepted  assumptions  are  in  error,  such  as: 

a.  The  amount  of  sediment  of  a  specific  size  that  is  trans¬ 
ported  is  directly  related  to  the  amount  of  that  sedi¬ 
ment  present  in  the  bed. 

b.  For  every  unit  discharge  of  water  on  a  given  slope, 
there  is  some  fixed  discharge  of  sediment  of  a  given 
size . 

c.  The  way  an  experiment  is  conducted  in  a  flume  has  no 
influence  on  the  results. 

d.  A  recirculating  flume  is  equivalent  to  a  natural 
channel . 

e.  Specific  sequences  of  bed  forms  follow  only  specific 
changes  in  discharge. 

f.  The  channel  friction  factor  is  deterministic  and  can  be 
evaluated  by  the  equivalent  of  a  Moody  diagram,  in  which 
only  velocity,  depth,  slope,  and  sediment  size  are 
parameters- -  that  is,  the  value  of  any  one  of  the  parame¬ 
ters  may  be  determined  from  the  other  three. 

g.  It  is  necessary  to  divide  the  roughness  into  grain  and 
form  roughness . 

The  reason  for  the  above  assumptions  is  that  stream  behavior 
is  controlled  by  constraints,  either  natural  or  artificial,  on  the 
response  of  dependent  variables  to  changes  in  independent  vari¬ 
ables.  As  long  as  the  constraints  are  consistent,  the  responses 
are  consistent.  Thus,  the  response  under  the  constraint  of  con¬ 
stant  depth  or  constant  slope  is  different  from  the  response  when 
depth  and  slope  are  free  to  vary. 
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Simons  et  al.  ( 1975 ) 


151.  Simons  et  al .  (1975)  stated  that  understanding  river  system  behav¬ 
ior  depends  upon  the  competent  knowledge  of: 

a.  Geologic  factors,  including  soil  conditions. 

b.  Hydrologic  factors,  including  possible  changes  in  streamflows, 
runoff,  and  the  hydrologic  effects  of  changes  in  land  use. 

c.  Geometric  characteristics  of  the  stream,  including  the  proba¬ 
ble  geometric  alterations  that  will  be  activated  by  the 
changes  that  development  will  impose  on  the  channel. 

d.  Hydraulic  characteristics  such  as  depth,  slope,  and  velocity 
of  streams  and  the  changes  that  may  be  expected  in  these  char¬ 
acteristics  in  space  and  time. 

e.  Ecological/biological  changes  that  will  result  from  physical 
change  and  in  turn  will  induce  or  modify  physical  changes. 

They  listed  seven  major  factors  affecting  alluvial  stream  channel  forms: 

a.  Stream  discharge. 

b.  Sediment  load. 

c.  Longitudinal  slope. 

d.  Bank  and  bed  resistance  to  flow. 

e.  Vegetation. 

f.  Geology,  including  types  of  sediments. 

g.  Works  of  man. 

152.  Although  it  "exercises  its  greatest  influence  on  river  form,  char¬ 
acter,  and  resistance,"  bed  load  is  very  difficult  to  estimate  (Simons  et  al . 
1975).  According  to  Simons  et  al.  (1975),  three  major  variables  that  affect 
the  amount  of  bed  load  a  stream  can  carry  are: 

a.  Size  of  the  bed  material. 

b.  Slope  of  the  stream  or  average  velocity. 

c.  Nature  of  the  channel  (depth,  size,  shape,  and  roughness  of 
bed  and  banks ) . 

As  an  aid  in  evaluating  the  effect  of  these  variables  on  bed  load,  Lane  and 
Borland  (1951)  had  provided  the  following  criteria: 

a.  Smaller  concentrations  of  suspended  material  usually  imply 
higher  percentages  of  bed  load. 

b.  The  ratio  of  bed  load  to  suspended  load  is  usually  larger  for 
low  or  medium  stages  than  for  high  stages. 

c.  Streams  with  wide  shallow  channels  carry  a  higher  proportion 
of  sediment  as  bed  load  than  do  streams  with  deep  narrow 
channels . 
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d.  Streams  with  a  high  degree  of  turbulence  tend  to  have  smaller 
amounts  of  bed  load. 

e.  The  nature  of  the  source  of  sediments  influences  the  magnitude 
of  the  bed- load  correction;  that  is,  the  occurrence  of  large 
quantities  of  coarse  material  in  the  watershed,  so  located 
that  they  can  be  moved  easily  into  the  channel,  is  indicative 
of  higher  percentages  of  bed  load. 

153.  Simons  et  al .  (1975)  defined  "graded  streams"  (see  paragraph  2)  as 
"those  that  have  slopes,  such  that  their  energy  is  just  sufficient  to  trans¬ 
port  the  material  through  the  system  that  is  delivered  to  the  streams."  They 
stated  that  "...  this  concept  can  only  be  applied  as  an  average  condition 
extending  over  a  period  of  years"  (Simons  et  al .  1975). 

154.  Lane  (1957),  who  investigated  the  interrelation  of  slope,  dis¬ 
charge,  and  channel  pattern  in  meandering  and  braided  streams,  showed  that 
wher> 


SbQ1/4  <  0.0017 


for  sand  bed  channels,*  they  will  tend  toward  a  meandering  pattern.  Simi¬ 
larly,  when 

SbQ1M  >0.01 

the  stream  tends  toward  a  braided  pattern.  The  region  between  these  two  val¬ 
ues  can  be  considered  a  transitional  range  where  streams  are  classified  as 
intermediate.  This  intermediate  category  describes  many  of  the  rivers  in  the 
United  States  (Simons  et  al .  1975). 

155.  Lane  (1955)  had  studied  the  changes  in  river  morphology  in 
response  to  varying  discharge  and  sediment  load.  Additionally,  Leopold  and 
Maddock  (1953),  Schumm  (1971),  and  Santos - Cayudo  and  Simons  (1972)  investi¬ 
gated  channel  response  to  natural  and  imposed  changes.  These  studies  support 
the  following  general  relations: 

a.  Depth  of  flow  is  directly  proportional  to  discharge  and 
inversely  proportional  to  sediment  load. 

b.  Channel  width  is  directly  proportional  to  both  discharge  and 
sediment  load. 

c.  Channel  shape,  expressed  as  a  wi dth- to-depth  ratio,  is 
directly  related  to  sediment  load. 


*  Simons  et  al.  (1975)  used  S  instead  of  Sb 
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d.  Channel  slope  is  inversely  proportional  to  discharge  and 
directly  related  to  both  sediment  load  and  particle  size 

D50 

e.  Sinuosity  is  directly  proportional  to  valley  slope  and 
inversely  proportional  to  sediment  load. 

f .  Transport  of  bed  material  is  directly  related  to  stream  power 
and  concentration  of  fine  material  and  inversely  related  to 
the  fall  diameter  of  the  bed  material  D50  (Simons  et  al . 
1975). 

From  this,  evolved  Lane's  qualitative  relation 


a  Qb^SO 


The  above  relations  were  used  to  develop  Table  1,  in  which  Simons  et  al . 

(1975)  qualitatively  described  the  response  of  alluvial  channels  to  changes  in 
a  number  of  variables.* 


Maddock  (1976) 

156.  According  to  Maddock  (1976),  unique  relations  governing  the  inter¬ 
actions  of  pertinent  variables  describing  the  flow  of  water  in  alluvial  chan¬ 
nels  were  difficult  (if  not  impossible)  to  state.  He  did,  however,  recognize 
that  "values  of  certain  combinations  of  the  variables  are  constrained  within 
limits  that  are  defined  by  whether  the  variables  are  dependent  or  independent. 
Furthermore,  variables  that  are  independent  in  some  situations  are  dependent 
in  others"  (Maddock  1976).  Therefore,  in  the  short-term  field  situation, 
slope,  discharge,  suspended- sediment  load,  and  width  are  independent  vari¬ 
ables,  while  mean  velocity  and  mean  depth  are  dependent  variables.  However, 

In  a  long-term  field  situation,  discharge  and  suspended- sediment  load  are 
independent  variables,  whereas  slope,  width,  mean  velocity,  and  mean  depth  are 
dependent  variables  (Maddock  1976) . 

157.  Maddock  (1976)  stated  that  the  pertinent  variables  in  an  alluvial 
stream  flowing  with  uniform  depth  and  velocity  are  discharge,  sediment  concen¬ 
tration,  depth,  mean  velocity,  and  slope  of  the  energy  grade  line  (usually 
assumed  to  be  the  same  as  the  slope  of  the  water  surface),  the  kinematic  vis¬ 
cosity  of  the  fluid,  the  mass  density  of  the  fluid,  the  mass  density  of  the 


*  Simons  et  al.  (1975)  used  Qs  instead  of  Qb 
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Table  1 

ualitative  Response  of  Alluvial  Channels* 


bO  <u 

s-i  a 
C  0)  o 
o  e  >-1 

W  CO 


■  +  ■  + 


I  +  ■  + 


+  '  +  ■ 


+  '  + 


+1  +1  +1  +1 


+  ■  +  ' 


+  >  + 


>  +  •  + 


•  +  ■  + 


*  0Q 
U  -* 

<U  E 
>  l-i  t 
•H  O 

oi  u.  S 


3 

a)  o 
E  ^ 

Ic  + 

bO 

(U  c;-i 

06  O 


ai 

C  HI  H 
•rl  TJ  J3 

3  <0 

(U  cl  -H 
CO  >-i  + 

ECU 
tO  bO  > 

x.  tO 

U  £  ‘CJ 

o 


s  m  £  «  n  E 

r  t  r  t  t  t 

«  £  OQ  s  £  P3 


+  •  +  • 


> 

0d 

e 

o 

> 

0) 

co 

r— -4 

o 

Vj 

03 

4-1 

r— <  + 

•  4-> 

3-4 

U 

u-< 

cn 

■u 

3 

c 

c 

M 

O 

3 

a. 

O 

P 

Q 

f— 4 

C 

w 

0) 

cO 

c0 

o 

>— 1 

bO 

•H 

•r-4 

4-> 

•i-4 

XI 

3-4 

3-4 

3-4 

•«-4 

03  O) 

4J 

CO 

c0 

a> 

03 

03 

T3 

TJ 

C/3 

bO  U 

CO 

■  *— 1 

X 

u 

N 

U 

c0 

C0 

0 

cO  5-4 

•U 

V4 

o 

• 

c0 

•r-4 

i 

c0 

O 

X 

O 

u 

a  o 

03 

T3 

CO 

cn 

e 

U) 

TJ 

£ 

r— 1 

V) 

t— 4 

w 

03  X 

bO 

c 

> 

•i—4 

a; 

03 

c0 

•r-4 

03 

03 

•r4 

Q 

CQ 

CQ 

> 

CO 

> 

5 

*  Adapted  from  Simons  et  al .  (1975). 
**  M  =  meandering,  B  =  braided. 


sediment  particles,  the  median  size  of  the  sediment  particles,  the  standard 
deviation  of  the  sediment  sizes,  the  characteristic  fall  velocity  of  the  par¬ 
ticles,  and  the  acceleration  due  to  gravity. 

158.  Maddock  (1976)  realized  the  effect  of  bed  forms  were  very 
significant : 

There  are  any  number  of  observations  on  natural  streams  that, 
slope  being  constant,  a  range  of  velocities  may  be  observed  for 
the  same  mean  depth.  These  different  velocities  are  accompanied 
by  different  bed  forms  and  different  sediment  concentrations. 

He  introduced  a  resistance  equation: 


« .  ^  K  -  ipr 

(Wb  - 1)^1 

(gdSe)1/2  (dSe)1/z 

V2 

(119)* 


where  Vb  =  weight  of  bed  sediment.  He  stated  that  "...  the  near  constancy 

of  ft  under  unconstrained,  r r  self- formed,  conditions  is  the  only  reason  why 

presently  used  equations  have  any  engineering  application"  (Maddock  1976)  .  In 

addition,  Maddock  (1976)  pointed  out 

while  ft  is  essentially  constant  for  a  self-formed  channel,  it 
may  vary  between  channels.  For  example,  in  natural  channels  or 
canals  the  movement  of  sediment  is  not  uniform.  Where  perturba¬ 
tions  exist,  the  value  of  ft  for  a  bank- full  discharge  tends  to 
be  low.  However,  if  there  are  long  distances  between  tributaries 
of  major  streams  or  in  a  long  canal  of  uniform  section,  these 
perturbations  tend  to  be  dampened  and  the  value  of  ft  is  higher 
if  conditions  are  propitious  for  the  formation  of  a  flatter  slope 

159.  At  a  specific  location  on  a  natural  stream,  ft  is  a  variable 
depending  on  the  relative  discharge  of  water  and  sediment.  Under  these 
conditions,  high  values  of  ft  are  associated  with  channel  aggradation  and  low 
values  are  associated  with  channel  degradation  and  scour.  Values  of  V0  , 
d  ,  and  Se  can  be  used  to  compute  ft  and  predict  the  potential  degree  of 
channel  instability.  Analysis  of  experimental  data  led  to  the  following  con¬ 
centration  equation  (Maddock  1976)  : 


*  Maddock  (1976)  used  V  instead  of  V0  ,  D  instead  of  d  ,  s  instead  of 
f*'b  ,  d  instead  of  D  ,  o  instead  of  v  . 
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(120) 


1/4 

103 

V0Se  _  K(Wb  -  l)g1/2D 

K,  - 1)  gD 

HD)  f(D)d1/z 

Maddock  (1976)  plotted  several  sets  of  sample  data  to  demonstrate  the  excel¬ 
lent  correlation  of  Equation  120.  He  commented  that  the  relation  between  f(D) 
and  D  was 

not  that  of  particles  falling  in  still  water  but  represents  the 
summation  of  what  goes  in  the  stream,  and  there  is  no  reason  to 
suppose  a  still  water  fall  velocity  should  not  be  modified  by  the 
effect  of  turbulence,  particle  interference,  sorting  and  even 
variation  in  velocity  distributions  in  the  fluid.  Like  particle 
fall  velocities,  f(D)  responds  to  changes  in  temperature  (Mad- 
dock  1976). 

160.  Maddock  (1976)  stated  that 

K  is  a  dimensionless  parameter  that  has  a  normal  value  of  about 
60.  However,  when  the  mean  values  of  V0  or  d  are  used  in  an 
equation,  certain  velocity  distributions  are  implicit.  As  long  as 
there  are  a  large  number  of  particles  moving  over  the  bed  or 
within  the  fluid,  a  statistically  characteristic  velocity  distri¬ 
bution  may  be  achieved,  and  K  is  unimportant.  When  the  number 
of  moving  particles  declines,  either  from  a  decrease  in  concentra¬ 
tion  or  from  an  increase  in  size  at  the  same  concentration,  the 
moving  sediment  follows  those  locations  of  high  velocity  and  con¬ 
sequently,  sediment  moves  at  relatively  lower  mean  velocities. 

The  effect  of  this  is  to  reduce  K  to  as  little  as  30.  The 
skewed  velocity  distribution  may  be  caused  by  a  considerable  por¬ 
tion  of  the  bed  being  covered  by  stationary  sediment  particles,  in 
which  case  grain  roughness  controls.  However,  such  a  velocity 
distribution  can  also  result  from  bed  configuration.  Therefore 
with  few  sediment  particles  in  movement,  an  accurate  knowledge  of 
the  velocity  distribution  is  necessary  to  determine  the  rate  of 
sediment  transport.  Evaluation  of  K  is  dependent  on  estimates 
of  fall  velocity  and  sediment  size  which  are  hard  to  make. 

Consequently,  although  a  K  value  of  60  had  been  used  in  Maddock' s  paper, 
values  as  low  as  50  should  not  be  surprising.  He  arrived  at  two  conclusions 
(Maddock  1976)  : 

a.  Every  sediment  transport  equation  that  has  been  based  on  field  or 
laboratory  observations  has  an  area  of  application. 

b.  The  sorting  of  transported  sediment  is  an  important  method  of 
adjustment  to  changes  in  independent  variables. 
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Hev  (1978') 


161.  Researchers  have  attempted  "to  predict  the  hydraulic  geometry  of 
stable  alluvial  channels"  (Hey  1978) .  The  various  approaches  have  led  to  many 
different  equations,  and,  therefore,  different  solutions  to  stable  channel 
design.  Some  researchers  have  suggested  that  the  problem  is  indeterminate 
despite  the  observed  regularity  of  channel  shapes  and  patterns.  Richard  Pavid 
Hey  examined  the  reasons  for  the  diversity  of  opinion  and  presented  a  frame¬ 
work  that 

will  enable  a  general  deterministic  model  to  be  established.  Thi^ 
model,  based  on  an  understanding  of  the  processes  that  operate  in 
alluvial  channels,  indicates  that  it  is  possible  to  predict  the 
bank-full  hydraulic  geometry  of  unstable  a r,  well  as  stable  chan¬ 
nels  (Hey  1978)  . 

162.  The  early  empirical  work  by  Kennedy  (1895)  and  Lacey  (1930,  1934, 
1946)  resulted  in  a  set  of  equations  that  were  applicable  only  to  the  area  for 
which  they  were  derived.  Henderson  (1963)  and  Simons  and  Albertson  (1963) 
later  proved  that  bank  sediment  size  is  an  important  independent  variable  that 
should  have  been  included  in  those  equations. 

163.  Hey  (1978)  remarked  that  "the  uniqueness  of  these  empirical  regime 
equations  led  American  engineers  to  develop  a  more  rigorous  semitheoretical 
procedure  for  channel  design  process"  (e.g.,  Lane  1955).  This  approach  has 
been  referred  to  as  tractive  force  design  because  "every  particle  lining  (the 
channel)  perimeter  is  at  or  below  the  threshold  of  motion"  (Hey  1978)  .  In 
such  "fixed  bed  models,"  there  is  no  provision  for  bed  load. 

164.  Both  engineers  and  geologists  later  recognized  that  a  self- 
adjusting  channel  has  5  degrees  of  freedom  because  it  can  adjust  its  overall 
width,  depth,  velocity,  slope,  and  sinuosity  in  response  to  erosion  and  depo¬ 
sition  (Hey  1974,  Kellerhals  1967);  however,  Hey  (1978)  stated  that  "these 
five  variables  are  insufficient  to  uniquely  define  the  hydraulic  geometry  of 
alluvial  channels."  He  provided  three  reasons  to  support  his  case: 

a.  Width  and  average  depth  are  not  sufficient  to  describe  cross- 
sectional  geometry.  Hey  (1978)  stated  that  wetted  perimeter, 
hydraulic  radius,  and  maximum  flow  depth  d  provide  "a  unique 
definition  of  cross-sectional  shape"  and,  significantly,  they 
have  a  greater  rt  levance  than  width  and  depth. 

b.  Bed  forms  must  be  taken  into  account  and  their  size  and  shape 
must  be  predicted  to  obtain  a  determinant  solution  (Hey  1978, 
Yalin  1965)  . 
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c.  Sinuosity,  or  the  ratio  of  channel  length  to  valley  length, 

"does  not  uniquely  define  the  plan  geometry"  (Hey  1978) .  Many 
patterns  are  possible  for  a  gi/en  value  of  sinuosity.  Plan 
geometry  can,  however,  be  defined  "provided  arc  length  (i.e., 
the  channel  distance  between  successive  points)  is  also  speci¬ 
fied"  (Hey  1978)  . 

165.  Hey  (1978)  next  described  and  derived  the  following  process 
equat ions : 

a.  Continuity  equation 


V0  =  f(Q,  R,  P ) 


(121)* 


b.  Sediment  transport  equation 


5  =  f(P,  R,  dm 


D b, 


(122)** 


where 


dma  =  maximum  flow  depth 

Qsi  =  sediment  discharge  (input) 

Db  =  mean  bed  sediment  diameter 
ab  =  standard  deviation  of  bed  sediment  diameter 
Dr  =  mean  right -bank  sediment  diameter 

-  standard  deviation  of  right-bank  sediment  diameter 
Dl  =  mean  left-bank  sediment  diameter 
a i  =  standard  deviation  of  left -bank  sediment  diameter 
c.  Bed- form  equations 

•^a  =  ^(^>  &  >  ^tna  >  >  ab<  >  ° r  >  ^1> 

Aa  =  f(R,  P,  dma,  Sfr,  Db,  <7 b,  Dr,  or,  Dlt 


CT1.  Qsi) 


Qsi) 


(123) 


(124) 


*  Hey  (  1978)  used  Va  instead  of  ,  Ra  instead  of  R  ,  and  Pa 
instead  of  P. 

**  Hey  (1978)  used  a  instead  of  az  for  all  three  of  his  standard 
deviations . 
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where 


Aa  =  dune  wavelength 
Aa  =  dune  height 

d.  Bank  competence  equations 


P  -  f(R-£,  <4,  Sb,  Dlt  alt  Dfa,  ab ,  A,  A,  Qsi) 


(125)' 


<4  =  f(P£,  Pf,  5,  Dr,  a*.  Db,  a\t  A,  A,  Qsi) 


(126) 


where 

P£  =  bank- full  hydraulic  radius 
dm  =  maximum  bank- full  flow  depth 
A  =  bank- full  dune  wavelength 
A  =  bank- full  dune  height 
P£  =  bank- full  wetted  perimeter 
e.  Meander  equations 


P  =  Sv) 


(127) 


z  =  f(Pf,  <4) 


(128) 


P  f(P,  dma,  Sb,  Va,  Db,  crb  i  Dc,  CTr 


'!•  ul> 


P>  Z) 


(129) 


where  z  =  meander  arc  wavelength. 

166.  To  define  the  bank- full  hydraulic  geometry  of  alluvial  channels, 
Equations  121-129  must  be  solved  simultaneously  to  resolve  "the  interlinkage 
between  variables"  and  to  establish  "the  necessary  dependent  and  independent 
variables  to  give  a  unique  solution"  (Hey  1978).  The  following  relations 
describe  this  interlinkage : 


~  ney  (19/8;  used  S  instead  of  Sb  and  Db  instead  of  dm 
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V  =  f(Q,  Qsi,  Db,  al,  Dr,  o ?,  Dlt  ci\) 

Rf  =  f(Q ,  Qsi.  ^b.  Dr.  <^.  01.  Ol) 

Sb  =  f  (Q ,  Qsi ,  Db,  of,  £>r,  of,  Db,  oi) 

X  =  f  (Q ,  Qsi ,  Db,  of,  Dr>  4.  Oi.  «i) 

A  =  f(Q,  Qsi,  0b.  Ob,  0r.  of,  Di.  oi) 

P£  =  f(Q,  Qsi,  Db,  of,  0r,  of,  0L,  of) 

<4  =  f(Q,  Qsi,  Db>  of,  Dr,  Of ,  0i,  of) 

p  =  f ("Q ,  Qsi,  0b,  of.  0r.  of.  Db,  of,  Sv) 

z  -  f  CQ ,  Qsi,  0b.  of,  Dr ,  of,  dl,  of) 

167.  In  discussing  the  process  (Equations  121-129)  and  response  rela¬ 
tions  (paragraph  166),  Hey  (1978)  stated  that  because  the  process  equations 
define  the  interaction  between  variables  at  a  section  and  the  response  of  the 
system  to  the  imposed  conditions,  the  response  relations  have  been  developed 
for  unstable  as  well  as  stable  channels.  When  used  for  design  purposes,  an 
iterative  solution  would  have  to  be  adopted  to  ensure  channel  stability. 

After  the  channel  had  been  designed,  the  sediment  output  from  the  section 
would  be  calculated  and  compared  with  the  input  value.  Any  imbalance  between 
the  two  would  be  corrected  by  altering  the  bank  sediment  and  repeating  the 
process  until  the  two  values  were  in  equilibrium.  Although  bed  sediment  size 
could  also  be  adjusted  to  achieve  this  purpose,  this  would  be  unwise  because 
it  must  be  compatible  with  the  size  of  the  input  load.  Provided  that  dis¬ 
charge  and  sediment  input  are  invariant,  as  in  most  artificial  channels ,  this 
simple  model  can  be  used  for  design  purposes. 

168.  Commenting  further  on  the  response  relations,  and  citing  his 
earlier  work,  Hey  (1978)  continued 


Although  these  equations  (i.e.,  relations)  can  be  used  to  design  a 
stable  channel,  they  cannot  be  used  to  predict  channel  response 
consequent  upon  a  change  in  flow  regime  or  sediment  transport 
rates.  Under  these  conditions  upstream  and  downstream  feedback 
mechanisms  will  operate  and  all  the  variables  will  adjust  as  the 
channel  develops  towards  a  new  equilibrium.  For  example,  a  change 
in  flow  regime,  as  a  result  of  urbanization  or  river  regulation, 
will  affect  the  sediment  transport  characteristics  of  the  river 
and  the  size  of  the  bed  sediment.  As  the  channel  adjusts  its 
capacity  both  upstream  and  downstream  the  flood  regime  will  change 
and  this  will  modify  the  flow  variability  (Hey  1975a) .  By  widen¬ 
ing  the  perspective  to  incorporate  upstream  and  downstream 
responses,  casuality  can  no  longer  be  identified  and  all  the  vari¬ 
ables  are  interdependent.  Prediction  of  channel  response  in  these 
circumstances  is  not  only  depende’-*-  on  detailed  understanding  of 
equitf.  us  but  also  the  rates  and  mode  of  operation  of 
the  feedback  mechanisms.  Given  this  information,  it  will  then  be 
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possible  to  simulate  the  spatial  anil  temporal  response  of  the 
channel  to  a  given  stimulus  (Hey  1975b,  1976). 

169.  Hey  (1978)  believed  that  there  were  two  possible  methods  of  devel¬ 
oping  general  response  equations  to  define  the  bank- full  hydraulic  geometry  of 
stable  and  unstable  channels 

The  ideal  method  depends  on  establishing  precise  mathematical 
equations  for  each  adjustment  process  and  then  simultaneously 
solving  them  to  obtain  the  response  equations.  Unfortunately,  our 
lack  of  knowledge  regarding  these  processes  precludes  this  partic¬ 
ular  approach.  Alternatively,  they  (response  equations)  can  be 
established  by  submitting  field  data  to  multiple  regression  analy¬ 
sis  using  the  dependent  and  independent  variables  enumerated  in 
[paragraph  166].* 

Data  from  a  large  variety  of  rivers  would  have  to  be  used  to  ensure  the 
resulting  equations  have  general  application  (Hey  1978) . 


*  Paragraph  numbers  enclosed  in  brackets  within  quoted  material  refer  to  this 
report . 
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PART  III:  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Summary 


Kennedy  ( 18951 

170.  Kennedy  has  been  credited  as  being  the  first  person  to  develop 
equations  for  describing  canals  in  regime.  By  fitting  a  straight  line  through 
a  logarithmic  plot  of  depth  of  flow  versus  mean  velocity,  Kennedy  (1895) 
developed  the  following  relation  (Equation  2) : 


VQ  =  cdm 


(2  bis) 


where 

VQ  =  mean  velocity 
c  =  0.84 

d  =  average  depth,  ft 
m  =  0.64 

He  also  defined  U0  as  the  critical  velocity  or  the  velocity  at  which  silting 
is  just  prevented  (for  the  given  depth).  Due  to  the  consistency  of  Kennedy’s 
data,  he  originally  stated  that  bed  width  had  no  justification  for  being  cor¬ 
related  with  mean  velocity  and  depth  of  flow.  This  meant  that  either  a  deep 
narrow  channel  or  a  wide,  shallow  channel  could  be  designed  to  carry  the  same 
discharge.  Kennedy  (1904)  later  reassessed  his  research  and  developed  a 
"rough-rule"  relation  for  width- to -depth  ratio  (Graf  1971,  Lindley  1919). 
Lindlev  (1919) 

171.  Lindley  (1919)  retested  the  early  regime  theories  developed  by 
Kennedy  (1895)  and  Woods  (1917).  From  hydraulic  surveys  made  between  1915  and 
1917  of  the  Lower  Chenab  Canal,  Lindley  analyzed  widths  between  berms,  depths, 
and  gradients,  averaged  over  1-mile  reaches.  The  average  width-depth  relation 
could  be  expressed  as  (Equation  3): 

=  3  8  d1-61  (3  hi  s) 


where  bb  =  bed  width  developed  from  surface  width  under  the  assumption  that 
the  side  slopes  are  1V:2H  .  This  relation  was  valid  for  up  to  9  ft  of  depth, 
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with  no  indication  of  a  change  beyond  that  point.  Velocities  calculated  using 
Rutter's  relation,  with  N  =  0.9.25  ,  were  found  to  follow  the  exponential  law 
(Equation  4) : 


=  0.95  d0-57  bis) 

rather  than  Kennedy's  (1895)  Equation  2.  By  combining  the  two  expressions 
relating  width  and  mean  velocity  to  depth,  Lindley  (1919)  obtained  the  fol¬ 
lowing  relations  between  width  and  mean  velocity  (Equations  5  and  6) : 

for  Rutter  velocities, 


=  0.59  b^35 


(5  bis) 


for  Kennedy  velocities, 


=  0.50  fa,?'*0 


(6  bis) 


He  concluded  that  if  two  channels  carry  the  same  discharge,  the  one  with  the 
greater  "silt  charge"  (sediment  concentration)  must  be  wider  and  shallower  and 
have  steeper  banks.  Lindley  (1919)  had  advanced  the  quantitative  analysis  of 
regime  channels  by  defining  the  relation  between  velocity  and  width,  which 
added  to  the  already  established  relation  of  velocity  and  depth. 

Lacey  (1930) 

172.  Lacey  (1930),  suggesting  that  mean  velocity  is  a  function  of 
width,  depth,  and  a  "silt  factor"  f  ,  defined  the  silt  factor  as 
( Equation  7 ) : 


f  = 


(7  bis) 
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where  Va/Vok  =  ratio  of  the  regime  velocity  actually  obtained  with  a  given 
depth  to  the  velocity  for  the  same  depth  that  Kennedy's  (1895)  formula 
(Equation  2)  would  indicate  (Lacey  1930).  From  the  above  relation,  Lacey 
developed  his  first  general  equation  (Equation  8) : 


V0  =  1.17  JfR 


(8  bis) 


where  R  =  hydraulic  radius  .  To  test  the  relation  of  cross-section  area, 
silt  factor,  and  velocity,  Lacey  plotted  Af2  against  VQ  .  The  results  of 
this  plot  and  that  from  Equation  7  led  to  Equation  9,  which  relates  wetted 
perimeter  to  discharge: 


P  =  2.668  Q0-5 


(9  bis) 


where 

P  =  wetted  perimeter 
Q  -  discharge,  cfs 

This  equation  shows  that,  for  a  given  discharge,  the  wetted  perimeter  is  con¬ 
stant  and  independent  of  the  particle  size  transported,  with  particle  size 
merely  controlling  the  channel  shape.  Lacey  also  attempted  to  relate  the  silt 
factor  to  Manning’s  n  or  Kutter's  N  and  to  sediment  size.  The  general 
formulas  developed  by  Lacey  have  one  restriction.  For  velocities  less  than 
0.882  fps ,  the  resultant  wetted  perimeter  could  not  contain  the  cross- 
sectional  areas.  This  restriction  can  be  considered  either  a  limitation  of 
thi  formulas  or,  as  Lacey  (1930)  concluded,  an  argument  that  this  velocity 
must  be  one  of  physical  constants  of  water  because  Equations  7,  8,  and  9  do 
not  take  sediment  size  into  consideration. 

Lane  (1935) 

173.  Lane  (1935)  classified  the  formulas  previously  developed  for  sta¬ 
ble  channels  in  India  and  Egypt  into  two  types: 

a.  Those  giving  an  expression  of  velocity. 

b.  Those  giving  stable  channel  shapes. 

The  existing  regime  equations  had  been  developed  empirically,  and  Lane  con¬ 
tended  that  the  equations  omitted  important  factors  that  restricted  their  use 
to  the  channels  and  conditions  from  which  they  were  developed.  Lane  (1935) 
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contended  that  the  velocity  distribution,  mean  velocity,  and  quantity  of  mate¬ 
rial  in  motion  were  of  primary  importance.  He  concluded  that  the  capacity  of 
a  stream  to  transport  suspended  material  is  proportional  to  its  turbulence  and 
that  channels  carrying  bed  load  should  have  high  ratios  of  bed  width  to  depth. 
He  suggested  that  analysis  of  the  fundamental  interrelations  of  a  number  of 
factors  was  needed  to  understand  the  mechanisms  that  control  the  shape  of 
streams  in  alluvial  material: 

a.  Hydraulic  parameters  (bed  slope,  roughness,  hydraulic  radius, 
depth,  mean  velocity,  velocity  distribution,  and  temperature). 

b.  Shape  of  the  existing  channel  (width,  depth,  and  steepness  of 
side  slopes) . 

c.  Nature  of  material  transported  through  the  channel  (size, 
shape,  specific  gravity,  dispersion,  and  quantity,  as  well  as 
bank  and  subgrade  material) . 

d.  Other  factors  (channel  alignment,  uniformity  of  flow,  and 
"aging")  (Lane  1935). 

Lacev  (1936') 

174.  Lacey  (1936),  building  on  Lane's  (1935)  research,  correlated 
Lane's  factors  as  dimensioned  variables  using  the  Buckingham  n  theory.  This 
method  consisted  of  determining  dimensionless  arguments  from  all  the  known 
independent  variables  and  then  correlating  these  arguments.  A  poor  degree  of 
correlation  would  indicate  that  important  variables  had  been  overlooked.  From 
this  type  of  analysis,  Lacey  (1936)  developed  the  following  relation 
(Equation  17): 


P 

~R 


oc 


'RV0)ln 

v! 

,  —  J 

gR) 

(17  bis) 


where 

u  —  kinematic  viscosity  of  water 
g  =  acceleration  due  to  gravity 

This  is  similar  to  Lacey's  fundamental  shape  formula  (Equation  18): 


I  ■7-12V» 


(18  bis) 
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Lacey  (1936)  concluded  that  a  final  solution  to  the  problem  of  sediment  trans¬ 
port  would  be  found  by  investigations  using  dimensional  analysis. 

Shulits  (1936') 

175.  Shulits  (1936)  stated  that  the  extent  of  knowledge  of  fluvial 
morphology  had  generally  been  limited  to  qualitative  reasoning  and  some  site- 
specific  empirical  equations.  He  discussed  the  few  quantitative  principles 
that  had  been  developed  theoretically.  As  particles  move  downstream,  they  are 
reduced  in  size  by  abrasion  and  collision.  Based  on  Sternberg's  abrasion  law 
and  the  fact  that  the  slope  of  a  river  varies  directly  with  a  decrease  in 
bed-material  size,  Shulits  developed  an  exponential  equation  (Equation  19)  to 
approximate  the  slope  of  a  river: 


tfb  =  ^boe'C'  *  X 


(19  bis) 


where 

Wb  =  weight  of  bed  sediment 
I7b0  -  mean  weight  of  bed  sediment 
e  =  base  of  natural  logarithm 
c'  =  constant  of  proportionality 
0  =  coefficient  of  friction 
x  =  distance  traveled  downstream 

After  further  investigation  on  several  European  rivers,  Shulits  (1936)  con¬ 
cluded  with  reasonable  certainty  that  the  slope  of  a  river  could  be  expressed 
by  Equation  21: 


Z0  -  Z  =  K^l  -  e_Kx) 


(21  bis) 


where 

ZQ  =  mean  elevation 

Z  =  elevation  of  any  point  on  the  profile  above  a  given  datum 
KltK  =  constants 

He  modified  the  Schoklitsch  bed-load  formula  and  combined  it  with  the  abrasion 
and  slope  formulas  to  obtain  Equation  22: 
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KQ 

'  Kxb  ' 

e (4/3)ax 

gdrx/3 

(22  bis) 


where 

G  =  bed  load,  lb/sec 
a  =  coefficient  of  abrasion 

Therefore,  if  the  discharge  and  the  width  terms  increase  equally  downstream, 

'  r.ere  would  be  no  change  in  bed  load.  However,  if  the  discharge  term 
increases  more  rapidly  than  the  width  term,  there  will  be  an  increase  in  bed 
load  and  a  change  in  longitudinal  profile,  probably  degradation,  and  an 
increase  in  slope.  Should  the  width  term  increase  more  rapidly,  there  would 
be  a  decrease  in  bed  load  and  a  corresponding  change  in  channel  profile 
(Shulits  1936)  . 

Lacey  (1946) 

176.  Lacey  (1946)  attempted  to  advance  the  theory  of  silt  transported 
by  "stable"  channels  in  alluvium.  When  testing  the  earlier  regime  relations 
using  hydraulic  data,  he  discovered  two  curious  anomahv.^: 

a.  Despite  the  great  improbability  of  natural  alluvial  channels 
ever  attaining  regime,  alluvial  river  data  and  boulder  torrent 
data  (i.e.,  data  from  steeply  graded  streams)  confirmed  the 
regime  equations  in  many  respects. 

b.  Although  much  of  the  data  that  had  been  collected  over  a 
period  of  years  from  stable  channels  in  the  Punjab  displayed 
general  agreement,  certain  observations  from  the  channels 
thought  to  be  stable  proved  to  be  "discordant"  (i.e.,  con¬ 
flicting)  for  no  readily  ascertainable  reason. 

By  applying  the  principles  of  corresponding  speeds  and  dynamic  similarity  and 

by  introducing  the  terminal  velocities  of  silt  particles  into  his  physical 

theory,  Lacey  (1946)  succeeded  in  regrouping  the  fundamental  variables  of 

velocity,  depth,  and  water  surface  slope  to  obtain  two  new  parameters  from 

which  he  derived  Equation  25,  a  "normal  equation": 

(VS')  =  1.6  (d1/2S*r/3  (25  bls) 

where  ( VS *)  =  product  of  velocity  and  slope.  Lacey  stated  that  the  new  equa¬ 
tion  applies  to  all  active  channels  in  alluvium  and  is  irrespective  of  the 
precise  degree  of  scouring  or  silting,  or  variation  in  silt  charge.  All  of 
these  variables  are  implicit  in  the  actual  values  of  V  ,  d  ,  and  S  adopted 
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by  the  channel  when  associated  with  the  gradation  of  the  bed  material,  which 
must  be  sampled  and  measured  at  the  same  time  as  the  other  hydraulic  elements. 
In  addition,  Lacey  pointed  out  that,  for  every  stable  channel  transporting 
"regime  silt  charge,"  there  is,  for  a  fixed  diameter  D  of  the  bed  silt  par¬ 
ticles,  a  fixed  value  of  the  product  (VS*)  .  This  led  to  Lacey's  development 
of  Equation  26: 

D  -  0 . 73  (  VS’)  (26  bis) 

where  D  =  diameter,  millimetres,  and  Equation  27  (for  D  ,  inches): 

D  -  0.03 (VS’)  (27  bis) 

Lacey  (1936)  concluded  that  the  silt  charge  (sediment  concentration)  was 
clearly  of  importance. 

Inglis  (1947,1948) 

177.  Inglis,  like  those  who  preceded  him,  stated  that  alluvial  channels 
have  a  never-ceasing  tendency  to  attain  equilibrium,  a  condition  in  which  the 
load  is  carried  with  the  minimum  expenditure  of  energy  for  the  existing  condi¬ 
tions.  Small  adjustments  to  attain  equilibrium  are  brought  about  by  changes 
in  the  material  exposed  on  the  bed.  Such  changes  could  occur  either  by  depo¬ 
sition  "where  there  is  an  excess  charge"  (or  sediment  concentration) ,  or  by 
scour  "where  the  charge  is  in  defect"  (i.e.,  where  concentration  is  low) 
(Inglis  1947).  Inglis  (1947)  suggested  that  the  nearest  approach  to  equilib¬ 
rium  is  a  year-to-year  consistency  of  similar  discharge  conditions,  which  he 
called  an  "annual  stability."  He  developed  a  set  of  equations  to  describe 
river  equilibrium  and  then  compared  them  with  the  equations  developed  by 
Lacey.  In  Lacey's  formulas,  V  varies  from  P/R  only  in  the  coefficient. 

In  Inglis'  equations,  the  concentration  has  only  a  slight  effect  on  V  but 
has  a  very  great  effect  on  b/d  and,  hence,  P/R  ;  concentration  also  has  a 
marked  effect  on  slope  so  that  where  concentration  is  ignored,  wide  diver¬ 
gences  on  slope  are  found,  even  in  stable  channels.  For  the  purpose  of  sim¬ 
plification,  Inglis  (1947)  assumed  that  there  was  a  dominant  discharge  and  an 
associated  concentration  and  gradient  to  which  a  channel  returns  annually. 
Inglis  concluded  that  "for  any  discharge  less  than  the  dominant  discharge,  the 
slope  will  be  'in  defect1  (i.e.,  insufficient)  to  handle  the  dominant 
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concentration,  so  that  material  will  be  deposited,  and,  as  a  result,  there 
will  be  subdominant  charge  (.or  lower  concentration)  of  material  in  movement." 
On  the  other  hand,  when  the  discharge  exceeds  the  dominant  concentration,  the 
excess  energy  is  utilized  in  increasing  the  concentration  by  scouring  the  bed, 
eroding  the  banks,  and  developing  meanders  (Inglis  1947). 

Mackin  (1948) 

178.  Mackin  (1948)  defined  a  graded  (or  stable)  river  as  one  which, 
over  a  period  of  years,  delicately  readjusts  its  slope  "to  provide,  with 
available  discharge  and  the  prevailing  channel  characteristics,  just  the 
velocity  required  for  transportation  of  all  of  the  load  supplied  from  above." 
He  pointed  out  that  because  the  discharge,  channel  characteristics,  and  load 
do  not  vary  systematically  along  the  stream,  the  graded  profile  is  not  a 
simple  mathematical  curve.  Mackin  (1948)  also  believed  that  a  graded  stream 
responds  to  change  in  accordance  with  LeChatelier ’ s  general  law: 

If  a  stress  is  brought  to  bear  on  a  system  in  equilibrium  a  reac¬ 
tion  occurs,  displacing  the  equilibrium  in  a  direction  that  tends 

to  absorb  the  effect  of  the  stress. 

179.  Mackin  (1948)  stated  that  although  the  study  of  streams  as  agents 
of  transportation  can  be  approached  either  in  terms  of  relationship  between 
slope,  discharge,  channel  form,  and  the  size  of  grains  comprising  the  load,  or 
in  terms  of  energy  transformations,  the  two  approaches  should  not  be  combined. 
From  work  by  Gilbert  (1877)  and  Rubey  (1933),  Mackin  (1948)  emphasized  that 
the  share  of  total  energy  utilized  in  the  transportation  of  load  could  be  as 
small  as  3  percent.  Because  the  energy  dissipated  in  internal  and  external 
friction  is  overwhelmingly  greater  than  that  consumed  in  transportation,  the 
total  energy  determined  by  slope  and  discharge  remaining  the  same,  relatively 
slight  changes  in  the  channel  characteristics  cause  very  marked  changes  in  the 
transporting  power.  Many  researchers  had  developed  aquations  to  describe  the 
conditions  of  a  graded  stream,  but  Mackin  believed  that  this  approach  was 
inadequate  because  equations  are  transp usable .  As  set  up  in  an  equation,  load 
is  a  function  of  velocity.  However,  over  a  period  of  years,  the  load  supplied 
to  a  stream  is  actually  dependent  on  the  lithology,  relief,  vegetative  cover 
and  erosional  processes  in  operation  in  its  drainage  basin  and  in  the  graded 
stream.  The  most  efficient  cross-sectional  shape  for  transport  Taries  with 
the  slope,  the  amount  and  gradation  of  load,  the  resistance  of  the  banks  to 
erosion,  and  especially  with  the  proportions  of  the  total  load  carried  in 
suspension  and  moved  along  the  bed  (Mackin  1948).  Mackin  classified  changes 


in  regime  into  changes  in  downvalley  controls  and  changes  in  upvalley  con¬ 
trols.  In  general,  a  change  in  upvalley  controls  (load  and  discharge)  results 
in  a  change  in  slope,  whereas  a  change  in  downvalley  controls  (base  level) 
will  cause  a  change  in  the  level  of  the  longitudinal  profile,  but  not 
necessarily  in  the  slope.  Mackin  developed  many  general  statements  but 
avoided  manipulation  of  the  data  into  equation  form. 

Blench  (1951') 

180.  Blench  (1951)  stated  that  a  river  is  "in  regime"  if  its  mean  mea¬ 
surable  behavior  during  a  certain  time  interval  does  not  differ  significantly 
from  its  mean  measurable  behavior  during  comparable  times  before  or  after  the 
given  interval.  He  pointed  out  that  regime  theory  originated  on  the  great 
canal  systems  of  India  and  Pakistan  where  discharges  are  controlled  and  mean¬ 
ders  are  prevented;  however,  this  theory  is  not  generally  applicable  to  a 
"regime-type  river"  (Blench  1951).  Blench  found  the  exponents  of  Lacey's 
(1930)  regime  formula  (Equation  9)  to  be  consistent  with  the  laws  of  rigid 
boundary  hydraulics.  He  reasoned  that  as  width  goes  to  infinity,  depth 
approaches  the  hydraulic  radius  and  width  approaches  the  wetted  perimeter. 
Because  most  rivers  are  wide,  Blench  replaced  wetted  perimeter  with  width  b 
and  hydraulic  radius  R  with  depth  d  ,  and  leaving  Lacey's  exponents  the 
same,  arrived  at  Equation  42,  which  defines  a  bed  factor  fb  related  to  the 
nature  of  the  bed  load: 


=  4 


(42  bis) 


Blench' s  (1951)  second  regime  equation  (Equation  43)  defines  a  side  factor  fs 
that  measures  the  square  of  the  tractive  force  intensity  over  the  sides. 


43 
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-  4 


(43  bis) 


Blench  (1951)  developed  Equation  43  while  attempting  to  determine  the  reason 
for  the  universal  agreement  of  the  exponent  0.5  in  Lacey's  (1930)  formula 
(Equation  10) : 
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bt  =  2.67  Q0-5 


(10  bis) 


where  bt  =  Cop  width.  However,  he  contended  that  the  coefficient  2.67  should 
be  replaced  by  a  variable  dependent  upon  the  nature  of  the  bed  sediment.  He 
concluded  that  the  Lacey  formula  represents  an  averaging  of  the  relative 
importance  of  bed  to  sides  as  measured  by  the  ratio  fb/fs  meaning  that  all 
Lacey  channels  were  regime  channels  but  that  not  all  regime  channels  are  Lacey 
channels.  Equation  45,  developed  by  Blench,  became  his  third  regime  equation, 
thus  completing  the  theory: 


r,  VQb ^ 

1^1  -6  v 


(45  bis) 


where 

Sc  -  channel  slope  expressed  as  a  fraction 
C'  =  a  nondimens ional  constant 

Blench  emphasized  that  when  applying  regime  formulas  derived  on  canals  to 
rivers,  an  allowance  for  "disturbing  effects"  (i.e.,  turbulence)  and  variable 
discharges  must  be  made  "in  terms  of  engineering  common  sense  and  the  formulas 
regarded  with  suspicion  until  supported  by  results"  (Blench  1951)  . 

Leopold  and  Maddock  (1953) 

181.  Geomorphology  in  the  past  has  been  classically  treated  almost 
exclusively  in  a  qualitative  manner.  Leopold  and  Maddock  (1953)  believed  that 
the  qualitative  approach  had  been  constructive;  however,  they  proposed  that  a 
quantitative  analysis  should  be  attempted.  First,  they  investigated  changes 
in  velocity,  depth,  and  top  width  for  an  increase  in  discharge  at  a  particular 
cross  section  of  a  river.  Leopold  and  Maddock  (1953)  found  that  if  width, 
depth,  and  velocity  were  plotted  on  logarithmic  paper,  their  relations  to  dis¬ 
charge  could  be  approximated  by  straight  lines  and,  therefore,  described  by 
three  simple  power  functions  (Equations  50,  51,  and  52): 

b  =  ^Q"1’  (50  bis) 
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d  -  a2Qm’ 


(51  bis) 


WQ  =  a3Qm>  (52  bis) 

where  m1  ,  m2  ,  ra3  ,  ax  ,  a2  ,  a3  are  numerical  constants.  When  the 
width,  depth,  and  velocity  were  plotted  against  the  mean  annual  discharge  at 
various  cross  sections  along  the  river,  the  result  was  again  simple  power 
functions,  but  with  different  coefficients.  Leopold  and  Haddock  (1953)  also 
studied  the  effect  of  suspended  sediment  on  the  hydraulic  geometry.  Plots  of 
the  relation  of  suspended- sediment  load  to  discharge  were  characterized  by  a 
larger  scatter  of  points.  For  long  periods  of  record,  the  relation  can  be 
approximated  by  a  simple  power  function  (Equation  54)  (Leopold  and  Maddock 
1953)  : 

Qs  =  aQj  (54  bis) 


w}  ere 

Qs  -  suspended- sediment  load,  tons  per  day 
a,j  =  numerical  constants 

182.  They  concluded  that,  at  a  given  cross  section,  the  suspended- 
sediment  load  should  increase  more  rapidly  than  discharge  (Leopold  and  Maddock 
1953) .  Consideration  of  the  physical  characteristics  of  surface  runoff  and 
observations  of  bed  scour  led  to  the  tentative  conclusion  that  the  observed 
increase  in  sediment  concentration  resulted  primarily  from  erosion  of  the 
watershed  rather  than  from  scour  of  the  bed  of  the  main  stem  in  the  reach 
where  the  measurement  was  made.  For  each  station,  width  was  plotted  against 
suspended- sediment  load  on  a  semilogar ithmic  graph,  and  at  each  point  the  cor¬ 
responding  velocity  was  entered.  Isopleths  of  velocity  and  depth  were  then 
drawn.  This  graph  showed  that  there  was  an  increased  capacity  for  suspended- 
sediment  load  at  constant  discharge  by  decreasing  width  at  a  constant  velocity 
or  by  increasing  velocity  at  constant  width.  Some  authorities  had  contended 
that  despite  a  large  suspended  load  in  a  river,  only  the  bed  load  was  of  real 
significance  in  fluvial  morphology;  they  assumed  that  it  placed  the  greatest 
tax  on  the  energy  of  che  stream.  However,  Leopold  and  Maddock  (1953)  believed 
that  too  much  emphasis  had  been  placed  on  bed  load  and  that  the  role  of 
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suspended  load  in  determining  the  observed  characteristics  of  natural  streams 
had  been  underrated.  They  postulated  that  changes  in  channel  shape  occur  in 
response  to  a  change  in  suspended- sediment  load  brought  into  the  reach.  How¬ 
ever,  adequate  stream- gaging  data  were  not  available  in  a  sufficiently  close 
geographic  network,  or  over  a  long  enough  period  of  record,  to  allow 
differentiation  between  graded  and  nongraded  reaches  by  hydraulic  geometry 
alone,  though  such  differentiation  could  theoretically  be  possible. 

183.  The  effect  of  suspended- sediment  concentration  on  channel  rough¬ 
ness  has  been  discussed  by  many  researchers,  most  of  whom  concluded  that  an 
increase  in  concentration  resulted  in  lower  values  of  Manning's  n  .  Leopold 
and  Maddock  (1953)  found  a  change  in  velocity-depth  relations  that  coincided 
with  a  change  in  suspended  load-discharge  relations  at  a  discharge  well  below 
the  peak  on  the  Colorado  and  San  Juan  Rivers.  Such  local  adjustments  in 
velocity-depth  relations  result  primarily  from  changes  in  roughness  associated 
with  changes  in  sediment  load,  rather  than  from  changes  in  slope  of  the  water 
surface.  On  the  basis  of  their  work,  Leopold  and  Maddock  (1953)  agreed  with 
Mackin  (1948) ,  stating  that  the  wording  in  the  definition  of  regime  should  be 
changed  from  a  state  "in  which  slope  is  delicately  adjusted"  to  one  in  which 
slope  and  all  other  channel  characteristics  are  delicately  adjusted  (Leopold 
and  Maddock  1953). 

Lane 

184.  Lane  (1955)  was  a  proponent  of  the  energy  concept;  that  is,  ero¬ 
sion  of  the  bed  occurs  when  a  stream  is  capable  of  carrying  more  sediment  than 
is  supplied  to  it.  This  scour  increases  the  slope  of  the  tributary  streams, 
which  allow  them  to  introduce  more  material  to  the  main  stem.  He  stated  that 
a  stream  is  in  equilibrium  if,  "although  it  may  continually  fluctuate  between 
aggradation  and  degradation  over  a  long  period  of  years  in  terms  of  human  his¬ 
tory,  the  net  amount  of  change  is  not  sufficiently  large  to  be  detected  by 
quantitative  measurements"  (Lane  1955).  Lane  believed  that  most  alluvial 
streams  not  affected  by  the  works  of  man  were  in  equilibrium  for  most  practi¬ 
cal  engineering  purposes.  He  also  contended  that  if  equilibrium  existed  in  a 
reach  and  one  of  the  conditions  were  changed,  then,  over  time,  the  conditions 
of  equilibrium  would  change  throughout  the  entire  reach.  Because,  in  most 
cases,  data  are  insufficient  to  make  a  quantitative  analysis  of  equilibrium 
change,  Lane  developed  an  expression  relating  bed-material  load,  particle 
diameter,  discharge,  and  slope  of  the  stream  (Equation  55): 
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QdD  * 


(55  bis) 


where 

Qh  =  bed-material  load 
D  =  particle  diameter 
Sb  =  bed  slope 
Chien  (19551 

185.  The  formulas  proposed  by  Lacey  (1946)  involved  factors  based 
essentially  on  experience  alone,  making  use  of  the  theory  difficult  for  those 
who  lacked  experience.  Chien  (1955)  decided  to  evaluate  those  factors  quanti¬ 
tatively  in  terms  of  the  physical  conditions  imposed  on  the  channel  and  to 
show  the  possible  limitations  of  Lacey's  (1946)  regime  theory.  Chien  thought 
that  Lacey's  omission  of  sediment  load  as  an  explicit  variable  could  be 
explained  by  either  of  the  following  two  possibilities: 

a.  "The  silt  factors  .  .  .  actually  include  implicitly  the  sedi¬ 

ment  load  together  with  the  bed-material  size. 

b.  "Although  the  dimensions  and  the  slope  of  an  alluvial  channel 
must  depend  on  the  sediment  inflow,  a  certain  combination  of 
these  variables  may  depend  very  little  on  the  sediment  load" 
(Chien  1955). 

Chien' s  study  showed  that,  according  to  the  Einstein  bed-load  function,  the 
silt  factors  in  Lacey's  regime  theory  are  the  result  of  the  conditions  imposed 
on  the  channel  by  the  watershed.  In  addition,  he  found  that  for  conditions  in 
India  and  Pakistan,  the  silt  factor  fTO  depends  on  the  sediment  concentra¬ 
tion,  while  the  silt  factor  fRS  varies  primarily  with  the  bed-material  size. 
However,  for  other  conditions  fTO  also  depends  on  the  hydraulic 
characteristics  of  the  channel.  Therefore,  Chien  (1955)  concluded  that 
general  application  of  Lacey's  regime  theory  is  questionable. 

Schumm  (1960) 

186.  Schumm  (1960)  analyzed  the  effect  of  sediment  character  on  the 

shape  of  alluvial  streams.  He  pointed  out  ".  .  .to  relate  any  one  aspect  of 
stream  morphology  to  one  other  variable  .  .  was  an  oversimplification  of 

the  problem  and  the  reader  should  be  aware  of  the  other  important  factors  in 
stream  morphology.  Schumm  (1960)  described  the  sediment  comprising  the  perim¬ 
eter  of  each  channel  in  terms  of  its  weighted  mean  percent  silt-clay,  as  cal¬ 
culated  from  Equation  62: 
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(sc)(b£)  ♦  (sb)(2dm) 
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where 

M  =  weighted  mean  percent  silt-clay 
sc  =  percentage  of  silt-clay  found  in  the  channel  alluvium 
bf  =  bank- full  width 

sb  =  percentage  of  silt-clay  found  in  the  bank  alluvium 
The  shape  of  each  cross  section  was  expressed  as  a  width-depth  ratio  and  plot¬ 
ted  against  H  .  The  result  was  the  relation  expressed  in  Equation  63. 

=  255  M'108  (63  bis) 


This  relation  showed  that  channels  containing  little  silt-clay  are  relatively 
wide  and  shallow,  whereas  those  composed  predominantly  of  silt-clay  are  rela¬ 
tively  narrow  and  deep.  He  noted  that  the  shape  of  the  channels  seemed  to  be 
independent  of  the  discharge,  although  the  size  of  the  channel  was  related  to 
mean  discharge.  This  study  of  Schumm  (1960)  suggests  that  M  represents  the 
resistance  to  erosion  or  general  behavior  of  sediment  in  a  stream  channel  con¬ 
taining  only  small  amounts  of  gravel.  A  lack  of  correlation  between  channel 
depth  and  M  led  Schumm  to  conclude  that  perhaps  channel  width  is  more  sensi¬ 
tive  to  changes  in  M  than  to  changes  in  depth. 

Schumm  (1968) 

187.  Schumm  (1968)  conducted  a  study  of  the  Murrumbidgee  River  in  Aus¬ 
tralia  and  its  associated  paleochannels  in  the  Riverine  Plain.  The  percentage 
of  sediment  smaller  than  0.74  in  the  perimeter  of  a  channel  was  expressed  by 
Equation  62.  Schumm  (1968)  used  multiple  regression  analysis  to  obtain  equa¬ 
tions  for  the  following  relations,  whose  terms  are  transposable : 

a.  Bank-full  discharge,  percentage  of  silt  clay  in 
channel  perimeter,  and  meander  wavelength 
(Equation  65) : 


l 


438 


Qa“3 

W*7 


(65  bis) 


110 


where 

£  =  meander  wavelength,  ft 
Qs  =  bank- full  discharge,  cfs 

b.  Mean  annual  discharge,  mean  annual  flood,  silt-clay  percent¬ 
age,  and  meander  wavelength  (Equations  66  and  67): 


£  =  1,890 


M0-7* 


(66  bis) 
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(67  bis) 


where 

Qm  —  mean  annual  discharge,  cfs 

Qma  =  mean  annual  flood,  cfs 

c.  Discharge,  silt-clay  percentage,  and  channel  width 
(Equations  68  and  69): 

(68  bis) 


(69  bis) 
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d.  Discharge,  silt-clay  percentage,  and  channel  depth  (Equa¬ 
tions  70  and  71) : 


dm  =  0.59  M°-3*Q°-29 


(70  bis) 


dm  =  0.09 


0.42 

raa 


(71  bis) 


c.  Silt-clay  percentage  and  width-depth  ratio  (Equation  72): 
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=  106  W'°- 78  (72  bis) 


f.  Discharge,  silt-clay  percentage,  and  width-depth  ratio 
(Equations  73  and  74): 


*>£ 
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(74  bis) 


g.  Discharge,  channel  slope,  and  silt-clay  percentage 
(Equation  75) : 


sb 


=  59.5  M-°-38(?ro0'32 


(75  bis) 


where  Sb  -  bed  slope 

h.  Discharge,  channel  slope,  valley  slope,  and  silt-clay  percen¬ 
tage  (Equation  76) : 
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(76  bis) 


where  Sv  =  valley  slope 

i.  Sinuosity  and  silt-clay  percentage  (Equation  78): 


p  =  1.05  rt018 


(78  bis) 


where  p  =  sinuosity. 

Based  on  his  analysis,  Schumm  (1968)  concluded  that  the  morphology  of  a  stable 
alluvial  river  channel  reflects  the  hydrologic,  climatic,  and  geologic  charac¬ 
teristic  of  the  drainage  basin. 
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Schumm  (1969) 


188.  Schumm  (1969)  sought  to  demonstrate  that,  with  time,  river  regula¬ 
tion  could  cause  a  complete  metamorphosis  of  river  morphology.  He  found  no 
relation  between  channel  morphology  and  the  size  of  the  bed  and  bank  sedi¬ 
ments,  possibly  reflecting  the  small  range  in  bed-material  size  of  the 
channels  studied.  The  width-depth  ratio  and  the  sinuosity  were  found  to  be 
significantly  related  to  M  (Equations  80  and  81) . 


kf  =  225  AT108  (80  bis) 


p  =  0.94  M0-25 


(81  bis) 


He  revised  his  own  equations,  stating  that  channel  width,  depth,  and  meander 
wavelength  are  directly  related  to  discharge  (Equations  82-85  and  88-89, 
herein) : 
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c4  =  0.6rt°-8*Qm°-29 

dm  =  0.09  7f°-35QroV2 
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(82  bis) 


(83  bis) 


(84  bis) 


(85  bis) 


(88  bis) 
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flo.34 

t  =  1,890  v"a.  ( 8 9  bis) 


Equation  90,  on  the  other  hand,  indicates  that  the  slope  is  inversely  propor¬ 
tional  to  discharge: 


Sh  =  60  M-°-3BQm0-32  ( 9°  bis) 

In  addition,  Equations  82-83  and  88-90  demonstrate  that  channel  width,  meander 
wavelength,  and  slope  are  inversely  related  to  M  ,  while  Equations  81,  84, 
and  85  show  a  direct  relation  between  channel  depth,  sinuosity,  and  M  .  From 
these  ideas  suggested  by  the  empirical  equations,  Schumm  developed  10  propor¬ 
tions  to  qualitatively  describe  river  morphology.  These  proportions  related 
discharge,  bed-material  load,  total  load,  width,  depth,  meander  wavelength, 
slope,  width-depth  ratio,  and  sinuosity.  Schumm  (1969)  stated  that  his  inter¬ 
relations  qualitatively  suggest  that  channel  metamorphosis  could  occur  with 
changes  of  discharge  and  type  of  sediment. 

Blench  (1970) 

189.  Blench  (1970)  split  the  overall  silt  factor  in  Lacey's  equations 
into  parts  that  expressed  the  effects  of  bed  material,  cohesive  sides, 
concentration  of  bed  material,  and  effective  viscosity  of  water.  Blench  also 
replaced  width  and  depth  by  the  wetted  perimeter  and  hydraulic  radius  because 
they  gave  slightly  better  correlations  and  simplified  the  calculations.  From 
his  findings,  Blench  developed  a  set  of  equations  relating  width,  Froude 
number,  discharge,  slope,  acceleration  due  to  gravity,  and  kinematic  viscosity 
of  water.  He  emphasized  that  the  equations  were  restricted  to  the  conditions 
from  which  they  were  derived.  Blench  (1970)  concluded  that  there  was  a  strong 
parallel  between  determination  of  Mannings'  n  for  a  rigid  boundary  and  fb 
and  fs  for  the  noncohesive  mobile  sand  bed  and  erodible  cohesive  sides  of  a 
canal.  Like  Manning's  n  ,  fb  and  fs  are  evaluated  on  the  basis  of 
experience . 

Maddock  (1972) 

190.  Maddock  (1972)  introduced  two  equations  relating  regime  velocity 
to  hydraulic  roughness  and  to  transport  of  sediment  in  alluvial  channels  that 
move  appreciable  discharges  of  sediment  (Equations  110  and  111) : 


114 


PQS  =  10 -3f(D)C3/<’  (111  bis) 

where 

II  =  dimensionless  coefficient  related  to  a  characteristic  energy 
expenditure 

Y  =  unit  weight  of  water 
q  =  unit  discharge 

Ys=  submerged  unit  weight  of  sediment,  pcf 
p  =  mass  density  of  water 

v  =  fall  velocity  of  a  characteristic  particle  size  D 
f(D)  =  function  of  the  characteristic  particle  size  D 
C  =  concentration  of  the  total  sediment  discharge,  ppm 

191.  Maddock  (1972)  stated  neither  equation  can  be  considered  deter¬ 
ministic,  but  rather,  they  represent  tendencies  toward  minima.  He  noted  that 
Equation  111  was  not  applicable  to  shallow  flows  having  low  sediment  concen¬ 
trations.  Maddock  (1972)  stated  that  Equation  111  was  the  most  important 
equation  in  his  paper  but  that  its  limitations  should  be  understood.  There 
must  be  a  prescribed  amount  of  sediment  moving  on  the  bed  before  there  can  be 
a  maximum  amount  of  sediment  transported  for  a  given  rate  of  energy  expendi¬ 
ture  per  unit  mass;  Maddock  considered  this  to  be  a  basic,  yet  unrecognized, 
limitation  of  the  regime  theory.  In  Equation  110,  if  fl  is  a  constant,  then 
the  equation  is  a  regime  equation  and  includes  the  work  of  Einstein- 
Barbarossa,  Lacey,  Inglis,  Blench,  and  Engelund  (Maddock  1969). 

192.  Maddock  (1972)  combined  Equations  110  and  111  to  yield  two  other 
equations  (Equations  113  and  114) : 

(113  bis) 
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If  P  is  held  constant  in  Equation  114  and  the  equation  is  rewritten  as 
(Equation  117) : 
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(117  bis) 


Equation  117  confirms  Blench' s  statement  that  the  Froude  number  in  terms  of  a 
representative  depth  remains  the  same  for  all  channels  having  the  same  con¬ 
centration  and  particle  -  size  distribution  (Haddock  1972). 

193.  Maddock  (1972)  found  that  the  greatest  effect  of  a  change  in  the 
size  of  sediment  occurred  in  the  coarse  silt  and  fine  sand  size  and  that  tem¬ 
perature  affected  his  equations  to  a  degree  larger  than  expected  (Maddock 
1972)  . 

Simons  et  al .  (1975) 

194.  Simons  et  al.  (1975)  identified  seven  major  factors  affecting 
alluvial  stream  channel  forms: 

a.  Stream  discharge. 

b.  Sediment  load. 

c.  Longitudinal  slope. 

d.  Bank  and  bed  resistance  to  flow. 

e.  Vegetation. 

f.  Geology. 

g.  Works  of  man. 

Although  bed  load  "exercises  its  greatest  influence  on  river  form,  character, 
and  resistance,"  it  is  very  difficult  to  estimate.  The  authors  provide  the 
following  statements  (originally  developed  by  Lane  and  Borland  (1951))  about 
bed  load: 

a.  Smaller  concentrations  of  suspended  material  usually  imply 
higher  percentages  of  bed  load. 

b.  The  ratio  of  bed  load  to  suspended  load  is  usually  larger  for 
low  or  medium  stages  than  for  high  stages. 

c.  Streams  with  wide,  shallow  channels  carry  a  higher  proportion 
of  sediment  as  bed  load  than  do  streams  with  deep,  narrow 
channels . 
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d,  Streams  with  a  high  degree  of  turbulence  tend  to  have  smaller 
amounts  of  bed  load. 

e.  The  nature  of  the  source  of  sediments  influences  the  magnitude 
of  bed-load  correction. 

Simons  et  al .  (1975)  concluded  that  the  studies  by  Lane  (1955),  Leopold  and 
Maddock  (1953),  Schumm  (1971),  and  Santos - Cayudo  and  Simons  (1972)  supported 
the  following  general  relations: 

a.  Depth  of  flow  is  directly  proportional  to  discharge  and 
inversely  proportional  to  sediment  load. 

b.  Channel  width  is  directly  proportional  to  both  discharge  and 
sediment  discharge. 

c.  Channel  shape,  expressed  as  a  width- to- depth  ratio,  is 
directly  related  to  sediment  discharge. 

d.  Channel  slope  is  inversely  proportional  to  discharge  and 
directly  proportional  to  both  sediment  discharge  and  grain 
size  D5c  . 

e.  Sinuosity  is  directly  proportional  to  valley  slope  and 
inversely  proportional  to  sediment  load. 

f.  Transport  of  bed  material  is  directly  related  to  stream  power 
and  concentration  of  fine  material  and  inversely  related  to 
the  fall  diameter  of  the  bed  material  D50  . 

From  these  concepts,  Simons  et  al.  (1975)  developed  a  qualitative  analysis  of 

the  variables  affecting  alluvial  rivers  to  aid  in  the  assessment  of  reaction 

to  changes. 

Maddock  ('1976') 

195.  Maddock  (1976)  emphasized  that  unique  relations  governing  the 
interactions  of  pertinent  variables  describing  the  flow  of  water  in  alluvial 
channels  were  difficult  (if  not  impossible)  to  state.  His  realization  of  the 
important  influence  of  bed  forms  led  him  to  introduce  Equation  119,  a 
resistance  equation  defining  B  : 


R2  =  ^ 

11/2 
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He  stated  that  "...  the  near  constancy  of  B  under  unconstrained  or  self- 
iormed  conditions  is  the  only  reason  why  presently  used  equations  have  any 
engineering  application"  (Maddock  1976).  At  a  specific  location  on  a  natural 
stream,  B  is  a  variable  depending  on  the  relative  discharge  of  water  and 
sediment.  Under  these:  conditions,  high  values  of  B  are  associated  with 
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channel  aggradation,  and  low  values  are  associated  with  channel  degradation 
and  scour.  Values  of  Va  ,  d  ,  and  S  can  be  used  to  compute  ft  and  predict 
the  potential  degree  of  channel  instability.  After  Maddock  (1976)  analyzed 
his  experimental  data,  he  developed  Equation  120: 
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This  equation  shows  an  excellent  correlation  when  plotted  with  several  sets  of 
data.  From  his  research,  Maddock  (1976)  arrived  at  two  conclusions: 

a.  Every  sediment  transport  equation  based  on  field  or  laboratory 
observation  has  an  area  of  application. 

b.  The  sorting  of  transported  sediment  is  an  important  method  of 
adjustment  to  changes  in  independent  variables. 

He y  (  1978) 

196.  Prior  to  1978,  a  number  of  researchers  attempted  "to  predict  the 
geometry  of  stable  alluvial  channels"  (Hey  1978).  The  various  approaches  led 
to  many  different  equations  and  different  solutions  to  stable  channel  design. 
Some  researchers  suggested  that  the  problem  is  indeterminate  despite  the 
observed  regularity  of  channel  shapes  and  patterns.  Hey  (1978)  pointed  out 
that  uniqueness  of  the  Indian  empirical  regime  equations  resulted  in  American 
engineers'  developing  a  semitheoretical  procedure  for  channel  design.  How¬ 
ever,  the  American  theory  required  the  bed  to  be  on  the  threshold  of  motion; 
therefore,  no  bed  load  is  present.  Contrary  to  previous  theories,  Hey  con¬ 
tended  that  a  river  possessed  nine  degrees  of  freedom;  he,  therefore,  derived 
a  set  of  nine  process  equations  and  nine  response  relations  (which  show  the 
interrelation  between  a  number  of  variables) . 

197.  To  define  the  bank- full  hydraulic  geometry  of  alluvial  channels, 
the  process  equations  must  be  solved  simultaneously  to  resolve  "the  interlink¬ 
age  between  variables  and  establish  the  necessary  dependent  and  independent 
variables  to  give  a  unique  solution"  (Hey  1978).  Because  the  process  equa¬ 
tions  define  the  interaction  between  variables  at  a  section  and  the  response 
of  the  system  to  the  imposed  conditions,  the  response  relations  have  been 
developed  for  unstable  as  well  as  stable  channels.  These  equations  and  rela¬ 
tions  can  be  used  to  design  stable  channel.*;  however,  because  of  upstream  and 
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downstream  feedback  mechanisms,  they  cannot  be  used  to  predict  channel 
response  dependent  upon  a  change  in  flow  regime  or  sediment  transport  rates. 
Hey  (1978)  believed  that  there  were  two  possible  methods  of  developing  general 
response  equations  to  define  the  bank-full  hydraulic  geometry  of  stable  and 
unstable  channels.  The  ideal  method  would  be  to  establish  precise  mathemati¬ 
cal  equations  for  each  adjustment  process  and  then  solve  them  simultaneously. 
However,  at  the  current  state  of  knowledge,  this  is  nearly  impossible.  There¬ 
fore,  multiple  regression  analysis  of  field  data  from  a  large  variety  of 
rivers  is  a  practical  solution  to  developing  general  equations  (Hey  1978). 

Conclusions 


198,  A  number  of  researchers  who  have  attempted  to  quantitatively 
describe  regime  channels  have  found  that  the  lack  of  data,  variation  in  dis¬ 
charge,  and  heterogeneous  nature  of  the  banks  make  analysis  of  large  alluvial 
rivers  difficult.  Therefore,  most  studies  have  resulted  in  empirical  equa¬ 
tions  derived  from  a  limited  amount  of  data  obtained  on  small  streams  or 
canals  where  the  variables  affecting  regime  characteristics  are  simplified. 

As  early  as  1935,  Lane  made  a  statement  that  demonstrated  the  limitation  of 
empirical  equations  in  describing  regime: 

Although  these  formulas,  no  doubt,  provide  workable  relations  for 
the  conditions  for  which  they  were  developed,  these  conditions 
have  not  been  delineated  sufficiently  to  enable  them  to  be  applied 
elsewhere.  In  general,  also,  they  were  developed  empirically  from 
a  very  limited  range  of  conditions,  and  in  most  cases  they  omit 
important  factors  from  consideration  (Lane  1935) . 

Later,  Blench  (1951)  cautioned  the  engineer  about  applying  regime  formulas  for 

conditions  that  vary  frcm  those  for  which  they  were  derived.  When  applying 

regime  formulas  derived  on  canals  to  rivers,  an  allowance  for  disturbing 

effects  and  variable  discharges  must  be  made  "in  terms  of  engineering  common 

sense  and  the  formulas  regarded  with  suspicion  until  supported  by  results" 

(Blench  1951). 

199.  Due  to  the  limited  use  of  empirical  equations,  several 
researchers,  including  Lane  (1935,  1955,  1957);  Mackin  (1948);  Schumm  (1960, 
1963,  1968);  and  Simons  et  al .  (1975)  developed  qualitative  equations;  how¬ 
ever,  equations  are  transposable ,  meaning  that  an  equation  relating  sediment 
load  and  velocity  can  be  restated  to  show  that  load  is  dependent  on  velocity. 
The  sediment  load  actually  is  a  function  of  the  lithology,  relief,  vegetative 
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cover,  and  erosional  processes,  so  that,  over  a  period  of  years,  the  velocity 
is  adjusted  to  carry  the  load  supplied  to  the  stream.  This  relation  is  not 
transposable .  The  qualitative  relations  suggest  a  direction  of  response  of 
the  river  due  to  a  change  in  one  or  more  of  the  variables.  These  suggested 
responses  should  be  treated  only  as  guides . 

2u0.  The  studies  presented  in  Part  II  show  that  researchers  cannot 
agree  on  the  exact  nature  of  interrelations  between  the  variables  affecting 
the  regime  of  alluvial  rivers;  however,  most  of  them  do  agree  that  the  follow 
ing  variables  are  involved  and  that  they  are  interrelated: 

a.  Geology- -sources  of  erodible  materials. 

b.  Hydrology- - including  precipitation,  streamflow,  and  runoff. 

c.  Geometry- - including  channel  cross-sectional  shape,  alignment, 
sinuosity,  bed  slope,  and  aging. 

d.  Hydraulics- -including  water  and  energy  slopes,  depth  of  flow, 
roughness,  velocity  distribution,  and  temperature. 

e.  Nature  of  material  transported- -  including  gradation,  mineral¬ 
ogy,  shape,  specific  gravity,  and  quantity  of  bank  and  bed 
material  available. 

201.  In  addition,  changes  in  land  use  and  man's  activities  (e.g.,  con¬ 
struction)  alter  the  natural  rate  of  runoff  and  soil  loss,  which,  in  turn, 
affects  the  regime  of  alluvial  rivers.  Even  though  investigators  agree  that 
changes  in  suspended-sediment  load  influence  the  regime  of  alluvial  rivers, 
the  exact  role  that  this  parameter  plays  remains  to  be  determined.  More 
research  is  needed  to  study  this  influence  as  well  as  the  interrelations  of 
suspended-sediment  loads  and  the  other  variables  affecting  the  regime  of  allu 
vial  rivers. 


Recommendations 


202.  There  are  no  generally  accepted  quantitative  relations  describing 
the  response  of  river  regime  to  changes  in  sediment  load.  The  future  develop 
ment  of  quantitative  relations  will  most  probably  be  an  approach  suggested  by 
Hey  (1978) : 

The  ideal  method  depends  on  establishing  precise  mathematical 
equations  for  each  adjustment  process  and  then  simultaneously 
solving  them  to  obtain  the  response  equations.  Unfortunately,  our 
lack  of  knowledge  regarding  these  processes  precludes  this  partic¬ 
ular  approach.  Alternately,  they  can  be  established  by  submitting 
fiel^  data  to  multiple  regression  analysis.  .  .  . 


120 


At  present,  discharge,  sediment,  and  channel  geometry  data  are  widely  avail¬ 
able  for  many  streams,  including  the  Mississippi  River  (Keown,  Dardeau,  and 
Causey  1981).  In  addition,  improvements  for  flood  control,  navigation,  water 
supply,  erosion  control,  and  river  training  are  well  documented.  Hydraulic 
and  meteorologic  data  are  now  routinely  collected  and  transmitted  via  geosta¬ 
tionary  satellites  to  ground  stations  for  rapid  processing.  Further,  Land- 
sat  4  provides  continuous  sun- synchronous  coverage  of  the  Earth  every  16  days. 
Band  7  of  Landsat  4's  Multispectral  Scanner,  which  is  sensitive  to  radiation 
in  the  near- inf rared  portion  of  the  electromagnetic  (EM)  spectrum  (0.8  to 
1.1  pm),  is  useful  for  mapping  water  bodies  and  delineating  land-water  inter¬ 
faces.  Landsat  4  also  has  a  seven-band  Thematic  Mapper  (TM) .  Band  1  of  the 
spacecraft's  TM,  covering  part  of  the  visible  EM  spectrum  (0.45  to  0.52  pm), 
is  designed  for  water  body  penetration,  making  it  useful  for  estimating 
suspended- sediment  surface  concentrations.  The  French  satellite  SPOT  also  has 
potential  for  sediment  studies  in  either  the  multispectral  or  panchromatic 
mode.  In  the  multispectral  mode,  the  reflected  infrared  band  (0.79  to 
0.89  pm)  is  most  useful  for  mapping  land-water  interfaces,  while  the  green 
band  (0.50  to  0.59  pm)  shows  the  most  promise  for  water  body  penetration. 

203.  Because  all  of  these  information  sources  are  potentially  adaptable 
to  stream  regime  analysis  and  because  automated  processing  techniques  can 
provide  the  means  to  handle  the  large  quantities  of  available  data,  the  fol¬ 
lowing  process  is  recommended: 

a.  A  multiple  regression  approach  in  conjunction  with  a 
qualitative  analysis  be  used  to  assess  the  interrelations  of 
the  identified  parameters  that  influence  stream  regime. 

b.  The  results  of  this  sensitivity  analysis,  in  turn,  be  used  to 
evaluate  the  impact  of  changing  sediment  loads  on  a  basin's 
regime  or  on  the  character  of  a  specific  stream. 
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APPENDIX  A:  NOTATION 


a>  al .  a2<  a3 

A,  A 
b 

hb 

bf 

fat 

c 

c ' 

C 


C1 


CB 

d,  d'  ,  d" 

do 

dm 

dma 

D 

Db 

De 

Di 


ft 

fs 

f’  fRS’  fvR 
f(  ) 

F 


Coeff icients 
Cross-sectional  area 

Width  of  a  rectangular  channel  having  depth  d,  ft 

Bed  width  developed  from  water  surface  width  under  the 
assumption  that  side  slopes  are  IV: 2H 

Bank- full  width 

Top  width 

Coefficient 

Constant  of  proportionality 
Concentration  or  charge,  ppm 
Constant 

Bed- load  concentration  (or  charge) 

Average  depths  of  flow,  ft 
Constant  depth 

Maximum  bank- full  flow  depth 

Maximum  flow  depth 

Particle  diameter 

Mean  diameter  of  bed  material,  mm 

Effective  diameter  of  bed  material,  mm 

Mean  left-bank  sediment  diameter 

Mean  right-bank  sediment  diameter 

Base  of  natural  logarithm 

Bed  factor 

Side  factor 

Silt  factors 

Function  of  (  ) 

Froude  number 

A1 


s 


G 

H 

Hb 

P  e 

Hn 

Hs 

fit 


J 

K ,  Kx 
i 

L,  L' ,  L" 


m,  mi,  mz  >  m3  >  m ' 
M 

n 

N 

P 

P,  P' 

Pi 

q 

qT 

Q 

Qa 

Qb 

Qm 


Acceleration  due  to  gravity 
Bed  load,  lb/sec 

Total  energy  available  after  precipitation  (=  elevation 
above  mean  sea  level) 

Energy  dissipated  as  heat  through  friction  of  flow 

Energy  expended  in  erosion  of  rocks,  soil,  etc. 

Energy  extracted  by  artificial  means  (e.g.,  turbines) 

Energy  used  to  transport  sediment  and  debris  loads 

Energy  lost  by  turbulence  (e.g.,  rapids,  stream 
expansions,  contractions,  and  structures) 

Energy  gradient 

Integer 

Constants 

Meander  wavelength 

Lengths  of  flow  or  longitude  scalar  distance,  ft 
Constants 

Percent  silt-clay  fraction  in  perimeter  of  channel 
Manning's  n 

Rutter's  rugosity  coefficient 
Sinuosity 
Wetted  perimeter 
Bank-full  wetted  perimeter 
Discharge  per  unit  width 

Total  sediment  transport  rate  per  unit  width 

Discharge,  cfs 

Annual  discharge 

Bed-material  load 

Mean  annual  discharge,  cfs 


A2 


Qma 

<?s 

Qsi 

Qt 

Qo 

R 

Rt 

S,  S' ,  S" 
Sb 

Sc 

St 

Sc 

Se 

So 

Sv 

Sw 

s* 

T,  T'  ,  T" 
v,  v' ,  v" 
V,  V ' ,  V" 
(VS) 

Vo 


v0/vo k 


^bo 


Mean  annual  flood,  cfs 
Suspended- sediment  load 
Sediment  discharge  (input) 

Percentage  of  total  sediment  load  that  is  bed  load 

Bank- full  discharge,  cfs 

Hydraulic  radius 

Bank- full  hydraulic  radius 

Slopes 

Percentage  of  silt-clay  in  bank  alluvium 
Percentage  of  silt-clay  in  channel  alluvium 
Bed  slope 

Channel  slope  expressed  as  a  fraction 
Slope  of  the  energy  grade  line 
Slope  at  the  head  of  the  reach 
Valley  slope 

Slope  of  the  water  surface 
Slope,  parts  per  thousand 
Scalar  times 
Terminal  velocities 
Velocities,  fps 

Lacey's  (1946)  product  of  velocity  times  slope 

Mean  velocity  or  regime  velocity,  fps 

Ratio  of  the  regime  velocity  actually  obtained  with  a 
given  depth  to  the  velocity  for  the  same  depth  that 
Kennedy's  formula  (Equation  2)  would  predict 

Weight  of  bed  sediment 

Mean  weight  of  bed  sediment 

Distance  traveled  downstream 


A3 


N.Q  CVJ*-*  CM  L| 


X 

z 

Z 

Zo 

a 

al >  a3 .  q5>  a7 

B 

Y 

Ys 

A 

Aa 

X 

*a 

o 

p 

Pf 

Ps 

a 


</> 


Integer 

Meander  arc  wavelength 

Elevation  of  any  point  on  the  profile  above  a  given  datum 
Mean  elevation 
Coefficient  of  abrasion 
Coefficients 

Dimensionless  coefficient  related  to  a  characteristic 
energy  expenditure 

Wt ight  of  water 

Submerged  unit  weight  of  sediment,  pounds  per  cubic  bolt 
Bank- full  dune  height 
Dune  height 

Bank- full  dune  wavelength 

Dune  wavelength 

Kinematic  viscosity  of  water 

Density  of  water 

Fluid  mass  density 

Mean  particle  mass  density 

Constant 

Standard  deviation  of  bed  sediment  diameter 
Standard  deviation  of  left-bank  sediment  diameter 
Standard  deviation  of  right-bank  sediment  diameter 
Coefficient  of  friction 


A4 


